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This manual describes the procedures and 
data required to run the FOC version of the 
LiACIE Performance Predictor produced 
under Contract NAS9-14547. NASA/JSC 
Contract Technical Monitor is I. D. Browne, 
Earth Observations Division. 
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l.O GENERAL INFORMATION 
1. I SUMMARY 


The LACIE Performance Predictor (LPP) is a computer simulation 
of the LACIE process for predicting worldwide wheat prcduetion. The 
eirhulation provides for the introduction of various errors into the system 
and provides estimates based on these errors, thus allowing the user to 
determine the impact of selected error sources. 

The FOC LPP simulates the acquisition of the sample segment data 
by the LANDSAT Satellite (DAPTS), the classification of the agricultural 
area within the sample segment (CAMS), the estimation of the wheat yield 
(YES), and the production estimation and aggregation (CAS). 

These elements include data acquisition characteristics, environmental 
conditions, the classification algorithms, the LACIE aggregation and data 
adjustment prOGedures, The operational structure for simulating these ele- 
ments is shown in Figure 1. 1-1 and consists of seven key programs as 
follows: 

(1) LACIE Utility Maintenance Process 

(2) System Error Executive 

(3) Ephemeris Generator 

(4) Access Generator 

(5) Acquisition Selector 

(6) LACIE Error Model (LEM) 

(7) Post Processor 

Input data is processed by two programs --the LACIE Utility Main- 
tenance Processor and the System Error Executive. These programs 
process the input data, perform data compatibility checks, and form the 
basic data files for the LPP, 

The DAPTS simulation is performed by three programs --the 
Ephemeris Generator , the Access GenexatGr, and the Acquisition Selector 
which perform the data acquisition and seleetioH functioHS, 


BUPmODUCBIUTY 

OEIGINAL PAGE IS POOR 





Fig axe 1. 1-1, FOC l^ACtS Performance Predictor 
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The LACIE Error Model (LEM) performs the prime analysis task 
and is a single program that contains CAMS, YES, CAS, and the 
Segment Truth Generator. This phase contains the Monte Carlo con- 
trol and allows a number of operating modes. Each module within this 
phase can be Monte Carlo'ed separately without changing the others and 
the results aggregated in CAS. Modules may also be combined for Monte 
Carlo runi. The CAMS Simulator may be run With or without the segment 
acquisition statistics to give the sampling error effects with or without 
acquisition constraints. The CAS Simulator accumulates the results 
of the repetitive Monte Carlo trials and output the statistics. 

The post processor is a program that interfaces with several output 
files. It has the capability tP output special reports from its library 
of reports tp assist the user in post run analysis. 

1.2 ENVIRONMENT 

This program is designed to run on the UNIVAC 1110 utilizing the 
EXEC Vni operating system as installed at NASA, Houston,, Texas. The 
program has been coded entirely in FORTRAN V. The only special 
features utilized are the N TRAN package and the define file capability for 
direct access l/O. Most of the programs will utilize core storage of 24K 
or under. The only exceptions to this are LUMP which uses 26K, and 
LEM which uses about 48K. The most disk and drum files required at any 
time (not card and printed output) is 15 as utilized in the LEM program. 

The most file space required at one time (worst case is LEM assuming 
4000 segments) is 4, 500, 000 words. Only one 7-track taps unit is required 
for SAGE, on option, and for GRID. The largest disk file is the TACQ inter- 
mediate LEM work file with 2,040,000 words (assuming 2000 training 
segments). 

It is assumed that the users of this program know how to unload and 
load disk files to and from save tapes since most of the files are too large 
to be cataloged. The EXEC VUi job control cards for running each pro- 
gram on the UNIVAC IIIO will have to be designed and prepared. 


28234-6025-RU-00 
Page 4 


£ach program processes one country at a time so that each country 
must be handled as a separate case. A four digit case number is provided 
to identify all data files and output reports. A system of case number 
identification is required to uniquely identify a study, country and program 
user, 

It is assumed that all production usage of the programs will be per- 
formed in the batch mode even though some of the programs are small 
enough to be run from a terminal. 
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2.0 APPLICATION DESCRIPTION 

2. 1 LIST OF ALL PROGRAMS 
2. 1. 1 EPHEMS 

This program will compute orbital parameters for up to two vehicles 
orbiting about the earth for up to 549 days. This data is stored on disk 
files and represents a continuous swath path about the earth. 

2. 1. 2 GRID 

This program reads the NASA global weather tape. It processes 
the latitude and longitude in order to compute an associated grid row and 
column number representing a stereographic projection. It then writes 
the associated index number from tape into the grid matrix. 

2. 1. 3 LUMP 

This program processes substrata statistical and historical data 
cards, segment ID data cards and crop calendar data cards in order to 
generate the basic data files, CROPW, SUBHST and SEGID, for use by 
other programs. The various data fields are checked for valid entries. 

The identification for each card set is checked for compatibility with each 
other. If there are any errors, an error report is produced. 

2.1.4 SEE 

This program processes three input card sets to produce the three 
data files, YESERR, GAMERR and SIGEXT for use by the LEM program. 
Each data field is checked for validity and the ID of each card set is 
checked against the substrata historical file to detect missing data and 
incompatibility. If there are any errors, an error report is prodUeed. 

2. 1. 5 SA GE 

Given the swath table and the segment ID file, this program deter- 
mines how many times and under what conditions each segment is accessed 
by the satellites. The prog;ram generates the segment reference file, 
SEGREF. On optipn, a utility funetiGn can be performed to generate the 
weather data file, WEATHR, from the NASA weather tape. 
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2. 1. 6 SACS 


This program reads the segment reference and input control data 
fils and processes each segment using the crop window file to determine 
which accesses are legal. It lists all valid access for a segment and 
writes this information out on the acquisition file, ACQUIS. 

2.1.7 LZM 


LEM is an executive program which controls the operation of 
several application subprograms. Through these subprograms it simu- 
lates sample segment classification, strata yield estimation and production 
aggregation. LEM controls repetitive Monte Carlo trials based on input 
error distributions to obtain statistical estimates of the wheat area, yield 
and production at different levels of aggregation. The subprograms under 
its control are segment truth generator, CAMS, YES and CAS. 

2. 1. 7. 1 Segment Truth Generator 

This subprogram generates the true proportion of wheat and the true 
proportion of mixed pixels for each sample segment from the substrata 
historical file and segment ID file. 

2. 1. 7. 2 GAMS 

This subprogram provides an estimate of the proportion of wheat in 
each segment selected by the SACS program. Up to twenty-five estimates 
are allowed for traininr;; segments and one estimate for ordinary segments 
in each of four crop windows. Four types of error are accounted for in 
obtaining the estimates. CAMS gives a Monte Carlo simulation of three 
types of errors; 

1. Crop Calendar Errors 

2. Signature Extension Errors 

3. Glassification Errors 

A fourth item considered is multi-temporal sampling effects. 

2. 1. 7. 3 YES 

TMs subprogram simulates the yield estimation process of the 
LACIE system. It generates the yield estimates at the strata level of 
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from one to six estimation points in a simulation season. The estimates 
are computed from the true yield given in the input data, taking intO' 
account the effects of various estimation errors. 

2. 1. 7. 4 CAS 

This subprogram models the LACIE aggregation technique including 
the aggregation of wheat area and production from the strata to the country 
level and the estimation of the accuracy of the aggregation. The aggrega- 
tion is performed at various prediction intervals during the growing season. 
It also compares estimates to the truth baseline to compute actual errors. 

2. 1. 8 POUT 

This program is an output report generator. Upon direction from 
the user via input control cards, the following four types of reports are 
produced. 

1. Substrata Reference Data Report 

2. Population, Standard Deviation and Histogram Reports 

3. Histograms of Monte Carlo Statistics Reports 

4. Frequency of Sample Segment Acquisitions Report 

2.2 PROGRAM EXECUTION BLOCH DIAGRAMS 

The following block diagrams are given for one program at a time 
and in the basic order of execution. File names and program names will 
be used. See Section 2.4 and 2. 1 respectively for explanations, 

2. 2. 1 EPHEMS 
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2.2.2 GRID 


Q?’"" 



2.2.3 Card to Disk Utility 

This can be a system utility or any routine which can copy card 
images to disk as formated 80 character records. 


n — — 

^ \ 

T - \ 


Substrata 


Substrata 


Segment 


Crop Cal. 

Hist. Data 


Statistical Data 


ID Data 


Data 



Card sets I l&rough 5 must be terminated with a 'ZZZZ' in C. G. 6-9 of 
the last card. 

Card sets 6-8 must be terminated with a ’ZZZZ' in C. C. 5-8 of the last card. 
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2.2.4 LUMP 



2.2.5 SEE 


fij 

o 

5 a, 

.CO 


2 6 ' 
■ S 








2.2.6 SAQE 
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2.2.7 SACS 
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2.3 PERFORMANCE 

The processing times and core usage given below for each LACIE 
program are for a UNIVAC 110 8 with an EXEC R operating system. The 
runs used to obtain the data are from the handoff case which is part of 
the delivery package. Some of the basic sizing parameters are as follows: 

1 Region 

3 Zones 
14 Strata 

73 Substrata 
65 Segments 
50 Weather grid points 

40 Day swath with one vehicle used for acquisitions (two 
vehicles for swath generation) 

6 Prediction dates (ineluding four bio -windows) 

4 Monte Carlo iterations 
2. 3. 1 EPHEMS 


Time; 

2 min. 32 second 

Core : 

24059 jQ words 

2. 3. 2 GRID 

Time : 

7 seconds 

Core : 

11763 jQ words 

2. 3, 3 LUMP 

Time : 

25 seconds 

Core: 

26257 jg words 

2. 3.4 SEE 

Time: 

10 seconds 

Core ; 

14439 1 r, words 
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2. 3. 5 SAGE 

Creation of weather file time: 11 seconds ^ 

Regular run time: 1 min. 44 seconds 

Core: 21758 words 

2.3.6 SAGS 

Time: 9 seconds 

Core: 12625 words 

2. 3. 7 ;LEM 

Time: 5 miia. 40 seconds 

Core; 43 867 words 

2. 3, 8 POUT 

Core: 23308 words 

For RPTYPE = 1, Substrata ReferenGe Report 
CPU time - 11 sec. 

For RPTYPE = 2, Yield and Production Reports 
GPU time - 28 see. 

For RPTYPE = 3, Monte Carlo Reports 
CPU time - 6 see. 

For RPTYPE = 4, ■Acquisition Data Report 
CPU time -!■ 9 see. 

2.4 DATA BASE DESGRIPTION 

All files Used by LACIE programs and providing interface with the 
user or other programs are listed in Figure ?, 4-1. This figure is followed 
by a detailed description of each bii^ry file providing interface between 
programs. 


Name 

Univac 

No. 

— — — 1^'— ^ " 1 

DesGription 

Input 

Case 

Name 

Source 

Destination 

SWATH 

8 

Swath. Table /Sear ch Files 

ICSESW 

EPHEM 

SAGE 

SWATHR 

9 

Swath Reference File 

ICSESW 

EPHEM 

SAGE, POUT 

SEGID 

1 

Segment ID File 

ICSESG 

LUMP 

SAGE, LEM 

CROPW 

2 

Crop Window File 

ICSECW 

LUMP 

SACS, CAMS 

SUBHST 

3 

Substrata Historical File 

ICSESH 

LUMP 

CAS, LEM, 






POUT, SEE 

WEATAP 

13 

NASA Weather Tape - Monthly 


NASA 

SAGE 

SEGREF 

10 

Segment Reference Data File 

ICSESR 

SAGE 

SACS 

WEATHR 

H 

Weather File /Daily 

ICSEWF 

SAGE 

SAGE 

ACQUIS 

12 

Data Aequisitibn File 

ICSEAC 

SAGS 

CAMS. POUT 

CAMSF 

7 

CAMS Output File 

ICSECO 

CAMS 

CAS 

CAMERR 

8 

CAMS Error File 

ICSECE 

SEE 

CAMS 

CASF 

14 

CAS Cum Output File 

ICSECU 

CAS 

POUT, CAS ’ 

YESOUT 

10 

YES Output File 

ICSEYS 

YES- 

CAS 

SiGEXt 

9 

Signature Extension File 

ICSESE 

SEE 

CAMS 

YESERR 

11 

YES Error Model File 

ICSEYM 


YES 

SEGTRU 

13 

Segment Truth File 

ICSES T 

LEM 

CAMS 

CASDIS 

4 

CAS Distribution Output 

ICSECD 

' CAS 

POUT 

CRPERR 

8 

Crop Calendar Error Data Cards 


CARDS 

LUMP 

. INDMAT 

4 

Index Matrix Location File 


GRID 

LUMP 

SHISTD 

9 

Substrata Historical Data Cards 


CARDS 

LUMP 

SSTA TD 

10 J 

Sub strata Statistical Data Cards- 


CARDS 

LUMP 

SEGIDD 

11 

Sample Segment Location Data Cards 


CARDS 

LUMP 5 

SEGCWD 

12 

Sample Sub strata Crop Calendar Data Cards 


CARDS 

LUMP *0 

YESECD 

- _ . 1 ; 

YES Error Data Cards 


CARDS 

SEE. Ci 

CAMEGD 

- 4 ■ 

CAMS Error Data Cards 

! 

CARDS 

SEE 

SIGECD 

2 

Signature Extension Data Cards 


CARDS 

SEE 


Figure E. 4-1. Logical File Assignments 
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Name 


Uni vac 
No. 


Description 


Input 

Case 


Source 


Destination 


CASDSF 

TACQ 


CRPINt 
SU BiNT 
SEGINT 


CAS intermediate Work File 

CAMS Intermediate Work File for Training 
Segments 

Crop Calendar Intermediate File 
Su'bsitrata Intermediate LUMP Data 
Segment Intermediate LUMP Data 


CAS 

CAMS 


LUMP 

LUMP 

LUMP 


CAS 

CAMS 


LUMP 

LUMP 

LUMP 



Figure 2. 4-1. Logical File Assignments (cont'd) 
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2.4.1 Swath Table - SWATH 

Vehicle orbital data is stored for one orbit per day around the earth. 
The data represents a swath or path around the earth. There is data for up 
to 2 vehicles and for 549 days. There is also a record for each vehicle which 
contains the minimuTn and maximum A longitude. 

Access Method: Direct, fixed length records uses FORTRAN V direct access 

routine s 

Status: Permanent 

Sort: Vehicle number, then day 

A detail number is computed by the equation: 

Ree. No. = (NODAY + 1) * (VEH NO - 1) + DAYNO + 1 

A trailer Re c. No. is computed by: ' 

Rec. No. = (NODAY + 1) * (VEH NO - 1) + NODAY + 2 
Media; Disk - FASTRAND 
Record Formats: 

Header Record 


NAME (2) 

- 

8 "bytes, 

file name = 

SWATH TB 

ICASE 


2 bytes 


* 

IVEH 


2 byte s 



NODAY 

- 

2 bytes 



NLAT 

' - 

2 byte s 

- No. of latitudes; max. of 100 in a 




Record 


INLAT(2) 

- 

4 bytes 



1SEAT(2) 

- 

4 byte s 


IV1TIM(1) 


2 bytes 

- Year 


IVITIM(2) 

- 

2 bytes 

- Month 

» Start date Vehicle 1 

IVrriM(3) 


2 bytes 

- Day 


IV2TIM(1) 

- 

2 byte s 

^ Year 


IV2TIM(2 ) 

- 

Z bytes 

Month 

Start date Vehicle 2 

IV2TIM(3) 

- 

2 byte s 

- Day 1 


ISOSTR 

- 

2 byte s 




Total - 38 bytes or 17 words + 2762 bytes filler 
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Detail Record 

Each item is stored as an array for all its latitudes: 

LATNO - Latitude No. in deg, 2 bytes, range is 
0 to ±65; + - Northern Hemisphere 
- - Southern Hemisphere 


LALT 


TIME (2) - 


D LONG (2) - 

TIME(l) and (3) - 


Vehicle alt. from ALT, 2 bytes, 
KILOM^iO, range is 1000-15000 
Vehicle time in sec, 4 bytes, range is 
0-86, 400 

Vehicle A longitude in radians, 4 bytes, 
range is 0 to 2ir 

Swath latitude crossing time extremes, 
4 bytes each 


DLONG(l) and (3) - Swath latitude crossing longitude 

extremes, 4 bytes each 

The data for TIME, then D LONG, is the order of storage. 
Total Maximum Length = 100*^28 = 2800 bytes or 800 words. 


Trailor Record 

Gontains data for all latitudes; for one latitude: 

LATNO - See above. 

DLONMN - Minimum A longitude, 4 bytes, 
radians, range is 0- 2tt 
BLGNMX - Maximum A longitude , 4 bytes, 
radians, range is 0-2tt 
Filler - 1800 bytes 

Block Factor: 1 

File Size: 878,405 words for 2 vehicles and 549 days 

Usage: This table is used by the segment access generator to determine 

and record for a given segment each access through a given period 
of time. "SAGE" will be given 1 start date which must be s the 
date in the Header Record. The LATNO entry along with ISLAT 
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t 

and INLAT are used to determine if a segment is in or at any ,■ 
latitudes in the band of latitudes on the file. ISOSTR is used 
to stagger the southern latitude band with respect to the 
northern latitude band. This table is used in conji nction with the 
swath reference table. 
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2.4.2 Swath Reference File - SWATHR 

Contains nodal data for each rev, day and vehicle. This data forms 
part of the swath table and is used in conjunction with it. 

Access Method: Sequential, 1099 fixed length records on a FASTRAND file 

Status: Permanent 

Sort: Vehicle number, day 

Media: Dikk FASTRAND 

Record Format 

Header Record 


NAME (2) 

- 

8 bytes, file name = SWATH RE 

ICASE 

- 

Z bytes 


IVEH 

- 

2 bytes 


NODAY 


2 byte s 


IVlTIM(l) 


2 bytes - Year 


IVITIM(2) 

- 

2 bytes - Month 

1 for Vehicle 1 

IV1T1M(3) 


2 bytes - Day 


IV?.TIM(1) 

- 

2 bytes - Year j 


IV2TIM(2) 


2 byte s - Month j 

for Vehicle 2 

IV2TIM(3) 

- 

2 byte s - Day 


GARUSd, 1, i) 

~ 

a- semi- major axis, 8 bytes, earth radii, 
range is t to 2; for 1 £ i ^ 7, (veh) 

CARDS (2, 1, i) 


e-eecentrieity, 8 bytes, 
. 15; for 1 < i < 2 (veh) 

range is .00001 to 

C AR US ( 3, 1, i) 


i- inclination, 8 bytes, radians, range is 
1 to 2 tt; for 1 e i ^ 2 (veh) 

XOMEGD(i) 

- 

•• 

Q, 8 bytes, radians/seci 

i for 1 ^ i e 2 (veh) 

XWD(i) 


ttJ. 8 bytes, radians/ sec; 

; for 1 ^ i ^ 2 (veh) 

SA(3) 

■ - 

See input definition for card data 

Total 

.. 

24 words + 31 filler 



Detail Record: For 17 revolutions: 
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KVEH 

KDAY 


wpu, -.) 


NREV 


TMNODE(l)-{1.7) 


CAR US (4, 1, -) 


Current vehicle number, 2 bytes, range is 
1 to 2 

Current day number, 2 bytes, range is 1 to 

549 

to- argument of perigee, 4 bytes, radians, 
range is 0-2.it 

Number of revolutions, 2 bytes, range is 
1 to 17 

Time ascending node crosses equator for 
this rev. In seconds, 4 bytes, range is 
0 to 86, 4D0 

Longitude of ascending node for this rev. 

In radians, 4 bytes, range is 0 to 2tt 


17 


Ground Station Indicators - 2 bytes, use undefined for DAPTS, 

17 entries 


Record Length: 55 words 

Recommended Block Factor: 4 

File Size.' 197, 820 bytes or 60, 445 words 

Usage: This file is used in conj unction with the swath table by SAGE 
to determine when and if a segment is accessed. The swath 
reference file is searched for a day and rev entry within a 
supplied delta of a given nodal logitude and time. The day 
obtained in this manner is used to read the desired record 
from, the swath table. The calculations on the data can then 
be performed to obtain a further check for access. 




— 
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2, 4. 3 Segment ID File - SEGID 

This file is generated from NASA data cards fox use in the SAGE 
and I»EM program. It is generated in the LUMP program. 

Access Method: Sequential with fixed length records. 

Siatas: Semi -permanent. Changed only when the segment data base is 

changed. 


Sort: Country, then region, then zonOi then strata, then substrata, then 
segment, 4000 records, 

Media; Disk - FA STRAND 

Record Formats: 


Header Record; . 

Name(2) - 8 char, , file name - 'SEGMENT', 4 char, /word 
ICASE - Case no. 


ITSFG - 1 word integer = 0 all segments are training segments 

4 0 both ordinary and training segments 
present 

13 word filler 


Detail Record: 

Country ID - 4 bytes, 4 alpha eharaeter s 

Region ID - 1 word integer , 2 digit no, , I to 10 

Zone ID - 1 word integer, 3 digit no. , 1 to 100 

Strata ID - 1 word integer, 3 digit no. , I to 500 

Substrata - 1 word integer, 1 to 3200 

Segment ID - 1 word integer , 4 digit no, , 1 to 4000 

Training Seg Ind - 1 word integer, 1 - normal, 0 - training 

Training Segment Priority List - 6 words, each entry is the 
, segment no. of a training 

segment. The segments are 
• listed highest priority first. 
Less than 6 entries terminated 
by a 0. 

Latitude 1 word fit. pt, , in radians, lir/2 
Longitude - 1 word fit. pt. , in radians, 

Grid No. - 1 word integer, 1-16000 

Spring/Winter Flag - 1 word integer, 0 for winter, 1 for spring 


Total; 17 words 
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Last detail record has *ZZZZ' in the first 4 bytes followed by 
,0's in the rentiainder of the record. 

Recommended Block Factor: 20 

File Size: 68, 034 words 

Usage: Used by SAGE to get Lat/Lon and grid no. Used by LEM to place 

Lat/Lon on segment truth file. 


2.4. 4 Crop Calendar File - CRQPW 
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This file obtained from NASA data cards and is rewritten on a disk 
file for Gonvenient processing by the LUMP program. It contains the crop 
window start- stop dates for 4 crop growing windows. 

Access Method: Sequential with fixed length records. 

Status: Updated very infrequently. 

Sort: Country, then region, then zone, then strata and then substrata. 

There are a maximum of 3200 substrata records. 

Media: Disk - FASTRAND 

Record Formats: 

Header Record: 

Name(2) - 8 char. , file name - ’CROPWIND’, 4 char, /word 

IGASE - 1 word integer, 4 digit number 

IWIN - I word integer, no. of crop windows = 4 

Total - 4 words + 29 words filler 

Detail . cecord: 

Country ID - 4 bytes, 4 alpha characters 
Region ID - 1 word integer, 2 digit no. , 1 to 10 
Zone ID - 1 word integer, 3 digit no. , 1 to 100 
Strata ID - 1 word integer, 3 digit no. , I to 500 
Substrata ID - 1 word integer, 4 digit no. , 1 to 3200 
Data for Winter Wheat 

4 sets of 2 items - 1 set for each window 

Predicted Start Date 1 word integer, Zulu date - no. of days 

since 1950 

Predicted End Date - 1 word integer, Zulu date - no, of days 

since 1950 ’ 

Std Deviation of Seg. Calendar Error - 1 word integer in days 

Range: 0 — 5* f99 

Crop Window Error - 5-1 word integers, A no. of days in error. 

The 1st entry is for ist window start. The 
2nd -4th entries are for end of i-1 window and 
start of ith window (2^ i^4). The last entry 
is for (he 4th window end date. Range: 0 4:99 

Data for Spring Wheat 

Repeat of same type of data as for winter wheat 
Total; 33 words 

Trailer Record; Has 'ZzZZ' in the first 4 bytes followed by O’s. 
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ReGommended Block Factor: 10 
File Size: 105, 666 words 

Usage: This file is used by SACS in conjanction with the segment reference 

data fUe to determine if an access has occurred. It is also used by 
CAMS to determine crop calendar errors. If the first entry of a 
opring or winter start date is 0, then that data set is missing. 
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2,4,5 gubstrata Historical File - SUBHST 

This fUe is generated in the LUMP program from the substrata 
historical card data supplemented by intermediate calculations, it is 
used in the CAS analysis program. Data for one country will be stored. 

Access Method: Sequential with fixed length records. 

Status: Permanent, off loaded on tape, 

Sort: Country, then region, then zone, then strata, and then substrata. 

A maximum of 3200 substrata data records. 

Media: Disk - FASTRAND 

Record Formats: 

Header Record: 

Name(2) - 8 char., file name 1 'SUB HiST’, 4 char, /word 
ICASE - 1 word integer, case or ID no. , 4 digits 
IMXSEG - 1 word integer, maximum no. of segments in any 
substrata - 1 to 150 

Total Length: {4 + IMXSEG + 14) words 

Detail Record - for each substrata; 

Country ID - 4 bytes, 4 alpha characters 

Region ID - 1 word integer, 2 digit no. , 1 to 10 

Zone ID - 1 word integer, 3 digit no. , 1 to 100 

Strata ID ^ 1 word integer, 3 digit no. , I to 500 

Substrata ID 1 word integer, 4 digit no. , 1 to 3200 

NSEG - Mo. of sample segments in this substrata 1 to 150 

IDSEG(IMXSEG) - List of sample segments in this substrata, 

each entry is a 4 digit integer. 

Group No, - 1 wor d integer 

1 - One or more sample segments in substrata 

2 - Associated with I or more sample segments in the strata 

3 - No samjple segments in substrata or associated With the 

strata. 

■ . ir v 

PW - Historical proportion of wheat, fit. pt. in % 
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2 

Land Area - Fit. pt, , land area oi the substrata in kilometers 
PWj^. “ True proportion of wheat in %, fit, pt. 

^AGR “ ^ word integer, no. of agriGuiturai segmenta in the substrata 

- 1 word integer, no. of alloGated segments in the substrata 

6PW - Fit, pt. , bias of true proportion of wheat 
6PM - Fit. pt* , ratio of true mixed pixels 

CVl - Fit. pt. , Goefficient of variation for year-to-year change in PW 

CV2 - Fit. pt. , coeffiGient of variation for within county variation of PW 

CT3 - Fit, pt, , coefficient of variation for within County variation of 

proportion of mixed pixels 

CV4 - Fit. pt. , Goeffieient of variation of multi-year historical wheat 
area. 

Total Record Length - Max is for USSR, 168 words 

Trailer Record Length: Has 'ZZZZ' in the first 4 bytes followed by 0's. 


Recommended Block Factor: 4 
File Size: Maxirtium 537, 936 words 
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2. 4. 6 NASA Weather Tape Forrpat - WEATAP 


There is one file on tape written in binary mode (odd parity) with 
sequential I/O routines. The tape is 7 track written at 800 BPI density. 

There will be 1600 physical records on tape and each record will be 500 
words long. There will be 10 iGgical index point records per physical 
record. Each logical record will be 50 words long and will contain 
BDC information in the following format: 

E 
or 

W YEARS % FREC MEAN% 
Al, 12, 12.(13,812; 1213) 


N 

or 

INDEX l-GRIEM LAT S |LONGj 

(15, 2X, 13, 2X, 13, 3X. 312, Al, IX, 13, 212, 
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2, 4. 7 Segment Reference File - SEGREF 

This file is generated by SAGE to provide a reference data base for 
all segments to be analyzed by the SACS program. This file can represent 
a maximum of up to a 426 day time period. 

Access Method; Sequential with fixed length records. 

Status: Semi-permanent, Regenerated normally less than once. per week. 

W ill normally be kept on tape. 

Sort: Country, then region, then zone, then strata, then substrata, and 
then segment, There are a maximum of 4000 segments. 

Media: Tape 

Record Formats! 

Header Record; 

Kame(2) - 8 char., file name - 'SEG RSFE', 4 char, /word 

ICASE - 1 word integer, case no, , 4 digits 

Reference Date - No. of days since 1950 (2 bytes). Used in 
conjunction with acquisition date. 

NMAX - 1 word integer, no, of days in study; range 1 to 426 

Dummy Word 

ICSESW - 1 word integer, swath files ease no. - 4 digits 
ICSESG - 1 word integer, segment ID file case no. , 4 digits 
NVEH - No. of vehicles - 1 word integer 
Total; 9 words 4 599 word filler 

Detail Record: 

Country ID - 1 word, 4 alpha Gharaeters 

Region ID - 1 word integer , 2 digit no. , 1 to 10 

Zone ID - 1 word integer, 3 digit no* , 1 to 100 

Strata ID - I word integer, 3 digit no. , 1 to 500 

Substrata ID - 1 word integer, 4 digit no. , 1 to 3200 

Segment ID - I word integer, 4 digit no. , 1 to 4000 

Training Seg. Ind. 1 word integer, 1 - normal, 0 - training 

Spring /printer Flag - 1 word mteger, 0 for Winter, 1 for spring 
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Acquisition Data Set - Repeated 75 times x NVEH 

Acquisition Day - 1 word integer, Zulu date (no. of days since 1950) 
Sun Elevation Angle - 1 word fit. pt. , radians, range +90® 

Cloud Cover % - 1 word integer, % * 10 . 

Time of Acquisition - 1 word fit. pt, , sec, 1 to 86, 400 
Total Length: 608 words 

Last detail record has 'ZZZZ' in the first word followed by O^s in 
the remainder of the record. 

Slock Factor: 2 

File Size: 2, 432, 008 words for 2 vehicles 
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2.4.8 Weather Data FUe - WEATHR 


This file contains cloud cover % for 366 days and 16000 grid points. 
This file is generated from a NASA supplied weather tape. 

Access I'lethod: Direct with fixed length records - uses FORTRAN V 

direct access routines 

Status: Semi-permancrit. Can be regenerated for each use of SAGF 

or kept as a permanent file for IOC studies. 

Sort: Grid point, then day 

To get to a particular record N which represents 5 grid points: 

N = 1 + + 1 (if remainder) 

D 

+ 0 (if no remainder) 

Media: Disk - FASTRAND 

Record Formats: 


Header Record: 

Name (2) - 8 bytes, file name = 'WEATHER' 

NMAX - No. of days in a record, 2 bytes, 1 to 366 
Total: 12 bytes or 4 words 
171 byte filler 


Deta.il Record 


For 1 grid - all days 

Each day entry is a value 0 to 8 representing 
8th' s of 100%. 8-4 bit entries are stored 


right justified in a word; for 366 days it 


would require 


36_^s^4 
8 


= 183 bytes or 46 words, 


Blocking Factors: 5 (Internal) 


File Size: 2, 9 28, 000 bytes, 736, 000 words 


Usage: This file is used by SAGE to obtain the cloud cover data. 
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2. 4, 9 Data Acquisition File - ACQUIS 

This file contains a list of accesses for each segment processed in 
the segment acquisition selector program (SACS). It is to be generated 
using the Crop Window File and Segment Reference Data File. A maxi- 
mum 426 day time period is represented. 


Access Method: Sequential with fixed length records. 

Status: Regenerated - either from rerunning SACS or loading from a 

saved tape. 

Sort: Country, then region, then zone, then strata, then substrata, and 
then segment. A maximum of 4000 records for one country. 

Media: Disk - FASTRAND 

Record Formats: 

Header Record; 

Name (2) - 8 char., file name - 'ACQUISI', 1 char. /word 
ICASF - 1 word integer, 4 digit case no. 

NMAX - 1 word integer, no. of days in study; range = 1-^426 

IWIN - 1 word integer, no. of crop windows used in study -4 

HEAD(4, 4) - 4, 4 word entries, each ehtry consists of a 16 
character crop window title 

ICCAS - 1 word integer, crop window case number 

ISCAS - 1 word integer, seg refer case number 

IFILL(84) - 84 words of filler 


Detail Record: 

Country ID - 4 bytes, 4 alpha characters 
Region ID - 1 word integer, 3 digit no. , 1 to 10 
Zone ID - 1 word integer, 3 digit no. , 1 to 100 
Strata ID - 1 word integer, 4 digit no. , 1 to 500 
Substrata ID ^ 1 word integer, 4 digit no. , 1 to 3200 
Segment ID - 1 word integer, 5 digit no. , 1 to 4000 
For each of 4 Crop Windows: 

25 entries for 

Acquisition Date - 1 word integer, Zulu date 
Total No. of Acce sses - 1 word integer, 3 digit no. 
Total Length: 107 words 
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Dummy data record with country = 'ZZZZ' trails last data record. 
Recommended Block Factor: 5 
File Size: IOC - 428, 214 words 

Usage: This file will be used by the CAMS module within the acquisition 

and analysis program. It is used to obtain the acquisition list by 
Crop Window for each segment in Uie study. 
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2.4,10 CAMS Oatpat File Description - CAMgP 

This file is ased by the CAS program and is generated by CAMS. 

Record Formats; 

Header Record: 

Name - 8 char., 3 alpha, file name - 'CAMS OUT’, 4 char, /word 

Case Number - 1 word integer, 4 digit case number 

Error Model - 1 word integer, = 0 - Model 1, = 1 - Model 2 

CAMD Error File Case Number 

ACIUISI File Case Number 

CROPW File Case Number 

SEG TRUTH FUe Case Number 

Signature Extension File Case Number 

IMULTI - Multi- temporal Option Flag 

TSCC - Crop Calendar Error Option Flag 

ICLASS - Input C las sifieation Indicator 

ISEXT - Signature Extension Flag 

lACQ - Acquisition File Usage Flag 

ICAMS - CAMS Option Flag 

Total: 4 words of filler 


Detail Record (one per segment): 

Country ID - 4 bytes, alphanumeric 
Region ID - 1 word integer, 2 digit number, J^IO 
Zone ID - 1 word integer, 3 digit number ( 1-100 
Stratum ID - I word integer, 3 digit number, 1-500 
Substratum ID - 1 word integer, 4 digit number, 1-3200 


Segment ID - 1 word integer, 4 digit number, 1-4000 
True proportion of wheat this segment, fit. pt, , % 0-100 
Zulu Acquisition Day - 1 word integer (zero for no acquisition)] 
Estimated Proportion of Wheat - Fit. pt. ; % 0-100 
Error in Proportion of Wheat Estimate, fit. pt. ; J^99. 99 
Total Length: 19 words 


Trailer Record: 


One 

ordered set 
for each of 
4 windows 


Name - 4 bytes, 'ZZZZ' 

Total: 1 Word with 18 word filler 


28234-6Q25 -RU-00 
Page 3 $ 


File Size: 78,038 words - assuming 4000 segments 

Access Method: Sequential binary file with fixed length records 

Recommended Blocking Factor: 20 

Sort: Country through segment ^ 

Media: Disk - FASTRAND 

Status: Temporary, will be regenerated for each study. 


\ 






28234-602S-RO-d0 
Page 37 


2.4.11 CAMS Error File - CAMERR 

t 

This file is generated by SEE and used by the CAMS program'. 
Sort: Country, Region, Zone, Strata, Substrata, Segment 
Media: Disk - FASTRAND 

i 

Record Formats: 


Header Record: 

Name(2), - 8 char. , file name ’CAMS ERR', 4 char, /word 
ICASE - Case no. , 1 word integer 
Total: 3 words + 47 words filler 


Detail Record: 

Country ID - 4 bytes, 4 alpha characters 

Region ID - 1 word integer, 3 digit no. , 1 to 10 

Zone ID 1 word integer, 3 digit no. , 1 to lOG 

Strata ID - 2 byte s, 4 digit no. , 1 to 500 

Substrata ID - 2 bytes, 4 digit no. , 1 to 3200 

Segment ID - 2 bytes, 4 digit no. , 1 to 4000 

Crop Window (4) - one set of data below for each crop window 

(see SEE Problem Description Input Data 
for definitions) 


P(W/W) 

- Fit. pt. 

, 0 

to 1 

P(W/lv(I) 

- Fit. pt. 

. 0 

to 1 

P(W/0) 

- Fit, pt. 

. 0 

to 1 

®w/w " 

Fit 

. pt. , 

+9. 

999 

^w/w " 

Fit. 

pt. , 

0-9 

.99 

®w/m" 

Fit. 

pt. , 

+9, 

999 

^w/m" 

Fit. 

pt. , 

0-9 

.99 

®w/o 

Fit 

• pt. , 

±9. 

999 

^w /q “ 

Fit. 

. pt. , 

0-9 

1.99 

- FW. 

pt. , : 

!l9.999 

®PW " ® 

^It. 

pt. , Cl 

1-9. 

99 
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Record Length: 50 words 

Last detail record has 'ZZZZ' in the first 4 bytes fpllowed’by O's 
in the remainder of the record. 

Block Factor: 10 

File Size: 200, 100 words, assuming a maximum of 4000 segments 
Usage: Used by CAMS 

Status: Semi-permanent, changes only when CAMS error model changes. 
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/ 

2.4. 12 CAS Cumulative OutTaut File - CASF , ' - - 

This file is generated by the CAS module and is used by the output 
processor {POUT) and by LEM for restart. It contains values accumulated 
over all Monte Carlo iterations. . ' 

Access Method: Random access with fixed length reGords, 

Status: Changed for each Monte Carlo iteration in LEM. The final 

LEM version will be kept only as needed to complete a current 
error model study. Required for restart of the LEM program. 


Format: 

Record 1 is the header record, ' 

Record 2 is the country record, 

Records 3-12 are the region records, 

ReGords 13-62 are the zone records, 

Records 63-387 are the strata records. 

Header Record Format: 

1. Nanae{l) - 6 character file name 'CASCUM’ 

2, ICASE - Case no. , 1 word integer 

3. COUNTR - Country, A6 format 

4 , NT - Current Monte Carlo iteration number 
5* NREGS - No. of regions 

.6. NZTOT - Total number of zones 

7. NSTRAT - Total number of strata 

8. NBW - No. of bio-windows; 1-4 

9. NPDATE - No. of prediction dates; 1-14 

10-13, - tip to 4 bio -window numbers 

14-27, PRDATE - Up to 14 prediction dates in Zulu time 
477 words of filler 
Total Lengthr 504 words 

Each data record contains data blocks for up to four bio -windows and 
up to 14 predietion dates. The blocks for the bio-windows come first followed 
by the bloeks for the prediction points. Missing bio-windows or prediction 
points are represented by zero fill. 


"I 


Country Data Record Format: 
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This record consists of 18 blocks, each containing the 28, quantities 
in Data Set 17. • ’ ‘ 


1. HWAG - Historical wheat area (WA) 

2. TWAC - True WA 

3. EWAC - Estimated WA 

4. AERRC - Error in WA 

5. AVARC - Variance in WA 

6. TPRODC - True production 

7. EPRODC - Estimated production 

8. PRERRC - Error in productioh 

9. PRVARC - Variance in production 

10. TYG r True yield 

11. EYC - Estimated yield 

12. YERRC - Error in yield 

13. MIC - No. of acquired Group I segments 

14. M2G - No. of acquired Group II segments 

15. CTIG - No. of Group I substrata with acquired segments 

16. CT2C - No. of Group II substra-ta with acquired segments 

17. CT3C - No. of Group HI substrata 

18. ANA VC - Analytic area variance 

19. ANPRVG - Analytic production variance 

C 2 

20. SQAERC - S (area error) 

C 2 

21. SQPERC - £ (production error) 

C 2 

22;. SQYERC - £ (yield error) 

23. GLEWA - Gonfidence level about estimated wheat area 

24. CDEPRD - Confidence level about estimated production 

25. CLATEC Confidence level about true wheat area using 

estimated variance 

24. GJbPTEC - Confidence level about true production using 
estimated variance 

27, CLATWC - Confidence level about true wheat ar ea using 

within , GO unty va r iane e 

28. CLPTWC - Confider i level about true production using 

within county variance 

Total length: 504 words 


Region Data Record Format; 
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This record consists of 18 blocks, each containing the 22 quantities 
from Data Set 16 plus identifying information. 


1 . 

2 . 

3. 


REGION - Region ID 

0 - 


:} 


2 words of filler 


4. 

NZONES 

I Number of zones in region 

5. 

HWAR 

\ 

9m 


6. 

TWAR 

• 



23, 

■ 

ANPRVR 


V Similar to country quantities 

24, 

SQAERR 

- 1 


25. 

SQPERR 

- 


26. 

SQYERR 

- > 


27. 

28, 

0 

> 2 words of filler 

0 -i 


Zone Data Record Format: 

This record consists of 18 blocks, each containing the 22 quantities 
from Data Set 15 plus identifying information, 

1. REGION - Region ID 

2. ZONE - Zone ID 

3. 0 - Filler 


4. 

NSTRAZ - No. of strata i 

n zone 

5. 

HWAZ - ' 



6. 

TWAZ - 

a 



23. 

*■ . 

ANPRVC - 

Similar to 

region and country quantities 

24. 

SQAERZ - 



25. 

SQPERZ - 



26. 

SQYERZ - ^ 



27. 

0-1 




V 2 words of filler 


00 

0-1 
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Strata Data Record Format: 

This record consists of 18 blocks, each containing the 22 quantities 
from Data Set 14 plus identifying information. 

DSET14 


1. 

REGION - Region ID 

2. 

ZONE “ Zo'jlIq ilD 

3. 

STRATA - Strata ID 

4. 

0 - Filler 

1-5. 

HWAS 

2-6. 

TWAS 

3-7. 

EWAS 

1 

00 

«’ 

AERRS 

5-9. 

AVARS 

6-10. 

TPRODS 

7-11* 

EPRODS 

8-12. 

PRERRS 

9-13. 

PRVARS 

10-14. 

YS 

11-15. 

ESTYS 

12-16. 

YERRS 

13-17. 

MIJS 

14-18. 

M2JS 

15-19. 

CTIS 

16-20. 

CT2S 

17-21. 

CT3S 

18-22. 

ANA VS 

19-23. 

ANPRVS 

20-24. 

SQAERS 

21-25. 

SQPERS 

22^26. 

sqYers 

27. 

0 

28. 

0 

Recommended Blocking Factor; 

File Size: 

195, 048 words 
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2.4,13 YES File DegeriptiQ-n * YEgOUT 

/ 

This file is created by YES and used by tue CAS program. 

Access Method: Sequential with fixed length records. 

Status: Semi-permanent - changed only when YES error model changes. 
Sojrt: Country through stratum - 500 strata s maximum 
Media: Disk - FASTRAND 
Header Record: 

Name (2) - 8 char., file name - 'YES' , 4 char, /word 
Case Number - 1 word integer, 4 digit case number 
Total: 20 words filler 

Detail Record 

Country ID - 4 bytes, 4 alpha characters 

Region ID - 1 word integer, 2 digit nurnber, I -10 

Zone ID - 1 word integer, 3 digit number. 1-100 

Stratum ID - 1 word integer, 3 digit number, 1-500 

YSTR - True yield for stratum, fit, pt, in quintals /HEC TAR 

One set for each of 6 prediction points 

IZFRDD - Zulu prediction date, 1 word integer 

YSCl - Estimated yield for this prediction date, 1 word integer; 
0-99. 99 in quintals /HEC TAR 

VSYCl - Standard deviation of yield error, fit. pt, ; 0-99. 99 in 
quintals /HEC TA R 

Total: 23 words 
Trailer Record: 

Name (2) - 4 bytes, file ender, 'ZZZZ' 

Total-: 1 word With 22 word filler 

Recommended Block Factor: 10 

File Size: 11,546 words 
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\ 


\ . y 




2, 4. 14 Signature Ej^tension FUq 

This file is generated from cards in SEE (System Error- Executive), 
It contains coefficients for signature extension processing in CAMS for 2 
error models. 


Access Method; Sequential with fixed length records. 

Status: Semi-permanent. Only changed when the error model changes. 

Sort: Country, then region, then zone. 

50 zone records maximum. 

Media: Disk - FASTRAND 


Record Formats: 

Header Record: 

Name (2) - 8 char., file name - 'SIGEXTEN', 4 char. /word 
ICASE - 1 word integer, 4 digit case no, 

56 words of fill 


Detail Record: 

Country ID - 4 bytes, 4 alpha characters 
Region ID - 1 word integer, 2 digit no. , 1 to 10 
Zone Id - 1 word integer, 3 digit no. , 1 to lOO 


Model I Data -- Classification error eoefficients 

Coefficient set to account for probability of classifying as wheat 
given wheat. Same set as for Model 2 - 14 fit. pt. words. 

Coefficient set to account for probability of classifyinig wheat 
given mixed, 14 fit. pt. words, 

®1M*. ®2M’ ^lij^’ ^ ^ ^ 

Coeificient set to account for probability of elassifying as wheat 
given other, 14 fit. pt. words. 

®10* ®20’ ^liQV^2iQ ^ 1^6) 
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Model 2 Data - Ignore Glassification errors 

- 2 word fit, pt. , bias eoefficients 

ff. . (l^i^6) - 6 word fit, pt. , std. dev. multiplier coeffieients 

^ for each of 6 training priority segments 

a_. (l^i^6) - 6 word fit. pt. , std. dev, adder coeffieients for 

^ each of 6 training priority segm.ents 

Total; 59 words 

Last detail record has 'ZZZZ' in country ID field followed by O’s. 
Recommended Blocking Factor: 10 
File Size; 6, 018 words maximum 
Data Ranges; 

0's are between -9. 999 and +9. 999 
B's are between 0 and +9. 999 
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2,4,15 YES Error Model File 

' ~ r 

This file is generated from cards in SEE. It contains reference 
yield data and yield error parameters. It is used by the YES module. 

Access Method; Sequential with fixed length records. 

Sta,tus; Semi-permanent. Changes when new error model is required 
or when new yield reference data is available. 

Sort: Country, then region, then zone, and then strata. A maximum of 
500 strata records. 

Media; Disk - FASTRAND 

Record Formats: 

Header Record: 

Name{2) - 8 char. , file name •“ 'YESERROR', 4 char, /word 
ICASE - 1 word integer, 4 digit case no. 

20 words of 0 fill 

Detail Record* 

Country ID - 4 bytes, 4 alpha characters 

Region ID - 1 word integer, 2 digit no. , 1-10 

Zone ID - 1 word integer, 3 digit no, , 1-100 

Strata ID - 1 word integer, 3 digit no. , 1-500 

^TRUE " yield in quintals /HEC TAR, 1 word fit. pt. 

Truncation Month Yield Error Data 

In asGending order, a zero entry terminates 
6 - 3 word sets of following data 

Truncation Date in Zulu - integer 
Bias error in quintals / HEC TAR - fit, pt. ; +99*99 
Standard Deviation - fit. pt. quintals /HEC TAR; 0-99*99 
Total: 23 words 

Trailer Record: 

Has 'ZZZZ' in country ID field followed by Q's, 

Recommended Blocking Factor: 20 
File Size: II, 546 words 
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2.4,16 Segment Triith File , 

This file is generated under control of the LEM processor in order 
to provide segment truth data for CAMS module. 

Access Method: Sequential with fixed length records. 

Status: Changed for each Monte Carlo iteration in LEM. The final 

LEM version will be kept only as needed to complete a current 
error model study. 

Sort: Country, then region, then zone, then strata, then substrata, then 
segment. A maximum, of 4000 records. 

Media: Disk - FASTRAND 

> 

Record Formats: 


Header Record: 

Name(2) - 8 char. , file name - ’SEGTRUTH*, 4 char, /word 
ICASE - Case no. , 1 word integer 

ITSFG - 1 word integer = 0 all segments are training segments 

4 0 both ordinary and training segments 
are present 

12 words filler 
Detail Record: 

Country ID ~ 4 bytes, 4 alpha characters 

Region ID - 1 word integer, 2 digit no. , 1 to 10 

Zone ID - I word integer, 3 digit no. , 1 to 100 

Strata ID - 1 word integer, 3 digit no. , I 'to 500 

Substrata ID - 1 word integer, 4 digit no., 1 to 32C0 

Segment ID - 1 word integer, 4 digit no. , 1 to 4000 

Training Seg. Ind. - 1 word integerj 1 - normal, 0 - training 

Training Seg. Priority List - 

6 words, each entry is the segment no. of a training segment. 
The segments are listed highest priority first. Less than 6 
entries are terminated by a 0. 

Spring/W inter Flag - 1 word integer, 0 for winter, 1 for spring 
True Proportion of Wheat - 1 word fit. pt. , PW in % 0-100 
True Proportion of Mixed - 1 word fit. pt. , PM in % 0-100 
Total Record Length: 16 words 
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Trailer Record: 

Has 'ZZZZ* in country ID field followed by O's. 
Recommended Blocking Factor; 20 


File Size: 64, 032 words 
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2. 4. 17 CAS Distribation Oatpat File - CASDIS 

This file is generated by the CAS module arid is used by the output 
processor (POUT). It contains certain quantities (errors and confidence 
levels) for each Monte Carlo Iteration. 

Access Method: Random access with fixed length records. 

Status: Changed for each Monte Carlo iteration in LEM. The final 

LEM version will be kept only as needed to complete a current 
error model study. 

Format: 

Record 1 is the header record, 

Records (2-4) + 63 (PN-1) are the country records (3 records are 
required). 

Records (5-14) + 63 (PM-1) are the region records. 

Records (15-64) + 63 (PN-1) are the zone records; where PN is 
prediction point no. 1 ^ PN ^ 18. 

There are 1 + 18 x 63 = 1135 records. Each set of 63 records after 
the first record belong to a prediction point. 

Header Record Format: 

1. Name(l) - 6 character filename 'CASDIS* 

2. ICASE - Case no. , 1 word integer 

3. COUNTR - Country, A6 format 

4. NT - Current Monte Carlo iteration number 

5. NREGS - No, of regions 

6. NZTOT - Total number of zones 

7. FILLER 

8. NBW •- No. of bio-windows; 1-4 

9. NPDATE - No. of prediction dates; 1-14 

10-13. - Up to 4 bio -window numbers 

14-27. PRDATE - Up to 14 prediction dates in Zulu time 

276 words of filler 


Total Length: 303 words 
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Format of First Country Data Record: 


W Qrd(s) Description 


1. 

AEREFC 

Reference value for area error 

2. 

PEREFC 

Reference value for production error 

3. 

YEREFC 

Reference value for yield error 

4-103. 

AERRG 

th 

Word n+3 specifies the area error for the n 
Monte Carlo iteration 

104-203. 

PRERRG 

Word n+103 specifies t&e production error for 
the n^^ Monte Carlo iteration 

204-303. 

YERRC 

th 

Word n+203 specifies the yield error for the n 
Monte Carlo iteration 

Format of 

Second Country Data Record; 

WQrd(s) 

Description 

1. 

CLEARF 

Reference value of the area Gonfidence level 
Est/Est 

2. 

GLEPRF 

Reference value of the production confidence level 
Eet/Est 

3. 

CLTARF 

Reference value of the area confidence level 
True/Est 

4-103. 

CLEWA 

Word n+3 specifies the area confidence level 
Est/Est for the n^^ iteration 

104-203. 

CLEPRD 

Word n+103 specifies the produetion confidenGe 
level Est/Est for the iteration 

204-303. 

CLATEC 

Word n+203 specifies the area confidence level 
True/Est for the n^^ iteration 
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Format of Third Country Data Record: 


Word(3) Description 


1. 

CLTPRF 

Reference value of the production confidence 
level True/Est 

2. 

CliAWCR 

Reference value of the area confidence level 
True/WC 

3. 

CLPWCR 

Reference value of the production confidence 
level True/WC 

4-103. 

CLPTEC 

Word n+3 specifies the production confidence level 
True/Est for the nt^ iteration 

104-203. 

CLATWC 

Word n+103 specifies the area confidence level 
True/WC for the n**^ iteration 

204-303. 

CLPTWC 

Word n+203 specifies the production confidence 
level True/WC for the ntia iteration 

Format of Region or Zone 

Data Records; 

Word(s) 

Description 

1. 

AEREFR 

Reference value for area error 


or 

AEREFZ 

• 

2. 

PEREFR 

Reference value for production error 


PEREFZ . 


3. 

YEREFR 

Reference value for yield error 


or 

YEREFZ 


4-103. 

AERRR 

or 

AERRZ 

Word n+3 specifies &e area error for the 
Monte Carlo iteration 

104-203. 

PRERRR 

or 

PRERRZ 

Word n+103 specifies the production error for 
tlie n^ iteration 

204-303. 

YERRR 

or 

YERRZ 

#■ 

fell 

Word n+203 specifies the yield error for the n 
iteration 


Recommended Blocking Factor: 1 
File Size: 343, 40‘5 words 
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2,4,18 Index Matrix Lcjcatioii. File - INDEXMAT 

This file contains grid numbers or INDEXS for global coverage. Each index is 
associated with an (I,J) element corresponding to the stereographic projection 
coordinate of a given latitude and longitude. . ' 


Access Method; 


Status : 

Sort; 

Media: 

Record Formats: 
Header Record: 

Detail Record: 
Blocking Factor; 


Direct with fixed length records - uses FORTRAN V direct 
access routines. 

Almost permanent, may never be regenerated. 

By I, J value 1<J^2§0 and for each >3, 1<I<50Q. Each record 
corresponds to a value of J. There are 251 records counting 
the header record. 

Disk-FASTRAND 


Name (2) - File Name = 'INDEXMAT' 499 words of filler 

For 1 value of J there are 500 index values, 1 for each I, 1 
to 500. The index value is an integer 1-16000. 

1 


File Size: 125,500 words 

Usage; This file is used as input to LUMP to obtain the index value 

for each segment. 






♦ 
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3. 0 PROCEDURES AND REQUIREMENTS FOR EACH 
PROGRAM IN LACIE 


3. 1 EPHEMS 
Operational Assamptions 

• There can be a maximum of two vehicles, 

• There can be a maximum of 549 days. 

• The latitude band for swath table generation is +_65° latitude 
and there will be a maxiTnum of 100 latitude points to process. 

• Only one case is run at a time. 

• The complete card data set must be entered for each run, 

• The existing swath tables will be destroyed each time this 
program is run. 

• The swath tables will be designated as permanent files since 
it is assumed that this program will be run about once per 
year. 

• F.^ed field card formats are utilized. 

• A sun synchronous inclination will be computed for each z'Vin 
condition. 

• Only first order oblateness terms are modeled in the orbit 
generation; i. e. , nodal regression and apsidal precession, 

• Orbital elements are maintained as Keplerian elements 
(a, e, i, a, 0 ), t) and are updated on a rev-by-rev basis, 

• Orbit initialization equations will insure passage through a 
reference point at a specified time on the first day. Active 
passage is assumed to be on the descending pass. 

3. 1. 1 Input Card Data 

3. 1. 1, 1 Input Data DescripMpn 

See Tables 3. 1-1 and 3. 1-2. 


Dimension 

1 


Table 3. 1-1. EPHEM Program Control Input 

( Nominal | t i 

Value 1 Kange j Units ) 


2 

Blanks 


F 

(Blank) 


1-9999 

1-549 

1-2 


T or F 


5,65 

0-65 

Deg 

15,45 i 

0-65 

Deg 

1,549 

1-549 


8*0 



6 

0-12 



Description 


A 4 digit case number to identify the 
printed output and the swath tables 


No. of days to generate ephemeris data 
for the swath table 


Number of vehicles to process 


Provision for an 72 character case header 

to print out at the top of every output print 
page 


Flag to allow printout of intermediate data 
for swath calculations. Data for the 
following rouitnes is printed based on 
subscript values: 

1. SWATH 2. Swath table and 
4. REV TAB reference file 


Northern Hemisphere latitude band for 
swath generation 


Southern Hemisphere latitude band for 
swath generation 


Ephemeris data display flag 
= [start day, end day] pairs 


No. of months delay to start the genera- 
tion of swath data on the swath table for 
the southern latitudes 


1 
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Table 3. 1-2, Orbital Determination Input 



Units 


Description 


66 5 G to 
7700 
10 - 5 <. 15 


0-360 
0-360 
+ 65° 


- 10-+10 


a 8 .'V - vfiwi'f 


Kilom. 


Deg. 

Deg. 


Vehicle 1 reference latitude passage time 
IVlTIM(l)-(Year-l9O0)j IVlTIM(2)-Month no< 
IV1TIM(3)-Day No. ; IV 1TIM(4) -Hours i 

IV 1 TIM (5) -Minutes; IVlTlM(6)-Seconds 1 


Vehicle 2 reference latitude passage time; 
same information as Vehicle 1. 


Orbital elements for Vehicle 1 
Semi -major axis 


Longitude of reference latitude 
Argument of perigee 
Vehicle reference latitude 


Orbital elements for Vehicle 2, 
Same data as for Vehicle 1. 


Heference scan angles defining swath 
field of view 


TO 

(I) 


to 00 

in I 
ui Ox 

o 

CO 

in 


s 

« 

O 

o 




■:a 


I 

:1 


J 
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See Table 3. 1-3. 

3, 1, 1, 3 Deck Sgtap 

The order of inpai; is: 

1. Header card 


2, Program. Control Card 

3. Orbital determination data (2 cards) 

3. 1. 1, 4 Rules for Entering Data on Cards 

3. 1. 1. 4. 1 General 

1.. Integers must be entered right Jastified. 

2. F -format numbers must have the decimal point present, i. e, , 
F5. 1-XXX.X, F4. 0-XXX. 

3. The card sequence numbers in the C. C, 79-8® be present 

in all data cards. 

3, 1, 1. 4. 2 Rules for Spe c ific FieMs 


IVEH - No more than 2 vehicles can be entered in this field. 

NODAY - A number larger than 549 will not be aecepted. 

lYR - An entry less lhan 64 will not be aeGepted. 

INLAT and ISLA^T - -A negative entry or an entry greater than 65 
will not be accepted. The total number of latitudes 
allowed, represented by the combined bands, is 100, 

IPPI5 - An entry larger than 549 will not be accepted. 

SA - An absolute value greater than 10^ is not allowed. 

IVITIM and IV2TDvI - The first subscript entry (year)' must be 
between 1964 and the present. 



W)i K S E-. h 

S:> i 3 3 


3.1-3. PROGRAM C 










/ 

i 




V 

1 

LOL CARD FORMAT 




;• ; 

34 37 40 43 

46 

49 

52 

79 ■ 


28 3 


Ej>hemeris Data Display Flag Pairs 


Day Day Day 



Start 

Day 



Start Stop g 

Day Day m 


SAU) 

SA(2) 

. Left 

Veh. 

Side 

Refer. 

Refer. 

Scan 

Scan 

Angle 

Angle 



I VI TIM - Veh, 1 


Y ear 


I Day 


U i O: 

^ d <u 
ffi S M 


IV2TIM - Veh. 2 


Y ear <u 

^ Day h I I 

9 S S 

■I J ^ 

A t/3 



F5. 1 14 12 13 12 12 12 14 12 13 12 12 12 


19 25 


ORBITl - Veh. 1 




CARD 2 

35 43 49 53 59 


ORBIT2 - Veh. 2 



::g 


AN - 

A 

F 

F 

L - 

L. 
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3, 1, 2 Sample Card Inputs 


TtSr CA5E FOR HAMUOFFiPyS, 

TF.S7 1 

- 40 UAt 

HEADER 

SWATH* 12/3U/7S SVAWTt2 VEHIGL 




TC#^6t iV'tH NnrjAY tmuT 
! 2 ■■^0 Tj 65 

TSLAT 
20 65 

LlOcPU 

ftff 

1PPI5 STaRI/STuP times 
1 ito 0 0 0 0 0 0 

0 

IJ3USIP 
0 1 

01 ^ 

■ SA 

.U R.R 7b 12 

IVITIm 
30 IP 

15 <Jp 

IV2TIM 

75 12 31 5 »0 “0 




oa 

nuFUTj 

72*>0.0 .OSSPo P. 3uS, 

SO. 

766«,4 

0RB1T2 

,02500 105. 200. 

20. 



03 



3. 1. 3 Oatpat Report Data Defmitions 
See Table 3.1-4. 


4 










^ . / 


Table 3. 1-4. EPHEM Output Report Data Definitions 


Report N'enue 

Symbol 1 

[ Range 

Units 

Report 

Description 

WP 

W ' I 

0 

1 

Radians 

1. 

1 

Argument of perigee 

NREV^REV 

- 

1-18 


1 3 

A • j «. 1 

Current orbital revolution number in a day 

TMNODE = (Node Time) 

,1 

. 1 

- 

G-86, 400 . 

Sec 

i 

1. , 3. 

Time daring a rev that ascending node 
crosses equator 

1 

CARDS = (Node Dongitude) ' 

* 1 

n 

0-27T 

Radians or 

1. 

Dongitude of aseending node 


0-360 ' 

deg, min. 

3. 


1 



see 



VEMCLE 

- 

, 102 

-- 

i* , Za , 3# 

Vehicle number 

DAY i 

1 

- 

1-549 

Day 

l.,2. 

Day no. of mission 

LA TNO 

X 

-65-+65 

Degrees 

2. , 4. 

latitude of a vehicle 

lALT 

- 

800-1800 

Kilom , 

2. 

Altitude abo ve the earth of a vehicle 

TIME ! 

1 

- 

0-86,400 

See. 

2. 

Swath latitude crossing time extremes 

DLONGl 

■ 

Q-27T 

■ 

Radians ^ 

2. 

Swath latitude crossing longitude extremes 

Semi- Major Axis i 

a 

6650-7700 

kilom. 

3. 


Eccenti’ieity i 

e 

-6 

10 15 


3. 

— 

Period 

P 

6OOO-^65O0 

- See. 

3. 


.1 

InGlinatiOn 

i 

70-120 

Deg, min 

3. 




. 

see 


Argument of Perigee 

0) 

0-360 

Deg, min. 

3. 

OTQ ^ 

- . o w 




sec 


m 1 

Rt. As,g. GreePwicli 

ffi 1 

0-360 

Deg, min 

3. 

■ nO 

O 




sec 


t\> 

DLONMX 

0 

0-277 

Radians , 

4. 

Maximum value of DDONGl at any latitude • 

DLONMN 


0-2# 

Radians i 

4. 

C 

Minimam value of DLONGl at any latitude i 

o 

o 

■ i ■ ■ 1 > 


1 

*: ... . 


i 
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3. 1, 4 Sample Pages from Each Output Report 

1. See Figure 3. 1-1 for Sample Swath Reference Record. 

2. Soe Figure 3. 1-2 for Sample Swath Table Record 

3. See Figure 3. 1-3 for Sample Ephemeris Dis’>lay 

4. See Figure 3. 1-4 for Sample Table Search Record 



TEST CASt FOR H*Nh0Pf»SV5. TeSI Ji - 40 OaY S»iAlH, 12/30/75 START»2 VEHlCL LPR STMULATrOR CASE 1 PAGt 1 


>r 


1 

'1 

S/,AIH REFfcPtNCE 

ptrURo 

vf:hjcL£ 1 

DAV 1 



U\i 

isii£.y 


GAkUS 



.*0191o*01 

1 


U J 





eiAlSil+Ob 

.446247+01 





^2?3S36+05 

.401 15C. + 01 





• } 4 OS 



1 


b 

,30/546+05 

.310V72+01 


' ■ 


4> „ 

_,JiOV*>5.1+05™ 

.?ej5Hf»l + 01 


'T 

1 


f 

.471557+05 

*?2n7HVtUl 


'i . 


a 

,5335t>2+0b 

ol/So<?ti + 0l 




V 

_«.5255t>7+05 

.I3060o+0l 


>1 


lU 

,05/5/2,05 

,abs.in/"00 




1 1 

. / t '/S/fl+05 



i ■ * 


12 

, /ai5o3+0b 

_ . ♦^2 36^0+ 11 1 


f \ 


li 

.643566+05 




Figure 3. 1-1. Sample Swath Reference Record 


KEPB.ODtJCIBILrrY OF THE 
OEIGINAL PAGE IS POOR 
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TEST CASt EOM HAMDoFf‘ 'STS , 


TtS I 1 o aO DAY SW 
Sh'ATH Table recori) 


40 day swath* U/ 50/75 START *2 VFHICL 
LE RtCttRU VFHTCI t I day 1 


tPP STMUUATtON 


65 



1893 , 

64 

1 1 ta 


1913 , 


1 IjJJL- 



h? 

1167 / 


r 952 . 

ft 1 

It 7 36 


1 '/ 72 . 

• 6:0 



_LV. 9 , 2 ... 


11650 

l^vq. 

2011 . 


1 1 9 116 


2051 . 


12.UL... 

• ‘>S 

SE 1?21“ 

bO U\4ll 
40 jSiSo 


.iLflA.A.: 

t< n 7 f' . 

2117. 

2ti7, 



iiWf’. 

21 Y'l, 


„?.05,04_ 
2Ui'0 • 
20oV. 

2 I 'I / , 
_? I <!&»-_ 
?1«0. 
?20'j. 


u n .. 

_-12'54 1- 



-- 

47 

12444 



Uh 

1 2. 9 P'6 



4.5^ 

-_L?.3 2'b. 


—Reins*. 

1 4« 

l236o 



I 

1 2*1 (1 4 

1 ? .* 



^iJ2 

41 l?bV/ 
40 12/12 

in l?/7/ 
i7 1?H0U 

i:6 

iS 

5« 1?U02 

3,5 l?Y.l.,«- 

i..? 1?V"2 

31 



20 1 ^oO f 

2« I5a2'j 
27_..i3u'ii.. 

rlf. 

2S rU'VU 

^2 4 l.J.U:!'/.. 

25 i5i.iOV 

22 15110 

21. _l3l.lY , 

20 151?7 ■ 

10 15153 




23.J2.- 

2531 , 
25/1 . 

2.5.V.0:._ 

,■'4 10. 

2020. 


?37>4, 
PiflO. 
_2300.,._ 
24 U'. 
?>\ 32. 


„.2 '1 ■■( O,.. 24.'i 7 j.„ 

2'toO, ?47/. 

24!ifl, 340 d, 

„ 531 

2327. 2333. 

2'.,4 7, 2334. 


JLiilEiSl. 

I0u3. 

T02S, 

1063. 
101)3. 
_2.0O5,_. 
2021 . 
2040. 

_i >0o0. , 
20 / 0 . 
209fl. 

2137, 

2136 , 

2 1 0 4 . 
2215. 
_,2 2-3.5 
2232. 
22 / 1 . 
_.229J.»._ 
2510. 
2320. 

_2i.4a._ 

256fl. 

2567. 

_- 2 U 0 , 6 ._ 

2«2n. 
24 45. 
_2.4r4._. 
24«4. 
2505. 
_. 252 ?,_. 
2542. 
2561 . 



2303. 

cibOO, 


261 3. 

ati9^ 

5 

?o52,. 



2 6 M4 • 

2b52. 


^6 b4 * 

26 71 . 



J=69J.- 

Rh9 1*— 

270 3,. 

2/10. 

^7iy^ 


2729. 

^ 73 / 1 , 


_274:‘/.. 

a7b'^s.^ 

A. 

/l7ti 1 . 

2/60. 

ri7/S* 

^761 * 

2/fUl. 


_i>00n a. 

^60./fc 



F igare 

3 . 1 - 2 . s 


^ULl/lliU-Ll 

,343352+01 

.341301+01 

.:3.30oio+.ai_ 

,33764(>t01 
.3i7'lO4 + 0l 

, 5 3 « 6 4 3 ♦ y_l_ 

.5351»'2-*U1 

.331W02+01 

.53,0,403+0.1._ 

,520254+01 

,52ao73+01 

..5260.3 3 A.0.t„ 

.525riaOHU 

.52««53+01 

,.52566./ + fli_ 

.322020+01 

,522007+01 

.321i?0iftL_ 

.520270+01 

.5194S8+01 

,31+'663+,u!_ 

.•317692 + 01 
,3l /l«640l 

.3l6‘il-0,.+.0L_ 

.513/1 .4+01 
,313023+01 

..31«354+,0J_ 

.313600+01 

,5l5o6U+01 

.3l?43‘4.+.y.1_ 

,511622+01 

.311222+01 

.■3l0b54 + l).t_ 

.5 (003/+ 01 
,300401+01 

,50H</5n+.0Ji_ 

,'5oa3ao< 01 
,3o7b3H 01 

,50'/.3?l + 0l_ 

,30M)00+01 
, 5o+)2fi3 ♦ 0 1 
,30Sy.7/ + t!.l- 

,3o3276+01 
,5(i4/«2 + 0l 

,.504 205+O.L 

,303o09H)i 
.303351+01 
,.3y26S/ + 4L 


.335223+01 

.3334a/+01 

._,.33LB/i/.+oL_ 

.330333+01 

,32895a+ul 

_.522B.03+.yl 

,526342+0 1 
.52514/+U1 

„,J2«0 ta+0.1 

,52203o+Ol 
.5210'>a+0l 
__.,52na u.t,v.,V — 
,510963+01 
.510U3/+01 

.,514184 + 0.1 

.517343+01 

,516333+01 

_..5lS./.4y.+.0..t 

,314901+01 

,314253+01 

_,513342+..0.1 

.5l264t) + ol 
.3121 /3+or 

,314674+01 
,5in24/+ul . 

,,.5.^06 53+0-1 

.5000 3/ + 01 
.504461+01 
_.5y7iM7.,+-0-I_ 
.•3.OT314+0 1 
,506/62 + lil 

.5,v62I9toJ 

,5v56'33 + 01 
.3US160+01 

.5-04:r,;i3t!l.l 

.304133+01 
.505611+01 
_,3o3l 3b + Ul_ 
.502643+01 
.302162 +01 

,5ot 6n4,+ yl 

,•30121 itOl 
,300/44+01 

.3U023 1+0,V_ 

,290632+0! 
,24*ri6.-) + 0 I 
,29^91 / + 0l 


,327176+01 

,525/23+01 

. 3 21* 3 66 + .0.1 

.323090+01 
. 5216P/+01 
_.,.i2.n-/no + O.L_ 
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.298349+01 
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, 2934SO+ 01 
,.295o2/,t01.. 
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TFST CASE FOH HANncFF »!>VS. TEST I - 40 OAV S^A TH, I 7b STAR1 t£ VKHICL LPP SIMULATION CASE 


1 PACE 


i-5 PkEOICTUR OiSPLAV FOR DATE 12/30/75 

ve hi cle 1 ^ 


-<1 


a. 


SEMT-NA.Iur axis 

_L.c.ceji.T_Ki.cj;xs: 


= 7290.000 

.s ,.O.‘55:0 O-flil-'L 


LFKTitr) s 

IHCl-XNATIflN = 99 2 2U 

jCiRlXiEMEAlJ- U.F..P.E Rl.G,fc£g-3/l.9 ,.52 ? 0 


RT. ASC. CKFENMCH = 9/ 48 41 


JiU. 




( 

‘ 

H 

?) 

M 

995? . bill 
1 <v; 53. ObTv 
2.2353j^iily 

201 30 59 
255 40 51 
229 So 43 


j 


fn 

!U00 

204 0 35 


i 


•ij 


I /8 10 27 


u ' 



'■V.'’55_._1.4'sij 

Vb2 ?0 39 


1 


7) 

4/i5S.^Jo7u 

126 3u 11 





5 i 3 5 A . 1 0 ’■> U 

loo 4o 3 


1 

• ■ 

- . 9) 

tj 5 A tj 1 2 0 

.. 74 /I9 5S 




in j 

65757. P340~ 

b9 *i7 


1 


in 

71957.7570 

V 39 





-Xe.i ‘j s.,.2 s.Qlu,_ 

3b7 19 3t 


1 



ri435n.flo?u 

331 ?9 





% 




Figure 3, 1-3. Sample Ephemeris Display 
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■■■ 


TEST CASt FOH HA#^r>oFT • SYS « TEST i - 40 OaY $1^ A 7 H ^ 1 cJ/ 30/ 75 START?ii VtHlCL SlttULA'llOM CASE 1 PAGE X2i 

■ TA^It: SEAKCH KFcm<D VEHICLE 1 


i\L n^ K>. 

,i0a0^«+0l 

__f JO / f\tSLl- 

,30/ / 4Si t>l 
,3.nA.i05%^0l 

,3iUfW0i0l 

.ii ifur>fu 

.31 3hO'5+fti 
,314431+01 
,315?/j7 + Ol 
.315993+Oi 
,31t)733 + 01 


“ 

so 

.Jf!.lai + ni 


49 



. 

.. ,,41V«pS + 0L 


4/ 

.i20OV9+Ol 

• 

4it) 

.3PUVIS+01 



.,4?.l,5.l?tfl.L 


4:4 

,3?lB9At0l 


4 3 

.iPcfOUvt n 1 




' 

4 1 

.i(^iSci<5+ni 

’I 

40 

,324V03i0l 

1 ■ 

39 

.3?4 6.p.A.+ n 1 


30 

,iPbt3ntnl 

.1 

3/ 



3o 

*ipL.l 3Bt0 l_ 

ii 

3 b 

, 3Po6 «;/■>+ 0 1 

■\ 

34 

, i? / 1 Uh + fl 1 

t 

33 



3c; 

,3P^bO=4n4nl 

"■j 

31 

,3p!.vb07 + 0l 


30 

, ypw9o 1 + 0 1 


;?9 


,1^ 

?b 

, n i 


?7 

.i5’W?9?+nl 

U_ - 



j 

?b 


‘I 

?4 

.i'<i5t«inl 

*r 

' ?3 

,332007+oi 

j 


,33.2425 + 01 

.i 

PI 

,332B4p.+ n L 

-i 

‘ 20 

,333Pbs’+Oi 


19 

, i33BoA+ 0 1 


IB 

, 3 340 74 til i 




.^dl\A7VMt 

, .jAb7b/Jtr 
,^<?.0So6 + ^ 
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”*;d9b?6/j+f 
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,3133/1 +( 
3 1 0 04 O-tJ 

■ .3146941 r 
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Figure 3. 1-4. Sample Table Seareh Record 
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3, 1, 5 File Requirements 
Input: None 

Output: SWATH - Swath Table File 

SWATHR - Swath Reference File 


3. 1, 6 grror and Recpyery 
3. .1. 6,1 General 

The program, will attempt to find as many sources of error during 
the input card processing as possible. The program will continue check- 
ing for further input errors upon detecting any input error. Since most 
of the computations are contained in subroutines which already exist, the 
philosophy of continuing proGessing after detection will be retained. There 
are two levels of error. These are: 

Level 1 - continue processing 

Level 2 - job fatal 

When a level 1 error occurs, the program will print an informative 
message and continue. If such an error occurs during the calculation 
phase, a printout of key data quantities is given in addition. When a 
level 2 error occurs, the program will print an informative message and 
return control back to the computer system. 

3, 1,6,2 Input Errors 

Level 2 

1. A check is made to see if NOSAY is between 1 and 549. 

Message: 

;C«*NOPAY IS OUT OF RANGE. IT IS NOT BETWEEN 1 
AND 549 

2. A check is made to see if INLAT and ISLAT are each between 

0 and 65. In addition, ISLAT(l) {INLAT(l)) must be ^ ISLAT(2), 
{INLAT(2)). Message: 

**«INLAT or ISLA T is OU T of RANGE, THEY ARE NOT 
between 0 AND 65 OR THE FIRST VALUE OF EITHER 
ONE IS NOT LESS THAN OR -■ TO THE, SECOND. 

3. A check is made to make sure the number of latitude points 
is less than 101 as specified by ISLAT and INLAT, Message: 


28E34-e>025-RU-OO 
Page 66 


***THE TOTAL NUMBER OF LATITUDE POINTS 
SPECIFIED BY ISLAT AND INLAT EXCEEDS. 100. 

4. A check is made to see if each entry pair for IPPI5 has 
entries between 1 and 549 and Uie first entry of a pair is 
^ to the second. Message: 

:}:*#IPPI5 (N) AND IPPI5{N+1) ARE OUT OF RANGE. THEY 
ARE NOT BETWEEN 1 AND 549 OR THE VALUE OF 
THE FIRST ENTRY IN A PAIR IS NOT LESS THAN OR = 

TO THE SECOND. 

5. IVEH is checked to be between 1 and 2. No more than 

2 vehicles can be processed by this program. Message: 

*«*IVEH IS NOT BETWEEN 1 AND 2. 

6. A check is made to see if each SA entry is in the range of 
-10 to +10 and is in ascending order and unequal to each 
other. Message: 

=S'**EITHER AN SA ENTRY IS NOT BETWEEN 0 AND 

ABS(IO) OR ENTRY(N) IS NOT LESS THAN ENTRY(N+1) 

7. A check is made to make sure IV1TIM(1) or IV2TIM(1) is greater 
than 1963 . Message: 

4c*5i«ErTHER.IVlTIM{l) OR IV2TIM{1) IS NOT GREATER THAN 1963. 

8. ORBITl(l) and ORBIT2(l)-a are checked to be in the range of 
6650 and 7700 kilometers. Message: 

**#EITHER ORBITl(l) OR ORBlT2(l) IS NOT BETWEEN 
6650 AND 7700 KILOMETERS. 

9 . A check is made to make certain that ORBITl(2) and 
ORBIT2(2)-e is in the range . 00001 and . 15. Message: 

:«cH:*EITHER ORBIT 1(2) OR ORBIT2(2) IS NOT BETWEEN 
.00001 AND . 15. 
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10. A eheck is made to make sure that ORBIT 1(4) and 
ORBIT2{4) -Qis in the range 0 and 360, Message: 

*=S'^EITHER ORBIT 1(4) OR ORBlT2(4) IS NOT BETWEEN 
0 AND 360. 

11. A check is made to make sure that ORBITl(5) and 
ORBIT2(5)-iO is in the range 0 and 360, Message: 

***EITHER ORBIT 1(5) OR ORBIT2(S) IS NOT BETWEEN 
0 AND 360. 

12. A check is made to make sure that ORBITl(6) and 
0RBIT2(6)-T is in the range +65. Message: 

❖::==HEITHER 0RBIT1(6) OR ORBIT2(6) IS NOT BETWEEN 
-65 AND +65. 

13. A check is made to make sure IVlTiM(4) or IV2TIM(4) 
specifies an hour such that the local vehicle passage time 
is between 700 to 1700 hours. Message; 

*4^'iVlTIM(4) ORIV2TIM(4) SPECIFIES A BOCAL TEHIC EE 
PASSAGE TIME NOT BETWEEN 700 AND 1700 HOURS. 



14. A check is made to make sure IV1TIM(3) or IV2TIM(3) is ;] 

i; 

in the range 1 and 31, Message: I 

*««iviTIM(3) OR IV2TIM{3) IS NOT BETWEEN 1 AND 31. • j 

15. A check is . made to make sure IV1TIM(,2) or IV2TIM(2) is i 

in the range 1 and 12, Message: 

sH:?::»IViTIM(2) OR IT2TiM(2) IS NOT BETWEEN 1 AND 12. 

16. If one or more of the above errors occur, the following 
nhessage will print: 

*** THIS JOB IS ABANDONED DUE TO THE FACT THAT 
I OR MORE FATAL INPUT ERRORS OCCURRED. 
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3, 1, 6. 3 Processing Errors ^ 

Level 1 

1. A check is made in subroutine HECTOR after attempting 
to solve for the eccentricity anomaly that the iteration 
loop converged on a solution. Message: 

=i=*=!'ERROR (KEPLER) SOLUTION FOR ECC. ANOMALY DID 
NOT CONVERGE AFTER 50 ITERATIONS. 

2. A check is made in subroutine REV TAB to determine if an 
anomaly oceurred in calculating crosstrack latitudes. If 
so« the following message is printed followed by data: 

REVTAB ERROR - ANOMALY IN CROSSTRACK LATITUDES 

Values for K, L, AA, JJ, DEL, XI, INT, NTRY.andLLL 
are then printed. 

3. A check is made in subroutine SWATH to make sure the 
iteration count of 5 is not exceeded in computing geocentric 
latitude and delta longitude. Message: 

Iteration limit on SWATH. Values for XLONGN, TEM, DXLON 
R(l), R{2), IFLG are then printed. 

3, 1, 6. 4 Innut/ Output Errors 

For sequential I/O the FORTRAN system on the UNIVAC takes 
control and prints a message identifying the problem and will either 
continue processing or abandon the job. If processing continues, the system 
counts the number of times this error re-occurs and if it happens a certain 
number of times, the system will abandon the job. 

For direct access I/O, the UNIVAC D.A, I/O package prints 
out an informative message, sets an error flag and allows processing to 
continue. In this program the swath table file is the only direct access file. 
Immediately after the informative message, the following message will 
print: 

^’i'^AN IRRECOVERABLE I/O ERROR HAS OCCURRED ON WRITING 
A RECORD TO THE SWATH TABLE. THE JOB IS BEING 

abandoned. 
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3 . 2 GRID 

3.2. 1 Inpait Card Data 

3. 2, I. 1 Inpat Data Deseription 

MINDEX - maximum number o£ index points to process on 
the NASA climatology tape. Range 1 to 16000, 
nominal value is 16000 (used to reduce run time 
for debugging cases). 

3. 2. 1.2 Card Format 


C.C. J 

1 5 


MINDEX 

5 

1 


• . * _ . 1 . . 1 . -1 


15 


Sample Card Inputs 



None. 

3. 2. 3 Output Report Data Definitions 

There is no normal printed output for this utility program. 

3. 2, 4 Samiple Pages from Each Output Report 
None. 

3. 2. 5 File Requirements 

Input: WEATAP - NASA Climatology Tape 

Output: INDMAT - Index Grid Matrix File 

3, 2, 6 Error and Recovery 

3. 2. 6. 1 General 

There are only three cheeks made for errors, one of which is fatal. 
3. 2. 6. 2 Error Checks and Mes sa ge g 
Fatal 

1, A cheek is made to make sure MINDEX on the input card 
Is between 0 and 16000. If not, the following message 
prints and the program stops. 
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Message: 

MINDEX IS NOT BETWEEN 0 AND 16000. PROGRAM IS 
TERMINA TED. 


Non-Fatal 

2. A check is made to make sure the index no. read from 
WEATAP file is between 1 and 16000. If not, the index 
no, is not incladed in the carrent INDMAT file record 
and the program reads the next WEATAP record. 

Message: 

THE INDEX NO. READ FROM WEATAP IS NOT BETWEEN 
1 AND 16000. THE DATA IS IGNORED. 

3, A cheek is made to see if an IGRID, JGRID grid pair has 
already been assigned an index nomber. If so, a message 
is printed and the current index no. replaces the one 
originally stored. 

Message: 

THE GRID PAIR IGRID nnn JGRID nnn HAS AN INDEX 
NUMBER ALREADY ASSIGNED, 
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3.3 LUMP 

Operational Assamtitiems 

• Only one case is run at a time, 

• Data Input Cards for the Sample Segments are sorted by 
country, region, zone, strata, substrata, segment, 

• Data Input Cards for the Substrata and Crop Calendar 'data 
are sorted by country, region, zone, strata, substrata. 

• Program is a one pass proeess in which cards are cheeked, 
written to disk if no fatal errors, and then LUMP is executed. 

• Maximum of 2000 training segments total. 

• The input cards for Substrata Historical data, Substrata 
StatistiGal data. Crop Calendar data, Crop Calendar Error 
data, and Sample Segment data are stored on separate input 
files. 

• The United States must be assigned the symbol 'USA ' and 
Canada must be assigned the symbol 'CAND'. 

3. 3. 1 Input Card Data 

3, 3, 1, I Input Data Description 

See Table 3. 3-1. 

3. 3. 1. 2 Curd F ormats 

See Tables 3. 3-2 and 3. 3-3. 


NAME 


Table 3. 3-1. LUMP Input Data Description 


SYMBOL 


NOMINAL 

VALUE 

RANGE 


i 

1 


0-9999 


1 


0-9999 


1 




1 


0-999 


1 


0-999 


1 


0-9999 


1 


0-9999 


2 





• 

0-999999999 

PW 

1 


0-100 


1 

• 

s, w 

fiPW 

1 


-9.999 to . 




9.999 

CV^ 

1 

* 

0-9.999 

CVg 



0-9.999 




0-9.999 

«PM 



9*.999 

j . 




0-3 

REPROD 


!* THE 


ORIGINS 


OOE 





±65* 


BH 








UNITS 


DESCRIPTION 


acres or 
hectares 


Deg, Min 


Identification case number for sample segment data 
(15) 


Segment number (14) 

Four character country name (i.e., USA, USSR) 
Region number (13) 

Zone number (13) 

Strata number (14) 

Substrata number (14) 

Substrata name (2A4) 

Substrata total area (USA and Canada acres, USSR 
and others-hectares) (floating point, FIO.O) 

% wheat in substrata (floating point, F7.3) 

S * Spring, 


Spring or winter wheat indicator 
W = Winter 


Deviation to true proportion wheat (floating point 
F6.3) 


Coefficient of variation for year-to-year change iri 
proportion wheat (floating point, F5.3) 


Coefficient of variation for within county variaticf 
of proportion wheat (floating point, F5.3) * 




ti 

JD 00 

fl) ^ 
-J I 

to 

o 

to 

\Jl 


Coefficient of variation for within county variatic 
of proportion mixed pixeis {floating point, F5.3) ' 

Deviation of true proportion mixed pixels (floatintl 
point, F5.3) 

Card type (li) 

=0 Substrata Historical Data 
=1 Substrata Statistical Data 
=2 Substrata Crop Calendar Data 
=3 Substrata Crop Calendar Error Data 
=4 Sample Segment Location Data 

Latitude with N or S indicating + or - 
respectively (Al, 12, 13) 


fo 

c 


O 

O 







Table 3. 3-1. LUMP Input Data Description (cont'd) 


NAME 


LONG 

IMOT 


IDAT 
lYRT 
■’ RSEED 


ILIST 




: ISHD 
ISSD 


ISID 


. 1 1 sew 

IlCWE 


ITSFG 


• NAGR ‘ 
NAL 
!ev4 
IGRP 
DELTWE 


SYMBOL 


DIMENSION 


IPLNTE 


6i-«4 


1 


8 


8 

8 

1 


1 

1 

1 

1 

4 


NOMINAL 

VALUE 


RANGE 


UNITS 


+180* 


1-12 


1-31 

64-99 

1.0-NNNNNNNN 

NN 


0.1 


0 , 1 


1-9999 

1-9999 

0- 9.999 

1- 3 

+ 1.0 


+99 


Deg, Min 


Days 


DESCRIPTION 


Longitude with E or W indicating + or 1 respective 
(Al, 12, 13) 


Date month, day, year of each of 4 phases for star^ 
phase and stop phase 4(1 X, 612) for winter or spri 
wheat 


Floating point odd whole number used as seed to 
random number generator (D12.0) 


List option to list all input data cards or only 
those in error (12) 

=0 Only those in error 
=1 List all input data cards 

Unit number of Substrata Historical Data (12) 

Unit number for Substrata Statistical Data (12) 

Unit number for Segment ID Data (12) 

Unit number for Segment Crop Window Data (12) 

Unit number for Crop Calendar Error File (12) - Se 
to 38 (8) 


I 


•da 


a 


Training Segment Flag (12) ^ 

=0 If all segments are training segments 
=1 If the training segment list is specified, ^ f | 
via segment location data cards. 

Number of agricultural segments in a substrata u 

Number of allocated segments in a substrata 

Multi -yr*Std. deviation of historical wheat area 

Group no assignment for a substrata 

Delta error for winter in predicting the bio windowl 
dates, expressed at fraction of the window interval 


eg 


Error in the mean planting date of the substrata 
relative to the true mean value. 











mm£ i UNITS I DESCRIPTION 


0-99 Days Within substrata variation for winter in crop 

calendar due to change in latitude, altitude, 
planting date, etc. 

0-1 0- training segment, ordinary segment 


0-9999 


S or W 
S or W 
ii • 0 

+99 

0-99 



List of segment numbers of training segments 
associated with the ordinary segment. These , 
segments are to be listed in order of decreas’ing 
priority. Highest priority is first. 

|Spring or winter v/heat indicator from 

|Crop calendar data 

Del ta error for spri ng - same as DELTWE 
Same as IPLNTE only for spring 
Same as for ISEGSD only for spring 
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Table 3. 3-2. LUMP Control Card 


1 9 21 26 31 36 41 46 51 56 


Case 

No. 

% 

Ran- 

dom 

No. 

Seed 

I 

List 
Opti on 


rrain- 

ing 

Seg 

Flaa 

IB 

Y 

Unit 

for 

Subst 

Stat 

)( 

Unit 

for 

Seg 

ID 

s 

Unit 

for 

Crop 

Calen. 

1 

Unit 

for 

Crop 

Error 

15 

D12.0 

12 


12 

^ ■ - — - 

12 

12 

12 

12 

12 


^S^f<^lLny 
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3. 3. 1.3 Deck Setup 

The order of the input is; 

1. Control Cards 

2. Substrata Historical and Statistical Data Cards 

3. Crop Window data cards and Crop Window Error Data 

4. Sample Segment ID Cards 

3, 3, 1, 4 Rules for Entering Data on Cards 
3. 3, 1, 4. 1 General 

• Integers must be right justified. 

• F-format or D-format numbers must have the decimal point 
present, 

• There is a crop window card necessary for spring and winter. 
If either spring or winter are not required, only one card need 
be input. 

3, 3, 1. 4. 2 Rules for Specific Fields 


See Section 3. 3. I. I. 
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3. 3. 2 Sample Card Inpats 


*! 

>1 

‘ f— 

lOi 

iy" 

I,.' 


LISTING OF INPUT CAWD IMAGES 


CONfPOL CARD 

1 .lonoooooooootOQi 


^ 10 11 12 6 


3. 3, 3 Outpat Report Definitions 

There are no normal report data items, 

3, 3, 4 Sample Pages from Each Qiatpat Report 

The following pages are samples of a card image printout and 
file contents dumps as controlled via ILIST, 



i 


SUaSTRATA historical AAP statistical OATa 


USA 


V 

IIS A 

i 

4 


lApACHte 0 o:o ooo ,_ao , ooo 

5*? 


USA 

1 

A 

^0 

1 ’.C50 .iso .lOU ,100 .bOO 

.100 


USA 

1 

a 

20 

SCOCONINO lOOOOOO. 20,000 



l-JiSA 

X 

4. 

20 _ 

5 fOSO si 00 ,100_,,100 ,&oo. 

,100 

\* 

IJ^SA 

i 

.4 


ISMoHAVt looioooo. aiiO.OOO 

39 

u 

USA 

i 

4 

20 

IS ,050 ,100 .100 .100 .soo 

0 10 0 

jC 

USA 

1 

4 


_i 7A* A V A J q 1 00 g o go . _ 2 o , o g o _ 

39 

f 

.U5A 

1 

4’ 

" 20 

17 '.o::6 ,ioo .ioo ,ioo ,soo 

,100 

s 

■USA 

1 

A 

20 

aSYAVAPAl lOOOOOO. 20.000 

39 

•t 

IJ^A ■ 

1 

4 

20 

2S .050 ,100 .100 ,100 .500 

,100 

IC 

' mA 

1 

4 

SO 

* 15MaH1c'upa locionob, so, boo 

39^ 


i|02_ 


A 02 


so 


15 ,0S0 


100 .100 ,J00 .SOtO .10 0 

20.0 0 0 ...SO A‘01_ 


1 

1 V rt 

USA 

1 

n. 

4 

ub 

21 .050 . 

too ,jO'i) ;ioo ,soo 

■'^ibb 


■ USA 

1 

4 

70 

2/YuUA 

10,00000* 20.000 

39 

401 


USA 

1„ 

4 

70 

27 .050 , 

100 .100 .100 ,500 

,100 

‘h- 

U^S^A 

i 

4 

90 

5C0CUTSE 

I'bO'OOO'o. 20,0 00 

39 

401 

IT 

U'SA 

X 

4 

90 

3 .050 . 

loo .100 ,100 ;soo 

.100 


USA 

1 

4 

9 0 

7GXLA 

lOOOOOO, 20.000 

39 

401 


U5-A 

... _ 

■ 4 

90 

7 oOSO , 

100 ,100 .ioo .500 

. 1 00 

f- 

•J^SA 

1 

4 

90 

<? graham 

iuooooo. 20,000 

39 

402 

t. 

'USA 


4 

90 

9 .050 , 

loo ,100 .100 .500 

,1 QO 

u 

V 

1 

4 

90 

UGf?ft;ULE£ 

ibb'oobb. 2b',ooo 

39 

403 


USA 

1 

4 

90 

11 ,050 , 

loo .100 .too .500 

• IOO 

hr 

• USA 

1 

4 

90 

19PIMA 

^ 1000000. 20.000 
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3. 3.5 File Reqairemenfes 


Input: SmSTD 

SSTATD 
SEGIDD 
SEGCWD 
CRPERR 
INDMAT 


Substrata Historical Data Cards 
Substrata Statistical Data Cards 
Sample Segment Location Data Cards 
Sample Substrata Crop Calendar Data Cards 
Crop Calendar Error Data Cards 
Index Matrbc Location File 


Output: SEGID 

SUBHST 


Crop Calendar (Window) File 
Substrata Historical File 


Inter me dia te : 

CRPINT 

SDBINT 

SEGINT 


Crop Calendar Intermediate File 
Substrata Intermediate File 
Segment Intermediate File 


3, 3, 6 Error and Recovery 


3. 3. 6. 1 General 


. The program will attempt to find as many sources of error during Input 
card processing as possible. The program will continue checking for further 
input errors upon detecting any input error. Proeessing will continue it 
possible if recovery from errors may be overcome. 

3. 3. 6, 2 Input Data Error s 

0 Compatibility checks (Do the country, region zone, strata, and sub- 
strata agree between the substrata statistical data and between the 
sample segment ID data and the crop window data?) 

• Checks for non-overlapping windows 

t Sequence checking of the card data 

• Input/output data l imit checks Crable 3. 3-4) 

• Training pfiority segment NNNN is either not a training segment or 
nonexistant-fatal. 

• A Group 1 substrata does not have at least 1 sample segment. NonfataT 

• A strata with at least 1 group 2 substrata does not have a sample segment 
itn any group 2 substrata . Nonfatal 

• A group 3 substrata has 1 or more sample segments . Fatal 
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Table 3. 3*4. LUMP Input Data Limit Checks 

t 


Parameters 

Range 

Substrata Area 

> 0 

Historical 'W 

0 to 10"^^ 

Deviation of True 

-9.999 to 9.999 

Yr-to-Yr CV of True PW 

0 to 9.999 

Within County CV of True PW 

G to 9.999 

Within County CV of True PM 

0 to 9,999 

Multi Year CV of Standard Dey 

0 to 9,999 

Ratio of PM to PW 

0 to 9,999 

Spring/Winter Wheat Indicator 

S or W 

Segment Latitude 

-65® to +65® 

Crop Vlindow Dates 


Year 

64 - 99 

Month 

1 - 12 

Day 

1-31 

Grid Index 

1 - 16000 

DELTWE , DELWEI 

*1.0, +1.0 
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• Region xxx zone xxx strata xxxx substrata xxxx beta distribution 
function could not converge on a solution 

t Region xxx zone xxx strata xxxx substrata xxxx beta function reset 

— — value since not between 0 and 1. Either or — have 

bad values 

• Region xxx zone xxx strata xxxx substrata xxxx seyment xxxx. No crop 
calendar data available for segment Intermediate file record, record 
skipped 

• Region xxx zone xxx strata xxxx substrata xxxx segment xxxx. No weather 
grid index no. could be found. The segment record is skipped. 

3. 3. 6, 3 Error Messages 
Card Field Errors 

For any error on a card the card image is printed prior to the message; 

1. Inconsistent substrata historical and statistical data 

2. Inconsistent Crop Window Error Data and Crop Window Data 

3. Substrata area must be positive 

4. Historical proportion of wheat must be non-negative 

5. DEVTPW must be less than 9.999 in magnitude 

6. COEFF, if variations must be between 0.0 and 9.999‘ 

7. DEVTPM must be between 0.0 and 9.999 

8. Latitude must be less than or equal to 65.0 deg, 

9. Longitude must be between -180 and 180 deg. 

10. Spring/winter wheat indicator must be W or S 

11. GV4 must be between 0.0 and 9.999 

12. Year must be greater than 64 

13. Month must be between 1 and 12 

14. Day must be between 1 and 31 

15. Overlapping crop window 

16. Group no must be equal to 1 , 2 or 3 

17. Delta errors in bio window prediction dates are not betv/een -l 9 iiid. +l 

18. Segment specified as ordinary must have a training segment priority 
list. 

19. Neither NAGR or NAL can be zero. 
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3.4 SEE 

Operational Assum'ptions 


• One case may be run at one time, 

• Only one country may be processed in a case. 

• The three data card sets must be on separate input files and 
in sort as follows: 

The major sort field is listed first. 

1. YES error data, set 5 

Country, region, zone, strata and card sequence number 

2. CAMS error data set 6 ' 

Country, region, zone, strata, substrata, segment and 
card sequence number 

3. Signature extension set 7 

Country, region, zone and card sequence number 

• All cards in a group with the same ID must be present even 
though a card may contain only blank data entries. 

• The program will require less than 20, 000 words decimal of 
storage in the CPU of the UNIVAC. 

o The program will run only if the specified substrata historical 
file is present. 

• Except in the case of an input card set being out of sort, the 
program will continue ehecking for field errors and for mis - 
matches vs. the substrata historical file, A fatal error will 
stop output files from being written. 

• Input data groups for which there is no ID match on the sub- 
strata historical file are considered to be extra data and are 
not written out on the output files. 
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• IDFRS and IDTOS specify ID limits for which certain errors 
are considered to be fatal. If a normally fatal error occurs 
on a card group outside these limits, it will be considered to 
be non-fatal. 

• A negative non-zero entry in any card input unit designator 
lYESR,. ISIGEX or ICAMER indicates that that data is not to 
be read or processed and the corresponding output file will 
not be created, 

• The United States must be assigned the symbol "USA " and 
Canada must be assigned the symbol "CAND".- 

• GAMS error data and signature extension data cards do not 
have to be entered at all levels. But there must be a YES 
data set for each strata entered since the yield value is unique 
for each strata. 

3.4.1 Input Card Data 

3.4, 1.1 Input Data Description 

See Table 3.4-1, 

3.4, 1.2 Card F ormats 

See Tables 3.4-2 and 3.4-3. 




the Table 3«4«1, SEE Input Data Description 


pa(;rjspV,k 


Dimension 



Nominal 


Value 


0,0 


0,0 


1 

2 


Range 


1-9999 


0-2 


Units 


Description 


Case no. to be assigned to identify the three 
output files. 


List option to list all input data cards or only 
those in error and all output file records. 

= 0 list cards in error 

- 1 list all input data cards 

- 2 list all input data cards and all output 

file records 


Specifies ID of ending zone for which errors 
are considered to be fatal. 1st item is i*eg., 
2nd - zone. 


File unit no. for YES input card data. Minus 
value means not to read or process this data. 


File unit no. for signature extension input 
card data. Minus value means not to reau 
or process this data. 


File unit no. for CAMS error input card 
data. Minus value means not to read or 
process. 


Specifies ID of starting zone for which errors 
are to be cGnsidered fatal. 1st item region, 
2nd - zone. Blank or 0 entry means entire 
ID range is active. 


Zero or blank means not to check on case no. 
of input substrata historical file. 


M 

OP 

M 

•tk 

OS 


§ 


Four eharacter country name (i. e* , USA 
USSR). 


(u 

o 


00 

nO 


3 

§ 

1'^ 

o 

o 




Table 3.4-1, SEE Inpvut Data Description (cGnt’d) 





Nominal 



• 

t 

Name 

Symbol 

Ditnei^sion 

Value 

Range 

Units 

Description 

IREG 

- 

1 

- 

0-999 

• 

Region no. 

IZONE 

- 

1 

- 

0-999 


Zone no. 

ISTRAT 

ISUB 

- 


1 

1 

- 

0-9999 

0-9999 


strata no. 

Substrata no. ti>. ' j 

Segment no. * a 

. ISEG 

- 

1 

- 

0-9999 

• 

YIEDDI 


1 


0-99.99 

Quin/ 

True yield of stratum (for USA , and CAND - 

■ 




hoc tar 

bushels /acre. 

: lERDTE 


6x3 

0 

See error 

Yr, Mon, 

Year, month and day for error truncation. 





checks 

Day 

H 1st value is 0, then there is no bias or std. 
dev. data. Any other 0 value terminates 
data entries. 


* BIASYI i 

®YI. 

X 

6 

0 . 

+99. 9 

Quin/ 

Bias error of yield. One value for each ^0 




he etar 

value of lERDTE for USA , CAND - bushels/ 







acre. 

SIGYIl 

°YI. 

1 

6 

0 

0-99.99 

Quin/ 

Standard deviation of yield error. One value I 





hectar 

for each ^0 value of lERDTE for USA , | 







CAND - bushels /acre. 

PWW 

P(W/W). 

4 

0 

0-1 

- 

Nominal frequency of labeling wheat given 
pure wheat pixels (conditional prob. ). One m 







for each window, ^ S ' 

a u» 

PWM 

P(W/M). 

4 

0 

0-1 

_ 

4). 

.Nominal frequency of labeling w'heat given ^ ^ “ 






mixed pixels (cond. prob. ). One for each o * 

window. m {: 

■ i 

I 

; PWO 

P{W/0), 

4 • 

0 

0-1 

. 

Nominal frequency of labeling wheat given G 


1 : 


« 

• 


pure other pixels (cond. prob.). One for o ^ 

each window, ® i 




• — ^ - r - 

‘ 



- 



a 




' ^ 




SEE Input Data Description (cont'd) 



• 


Nominal . 

1 1 




Name | 

Symbol | Dimension { Value | 

Range | Units I 

Peseription 




SIGWW 


SIGWM 


SIGWO 


SiGPW 


BIW 

B2W 

SIGIW 

SiGZW 


B(W/W) 
BIW /M) 


er(W/W) 

cr{W/M) 



+9.999 
£9.999 
+9.999 
0-9.99 
0-9. 99 
0-9. 99 
+9.999 
0-9.99 


Bias error in classifying pure wheat pixels. 
One for each winHow. 

Bias error in classifying mixed pixels. 

One for each window. 

Bias error in classifying other crop pixels. 
One for each window. 

Standard deviation in classifying wheat 
pixels. One for each window. 

Standard deviation in classifying mixed 
pixels. One for each window. 

Standard deviation in classifying other crop 
pixels. One for each window. 

Bias in proportion estimate (model #2). 

One for each window. 

Standard deviation in proportion estimate 
(model #2). 


FOR MODEL, 1 SIGNATURE EXTENSION 
+9.999 I - I Multip 


+9. 999 
0-9.99 

0-9. 99 


Multiplicative bias error W/W. 

Additive bias error W /W, 

Multiplicative std. dev. error (W/W). One 
for each training priority segment. 

Additive std. dev. error (W/W). One for 
each training priority segment. 
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Table 3.4-1. SEE Input Data Description (cent 'd) 


Name 

Symbol 

i ' 'i 

: BiM .j 

i I 

®IM 

B2M 

®2M 

V ^ SIGIM 

‘^liM 

' 'ii 

V ; SIG2M 

\' t ■ -i 

^21M 

1 ; 

: 1 BIO 

i 

>io 

,, | .B20 . 

®20 

1 1 SIGIO 

■ 

®UO 

|: SIG20 ' 

i 

t-''" 

‘^210 

h' B1 

% 

,Bz 


i?:.| $IGl 

'.'•I ’ 

% 

SIGZ 

^21 

ISET 

- 


Nonainal 

Dimension I Value 


Range 

±9* 

999 

i9. 

999 

0-9 

. 99 

0-9 

. 99 . 

+9. 

999 

+9. 

999 

0-9 

. 99 

0-9 

.99 

D E L, : 

?. SI G] 

+9. 

999 

+9. 

999 

0-9 

.99 

0-9 

. 99 


Units 


Description 

Multiplicative bias error (W /M), 
Additive bias error {Yv'/Mi, 


Multiplicative std. dev. error (W /M) , One 
for each training priority segment. 


Additive std. dev. error (W/M). One for 
each training priority segment. 


Multiplicative bias error (W /O) , 
Additive bias error {\V /O). 


Multiplicative std. dev. error (W /O). One 
for each training priority segment. 


Additive std. dev. error (W /O), One for 
each training priority segment.' 


Multiplicative bias error 


Additive bias error 


Multiplicative std. dev. error.. One for 
each training priority segment. * 


Additive std. dev, error. One for each 
training priority segment. 


Card set number in column 80. 5 for YES 

error data. 6 for CAMS error data. 7 for 
signature extension data. 
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Table 3,4-1, SEE Inpat Data Deiscription (eont’d) 


Nominal 



Hange 


Units 


Card sequence nurmber witiiin card set with 
the same ID. Used to identify each card of 
a group. YES error has a 2 card group. 
CAMS error and sig, ext. each have 4 card 
groups. 
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Table 3. 4 - 2. SEE Control Card 


1 

5 

6 

9 

12 

IS 

18 

20 

22 

24 



IDFRS 

IDTOS 


« 

o 

ICAMER 


ICASE 

H 

10 

3: 

M 

Region 

Zone 

Region 

Zone 

lYESU 

ICSESH 

14 II 

13 

13 

13 

13 

12 

12 

12 

14 




3 ’ 










A 
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Q 12 15 21 27 


jRdga jZone 


13 F6.3 



35 39 43 


Mult. Std, Dev. Error W /W 


Signature Extension Bata Card 1 
39 43 47 51 55 


°‘l2W ®13W *^14W 




1 AdditiV^e Std* Dev. 

Error W/W 

- 

*^21W 

^22W 

®^23W 

^24W 

®25W 

*^26V^ 



I 


6F4. 2 


6F4. 2 


15 21 27 



Signature Extensio^i Data Card 
35 39 #3 '47 51 


Mult. Std., Dev. Error Vf /M . 


55 69 63 67 


Additive Std. Dev. Error W /M 


^ ®11M 


®^13M *^14M 




®2 IM : *^22M ®23M *^24M ®^25M °26M’ 


■ 


F6.3 F6.3 



6F4. 2 

15 21 27 

31 

f 

35 

Signature Extension Data Card 3 
39 43 47 51 



Reg, Zone 


A4 13 13 F63 F6.3 


R 




Mult. Std. Dev. 

. Er ror W /Q 


®20 

; °ll!0 

*^120 °13'0 

*^140 

‘^150 

^160 


•k dditive ,St 


®210 '^220 


*^230 ®240 


‘^250 ! *^260 



6F4.2 


6F4.2 


9 12 15 21 27 


Reg. Zone Bj 


F6.3 Fa3 




Signature Illiitension Data Card 4 

3 i 35 39 43 47 51 55 59 ^3 ^ 

MaltipM-eatrv’a Std. Dev. I Additive Std. Q ev. Krror 














*0 

dT oo 

TO 

nO I 


6F4, 2 

Table 3.4^3. SEE Data Garde (conit'd) 


6F4. 2 


























































''REPRODUCIBILITY OF THE 
CltlGINAL PAGE IS POOR 


CiAMS Error Data Card 1 for Pure Wheat Pixels 


i 

] -5 

9V' 

12 

15 

1=9 23 

27 

31 

35 

39 45 51 57 ’ 

63 

67 

71 

75 

- — 1 

79J 

)0l 

i 



; 



NGminal Freci* Label Wht. 

Bias Error in Class ifvinir ’Vyiveat 

1 Std. Dey. Clas9l£tf„WJiL 



peig- 

hent 

Country; 

Reg.: 

Zone 

Strata 

Sub- 

strata 

p{wr/w}j. 

P(W/W)2 


PfW/W)^ 

n ^ I [ 

B(W/W)j| Bfwr/Wj^j BiW/Wl^jfiCWAV)^ 

aOV/V.Ui 

aCW/W)2 

tr(W/W)3 

ET(W/W)^ 

■ . 

1 

s 

Sl4 

^ A4 

13 

13 

14 

14 

4F4. 2 


4F6. 3 


4F4, 2 





CAMS Error Data Card Z for Mixed Pixela 


|eg- , 
jriant; 

■ 

Country! 

ba 

: 1 

Zone 

Strata 

Sub- : 
strata 

. Nominal Fre.q . Label Y/ht, 

Bias Error in. Classifying 

r Wheat 

Std, Dev. Classify Vfht. 

2 


P{\V/M)j 

PfW/M)^ 

P W/M)^ 

PfW/M^ 

;b (w/M) j 

B(W/M)2 

B(W/M)^ 


0:(W/M) ^ 

aCW/M)^ 

a(V//M)j 

cr{iy/M)^ 

ic4 A4 • 

s 

13 13 . 

9 12 

14 14 4F4.2 4F6. 3 4F4. 2 

CAMS Error Data Card 3 for Other Pixels 

15 .19 23 27 31 3S 39 45 51 57 63 67 71 75 79 80 

i 

;?eg“ 

i'Eieht 

!' 

Cfluntry 

Reg. 

Zone 

Strata 

Sub- 

strata 

Nominal Freq. Label Wht. 

Bias Error in CJasaifyinff Wheat , 

Std. Dev. Classify Wht. 

3 


PCw/o)j 

PCW/G!)^ 

P(W/0)3 

P(W/0)^ 

S{W/0)j 

B(W/0)2 

B(W/0>3 

BfW/O)^ 

CTCW/0)ji 

afW/O)^ 

aCW/0)2 

crr.\70)^ 




A4 


13 14 


14 


4F4* 2 


4F6, 3 


4F4.2 


GAMS Error Data Card 4 



c, 

9 12 

15 

19 

23 27 31 

35 

39 45 

5 1 

57 






Sub- 

sitra.ta 

Std. I 

)ev. Proportion Est, 

Bias in PronortiG-n Estimate 

?eg- 

jrient 

i\ ■ 

Country 

Reg. 

Zone 

Strata 

a(P/W) ^ 

s0(P/W)2 

ff(P/W)2 

crtP/Vf)^ 

i 

3(p/w)^|BtP/yr)2 

B(P/W)3 

^ J 

B(P/W)^ 


63 


A 4 


13 13 


14 


4F4. 2 4F6, 3 

Table 3.4-3. SEE Data Cards tcontM) 





14 
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3, 4, 1, 3 Deglt SetiiiP 

The input card data strearn will have one card, the program control 

card. 

The Y£S data will be on the file specified in lYES^J (IJnit 1 is assumed 
if this entry is 0). 

The CAMS error data will be on the file specified in ICAMER (Unit 4 
is assumed if this entry is 0). 

The signature extension data will be on the file specified in XSIGEX 
(Unit 2 is assumed if this entry is 0). 

In the case of the signature extension and CAMS error data cards, if 
the data is constant for a particular ID level such as zone, only one card 
group may be filled out for that zone and all ID levels below zone are left 
blank. , 

3. 4. 1. 4 Rule 8 for EnteriBg Data Item s on Cards 

3. 4. 1. 4. 1 General 

1. Integers must be right justified. 

2. Alphanumeric data must be left justified, 

3. F format numbers must have the decimal point present, i. e. , 

F6i 2 - NNN. NN. However, the user may override the specified 
field format as long as the total field width is not exceeded. 

4. Except for IYESU , ISIGEX and ICAMER, aU other nominal values 
are zero. 

3, 4. 1.4.2 Rules for Inputting Specific FLeids 

• ICSESH must mateh the ease no. of the substrata historical file. 

A If IDFR3 and IDTOS are left blank, then it wiU be assumed that the 

Occurrence of any fatal error will terminate the Writing of that 
output file with which it is associated. 
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If the first entry of lERDTE data is 0 or blank on card 1 of a YES 
data group, then the card 2 data is ignored (the card 2 must be 
present but can be blank except for ID), The lERDTE data and 
card 2 data from the last non<>zero lERDTE data group is used to 
obtain values for the current group. 

All standard devLationi bias and frequency data items are checked 
to be within valid ranges. See Section 3.4. 6. 2 for specific field 
cheeks. 






i 





I » -M tri J t-aa. 












3» 4, 2 Sajmple, Card. Data 
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IDFPS IDTOS 

I ICA$£ HIST KEb zone 0E6 ZONE 1Y£SU tSTGEX ICANR KSESH 

_i -2. > -fii ! ao 22 22 i_ 


3. 4, 3 Qatpat Report PefinifeiQpg 

There are no normal report data items, 

3,4,4 Sample Pages from Output Reports 

Tiie following page is a sample of a card image printout and 
file contents dumps as controlled by ILIST. 



V:£8 data liROUl* 

^ USA 1 a a 04 ^.gjft 75 t 2 m i ITO 2 t 7 ^> 3 176 31 ST 6 ms 

»B ,6 ,A .a .0 •.£ S.ob 4.00 3.00 2.00 1.00 

YtStRR FILE 

ICrBV a USA IRg R g I I?QRE s 4 18TRAT a 20 
Y = 6 .ftb CITRt>ffc(J).BIABVtij^sisriEnj«Iait 4 ) 

0465 .5 3.36 0496 .4 a,69 9S2T .3 2.02 9556 ,1 J,3b 95T0 .0 

— — . _ ___ ^ i , ^ 


iS_ 

.0025 


«OUiSA 1 a 

20 - 01,00 ,90 , 601.00 

1,00 «70 .50 .50 

,000 

,000 

,000 

.000 

,000 
, 0 U 0 

,000 

.000 

.10 

.10 

•to .10 
.10 .10 

.1016 

.1026 


1 , 0-0 . 5.0 .30 , 0'0 

. OiOO 

,000 

. n 00 

.000 

.To' 

.10 .10 

alOS 6 


,30 .ao ,10 ,00 

-.100 

•,050 

-.000 

-,000 



06 

I REG 1 iZONE 

4 ISTRaT 20 TSU 0 -0 

ISEG 

0 



• , 





tCiRT USA I'HgG I IZ 0 N£ 4 ISTRaT 20 TStlB ^O TSFG 0 

The camfhr input Data gpouk ddes not match with subhst file and is extra data 

— ^ ^ Y^OTT A ■ G R iilM " 

o 8 A a 4 SOiO.SoTSR 176 1 176 2 l 76 3 176 31576 415 
: ■ ■ . • .5 .3 ■■ 1 * . 1 . 3 ^.5 5.00 4. ,00 3.00 2 ei 00 l.flO 

v'esVrr pitF 

irtRY a USA IREG a 1 IZONE a 4 ISTRAT s 50 

Y a 7.06 HTR DT|;(H tBiASVtT l«STGTlEmrt al. 6 J a 

9465 .3 3,36 9496 .2 2,69 9527 ".1 ~ 2.02 9 SS 6 .,1 1.35 9570 -.2 


1.35 9570 


CAMS DATA GROUP 


mOt 

,00 

.90 

,601 

.00 

. 00 0 

, 0'00 

,000 

,000 

.10 

.10 

.10 

.1016 

1 

.00 

.70 

.50 

.50 

.00 0 

,000 

, 000 

,000 

,10 

e 10 

• 10 

.1026 

1 

.00 

.50 .30 

.00 

.000 

.000 

,004 

,000 

.10 

.10 

oio 

.1036 


, 3.0 

.20 

,10 

,o'o 

-rlblO 

-. 450 “ 

-‘roo’b 

-,ooo 




4 b 



ICTRY a USA iRfc'G a l I'ZflNt a 4 iSf«AT = 


GftMtRH 
'so IS^UB a 


13 IS EG a i 0 O 2 


CPWW( ri»'PKM(l)»PV:ntI) tOWRfl)rSIGWWCT 1 f 5 Wrtrn,SiGRHiCI) »0WOf z) »SlG«D(I),bPWtI) v.SIGPWII)tIal, 4 l a 

li ■ tf o 1 . 6 0 1 ,00 , OOP .IqO .OOP ,10 , 00 , 1 0 » , 100 , 30 ^ 

. O'l): . 7*0 . 5 ' 0 ‘ B O‘ 00 ‘ . 1 0 0 ^ 00 : 4 . .10 # 00 ' .10 ^.OSO # 2'0 

,ttO .50 ,30 , 00^0 , 1 oO ,O 0 'O .10 ,00 #10 -oUOO olO 


IGTkV a USA I REG = 


tZONE 


istrat 


gamep.r file 

50 ISOB a 


21 ISEG a 1066 


riP'-ii I ) Vpr® c I ) Fifirij ;s i Gvi w‘( I ) » ewHiriVVai g h m ( iT , a wo ri ) i si gr o rni 6 P rT 1 1 » sifepiifl IT t fa i ,.4) * 


1.00 

1 . 00 

1 .00 

.000 

, i 04 

. O'O 0 

.10 

.00 

.10 

-.100 

.10 

- VO . 

. 7 ^ 

. 5 i 0 

.OOiO 

. to 0 

.040 

.10 

.00 

.10 

- g 05 0 

.£0 

oOO 

• Sy . 

; 3 0 

.000 

.10 0 

.00 0 

.10 

.00 

10 " 

-aOflO 

Vio 

1.00 

.50 

.00 

,600 

,.loO 

.060 

olO 

.00 

• 10 

- .00 0 

.00 



YES ' 0 AT A GROUP " 

USA I 4 7010,707512 176 1 1/6 2 176 3 176 31576 415 ' 

,i -,l -,3 -.5 -,7 5,00 4,00 3,00 2.00 l.OO 


YESERR FILE 

I C TRY a USA iREG a 1 IZQ«E a 4 ISTRAT a 70 
7.20 tTTRuTE(l>, 8 lA 5 YCI)»SIGYXE(I}»Ia 1 , 6 ) a 

9465 .2 3.36 9696 .1 2,69 9527 -',1 2,0 2 9556 <-.2 1,35 9570 “.S 

. ^ \ ' . . . - , G AMS 041 A GROUP 


-.a 1,35 9570 
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3, 4, 5 File Requirements 


Input: YFSFCD 

CAMFCO 
SIGECD 
^BHST 


Yes Error Data Cards 
CAMS Error Data Cards 
Signature Extension Data Cards 
Substrata Historical File 


Output: CAMERR • 
SIGEXT - 
YESERR - 


CAMS Error Data File 
Signature Extension Data File 
YES Error Model File 


3.4,6 Error and Recovery 


With, the exception- of the wrong substrata historical file being mounted 
or any data out of sequencet the program will continue to check for all types 
of errors even though the writing of a given output file has stopped. There 
are three categories of errors; 

Fatal with immediate termination - level 3. 

Z. Fatal for a given input data type - error checking continues 
but the Writing of the corresponding output file stops - level 2 . 

3i. Nonfatal - normal processing continues - level 1. The card 
image of the card with field errors or card ID which triggered 
a general error will be printed prior to the message. The 
error checks, messages and fatal status are given in Sections 

3.4. 6.1 and 3.4. 6. 2 below. 

3 . 4. 6. 1 General 

Level 3 
Code No, 

A check is made that the correct substrata historical file has 
been mounted as speGified by IGSESH. 

Message: 

EITHER THE SUBSTRATA HISTORICAL FILE HAS NOT BEEN 
MOUNTED OR IT HAS THE WRONG CASE NO. PROCESSING 
IS TERMINA TED. 

2* A check on sequenGe is made from GOuntry down to the card 
sequence no. for each active card get, jn addition, a 

mi'ssing sequence no, in a card group with the same ID will 
cause the same action to be taken. 


I 
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Message: 

THE ... INPUT DATA SET IS OUT OF ORDER. PROCESSING 

IS TERMINA TED. 

3. A check is made to make sure the ID in IDTOS is ^ the ID 
in IDFRS. 

Message: 

ID IN IDTOS IS NOT GREATER OR = TO IDFRS 

The following errors are level 2 if the error occurs when the current 
ID is in the range specified by IDFRS and IDTOS. Otherwise, it is a 
level 1 error. 

10. A check is made to make sure that each SUBHST record has a cor. 
responding record for each of the active input card sets for YES and 
Sig. Ext. This situation will only be considered an error in case of 
Sig. Ext. input if there is at least one segment in the SUBHST record. 

Message: 

Optional line 

LEVEL 2 - WRITING OF FILE TERMINATED 

ICTRY ---- IREG IZONE — ISTRAT ---- ISUB ---- 

input data set DOES NOT HAVE A RECORD WHICH 

CORRESPONDS TO SUBHST ID. 

11. A check is made to make sure that for each record in SUBHST 
with one or more sample segments there is a CAMS input card 
group for each segment. 

Message: 

ICTRY IREG — - IZONE --- ISTRAT - — - ISUB - — - 

CAMS INPUT DATA SET DOES NOT EXIST FOR SEGMENT 

IN SUBHST RECORD. 

Level I Error 
Code No. 

15. A check is made to see if an extra input card group has been 
input with no match on the SUBHST file. 
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Message: 

THE INPUT GROUP DOES NOT MATCH WITH SUBHST 

FILE AND IS EXTRA DATA 


3, 4. 6, 2 Card Field Errors 

The occurrenGe of any of these errors is level 2 if the error occurs 
when the current ID is in the range specified by IDFRS a^nd IDTOS. Other- 
wise, it is considered to be a level 1 error. 


Code No. 

20. A check is made to make sure the error truncation dates in 
lERDTE on YES card 1 are in ascending order (if not all 0), 

Message: 

lERDTE ARRAY HAS DATES NOT IN ASCENDING ORDER 

21. A check is made to make sure YIELDI on YES card 1 is 
between .0 and 99. 99. 

Message: 

YIELDI IS NOT BETWEEN 0 AND 99. 99 

22. A check is made to make sure that BIASYI on YES card 2 
is between -99. 9 and +99. 9. 

Message: 

ABSOLUTE VALUE OF BIASYI IS NOT LESS THAN 99. 9 

23. A cheek is made to make sure that SIGYII on YES card 2 
is between 0 and +99. 9. 

Message: 

SIGYn IS NOT BETWEEN 0 AND 99. 9 

24. A check is made to make sure that the P{W /i) data for a 
CAMS card group is between 0 and i. 

Message: 

PW- IS NOT BETWEEN O AND 1 
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25. A check is made to make certain that B(W/i) and for a 

CAMS card group is between -9. 999 and +9. 999. 

Message: 

ABSOLUTE VALUE OF BW- OR BPW IS NOT LESS THAN 9. 999 


Code No. 

26. A cheek is made to make certain that for a CAMS card 
group is between 0 and 9. 99. 


Message: 


27 . 


Slow- OR SIGPW IS NOT BETWEEN 0 AND 
A check is made to make sure that B j^/and B 


9.99 

data on the 


signature extension card group is between -9. 999 and +9. 999. 
Message: 


ABSOLUTE VALUE OF Bl- OR B2- iS NOT LESS THAN 9. 999 

28, A check is made to be sure that cr or data on the signature 
extension card group is between 0 and 9. 99. 

Message: 

SIGl- OR SIG2- IS NOT BETWEEN 0 AND 9. 99 

29. A check is made to make sure the dates in lERDTE are in 
range. Year must be greater than 64, month- must be between 
1 and 12 and day must be between 1 and 31, 

Message: 

lERDTE ARRAY HAS A BAD YEAR, MONTH OR DAY NUMBER 






J-Lnr 


3.5 SAG£ 


- w 




I 
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OperatioBal As sumptions 

• A maximum of 426 days can be specified in NODAY--the number 
of days in the run. 


• It is assumed (hat only 1 year of data, starting January 1st, is 
to be generated on the weather data file (366 days). 

• The use of the weather file will be such that the look-up date will 
be determined modulo 1 year, i. e. , given a run day of 400, then 
(he look-up day will be 400 - 366 + 1 = 35, 

• Only 1 case is rtm at a time. 

• Only 1 country may be run at a time. 

• The segment reference data file and the weather file axe re- 
generated each time the program is run. There is no update 
capability. 


• A maximum of 150 acquisitions for any one segment is allowed. 


3. 5, 1 Input Card Data 


3. 5. 1. 1 Input Data Descriptipn 
See Table 3.5^1. 

3, 5. 1, 2 Card Formats 

"SAGE” is punGhed in card eolumns 75-78 of all eards. A sequence 
number is punched in card Golui^ins 79-80. 

See Table 3. 5-2 for card format. 


i- 





I 








REPRODXJGIBlIirrY OP THE Table 3. 5 -1. 

origmal page is poor 


Narrie 


Nominal 

Symbol Dimension Value 


ICSESW 


IGSESG 


NDAYS 


I lOPT 


I IVEH 


i iREPT 


DECR 


RAND 


IGRDN 




SAGE Input Data DesGription 


Range 


Units 


Description 



72 character case header which prints out 
at the top of every page. 

0-9999 - 4 digit case no. to identify the printed 

output and the segment reference data 
file. 

0-9999 - Case number identifying the swath table 

and reference files. 

0-9999 - Case number identifying the segment ID 

file. 

ISTIME(l) - (Year ^ 1900) 
iSTIME(2) - month no, 

ISTIME(3) - day no. 

Run start date 


1-426 


Days No. of days to process in the run. 

- Program run option, 

0 - run sage only 

1 - run weather file generation utility only 

2 - i^un both the utility and SAGE 

- List of vehicle number s to process in 
this run. 


T or F , - Flag to indicate whether access report 

is to be produced, T - yes, F - no. 

0- 100 Kilom. Swath deerement 

- - Random no, seed used to obtain daily 

weather data. Must be odd integer. 

1- 16000 - No. of grid points to be written on the 

weather file (debug only) 
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Figure 3.5-2. SAGF Data Card Formats 


Header Card 


T&e header ii^ormation is entered in C. G. 1-72; "SAGiE" is entered 
in C. C. 75-7'8. and 01 is entered in C. C. 79-80. 


5 h 


13 17 
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ta 
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ISTIME 
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o, 

h-l. 

0 

to 
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3» 5, 1. 3 Deck Setiip 

1« Header card . sequence 01 

2. Data card - sequence 02 
3, 5* 1. 4 Rules, for Entering Date, on gards 
3«5»1«4»1 General 

!• Integers tUusi be entered right justified. 

2. F format numbers must have the decimal point present* i. e. * 
F6.2 ‘ XXX. XX. 

3. The card sequence numbers in C.C, 79*80 must be present 
in all data cards. 

3* 5. 1.4. 2 Rules, for SnecifiC Field's. 

• 1CS£S\7 must 'match the case number on the swath reference 
and swath table files. 

• ICSESG must m.atch the case number on (he Segment ID file. 

• The start time in ISTlklE must not be less ihan the earliest 
vehicle start date on the swath table file. (Note, if only 1 Of 
the vehicles is to be processed, then ISTDvlE will be checked 
against that vehicle's start time" only. M ISTIME is not input, 
it will be assumed that the earliest vdtxicle start date is the run 
start date. 

o NDAYS must be a, I and ^ 426.' 

O MVEH must have entries of 0, 1 or 2. 




SA 6 F test case For HAiNi>Of F,SrSi test U» 3 S aAT»i/l/ 7 « START UP 

HEAUEHI 

SAUE TES'r CASE FOR H‘A«UOFF»Sr^^^^ TEST f#fs 0 AY»l/i/ 7 « SIAkT 


UPF SIMUbATION CASE 1 RARfc 


iCASt TOPr N 0 »YS lesESR iCSfSG ISTTMt LVEH IRCPT OECR 

■ 1 0* 3S' i 1 7r~l 1 i -0 1 STOo 


HAMW 

, I;000tt0o6u000)f Aoi 


IGROti ISWRTO NLA)Sh 
ifO *?T VT“ 



» 
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3. 5. 3 Qgfcpiit Report Dgi^ Definitioas 
See Table 3.5.3, 


3, 5,4 Sample Pages from Eaeh. Qgtpgt Report 

See Figure 3, S-l for Sample Segment Access Data Report. 


3, 5, 5 F Ue Requirements 


Input: SWATH 

SWATHR 
SEQID 
WEATAP 
WEATHR 


Swath Table 
Swath Reference 
Segment ID File ‘ 

NASA Climatology Tape (optional) 
Weather Data File (optional) 


Outpiut: SEGREF - 
WEATHR ^ 


Segment Reference File (optional) 
W eather Data File (optional) 


3. 5, 6 Error and Recgyery 
3. 5, 6. 1 General 


The program will attempt to find as many errors during the input 
card processing as possible. The program will continue checking for 
further input errors upon detecting any input error. There are 2 levels 


of error. These are; 


Xievel 1 . continue pro ees sing 
l^vel 2 - job fatal 

When a level I error occurs, the program will print an informative 
message and continue. When a level 2 error occurs, Sie program will 
print an informative message and return control back to jdie computer 
/System. 

3. 5. 6. 2 Input Data Errors 
.Eevel 2 

1. A check Is made to see if iCSESW matehes the case no. on ' 
the swath file Si Message: 

*«t«iCSES^V DOES NOT TCK THE CASE NO ON THE 

SWA TH INPUT FILES OR SWATH FILES NOT MOUNTED. 


' p’’* fc 








Report Napae 


COUNTRY 

Access count 

ACQUISIT 

MO/DA/YR 

REV 

VEH 

GTA 

EL 

CLOUD COVER 
ACQUISI. TIME 
ALTITUDE 





Segment 




Access 

- 

1-150 

-- 

if 

tr 

- 

1-18 

— 

it 

- 

1-2 

— 

M 


-10, +10 

Deg, min, 
sec 

11 


0-90 

Deg, min, 
sec 

II 


0-100 

% 

II 

- 

0-86400 

Sec 

fl 


800-1800 

Rilom. 

It 

II 

. 

s 

^ JL ^ -/a — 





• 


Description 


4 character country designation 

Current segment access count 

The month, day and year the segment 
was accessed by a vehicle. 

The revolution of the day that the 
access occurred 

The no. of the vehicle that accessed 
the segment 

The cros strack angle 
The sun ele vation angle 


The percent of cloud cover over the 
earth. 0 - none, 100 - full. 

The time during the day the access 
occurred 

The vehicle altitude above earth 
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St G hent acrts s data report 
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Figure 3. 5- 1., Samiple Segtnent Access Data jReport 
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2, A check is made to see if ICSESG matches the case no. on 
the segment ID file. Message: 

*sh*KICSESG DOES NOT MATCH THE CASE NO. ON THE 
SEGMENT ID FILE OR THIS FILE HAS NOT BEEN 
MOUNTED. 

3, A check is made to make sure that ISTIME is not less than '-• 
the earliest vehicle start time as specified on ‘the s\vath . 
table file. Message: 

^t^Si^ISTIME IS LESS THAN THE EARLIEST VEHICLE START 
l^E IVITIM OR 1V2T1M. 

4, A check is made to make sure NDAYS is between 1 and 426. 

Message: . " . 

***NDAYS IS NOT BETWEEN 1 AND 426. 

5* A check is made to make sure that the vehicle no. list in 
IVEH is between 0 and 2 and that tliere is at least one 
NONZERO entry. Message: 

««« lYEH HAS AN ENTRY NOT BETWEEN 0 AND 2 OR DOES 
NOT HAVE AT LEAST 1 NONZERO ENTRY. 

6, A check is made to make sure that a weather data file has 
been mounted. Message 

^❖❖THE WEATHER DATA FILE HAS NOT BEEN MOUNTED. 

A check is inade to make sure tliat the input quantity IVEH 
does not specify vehicles not on the swaUx files. Message. 

IVEH LIST IS NOT CGivIFATiBLE WITH NVEII AS 
•SPECIFIED ON THE SWATH FILES. 

8. A. check is made to make sure IGRDN is between 1 and 16000. 

«’!'*1GRDN DOES NOT HAVE A VALUE BETWEEN 1 AND 

16000 , 
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3, 5, 6, 3 Processing Errors 

Level 1 

20. A check is made while forming the acquisition list for a 
segment that no more than 150 acquisitions occur. 

Message: 

PROCESSING SEGMENT NNNN MORE THAN 150 

Acquisitions have occurred, no more 

ACQUISITIONS WILL BE PROCESSED. 

21« A check will be made to make sure at least 1 acquisition* 
of a segment occurs. • Message: 

***!'IN PROCESSING SEGMENT NNNN NO VALID 
ACQUISITION OCCURRED, 

22. A check is made to see if a segment's latitude is in the 
band about the equator which is not accounted fox in 
LAIT^O table. Message: 

»i'’!-.*SEGMENT NNNN LATITUDE NN NN NNIS NOT IN 

THE LATITUDE BAND SPECIFIED BY LATNO TABLE. 

23, A check is made in the weather file utility to make sure 

99 ^ Eg ^ Message. 

CUMULATIVE FREQUENCY DISTRIBUTION 
-F S^. 8 IS NOT BETWEEN 99 AND 101, 

3. 5. 6. 4 Input /Output Errors 

For sequential I/O the FORTRAN system on the UNIVAC takes control 
and a message identifying the problem and will either continue 

processing or abandon the job, If processing continues, the system counts 
the number of times this error re-occurs and if it happens a certain number 
of times, the system will abandon the job. 

For direct access l/O, the U NIVA G D, A. l/O package prints out an 
informative message, sets an error flag and allows processing to continue. 
In this program the swath table file is the only direct access file. Immedi- 
ately after the informative mes sage, the following message will print: 
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* 

40. IRRECOVERABLE I/O ERROR HAS OCCURRED 

IN READING/WRITING A RECORD ON THE 
FILE. THE JOB IS BEING ABANDONED, 

The blanks will be filled in with 'SWATH TB', or 'WEATHER' 
depending on the direct access file involved. 
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3.6 SACS 

Operational Assumptions 

• All detected errors cause the program to abort and output 
files and reports are terminated without con pletion. 

• Only one country may be selected for any one program run. •, 


• Only one case of input parameters will be processed against 
a single segment reference data file and crop window file, 

• A maximum of 426 days is allowed for the acquisition period, 

o A maximum of 4 windows is allowed. 

• A maximum of 25 acquisitions allowed per window for a train- 
ing segment and one acquisition allowed for an ordinary segment, 

3. 6. 1 Input Card Date. 

3. 6, 1. 1 Card Input Data Quantities 

Each parameter value entered via the input cards must be entered 
unless explicitly stated. Ail integers must be right justified, all hollerith 
fields must be left justified. Each parameter value is outlined as follows: 

a. XITTL - Card 01, Columns 1 through 4E. Contains the case title that 
is to be printed at the top of each output page. 

b. lACASE - Card 02, Columns 1 through 4. Contains the Integer ease 
number that identifies the AGQUISI file and each orinteci output 
page. 

c. ISGASE - Card 02* Columns 5 through 8. Contains the integer case 
number of the Segment Reference Data File (if 0 or blank ho cheek 
on input case No. will be made). 

d. IWCASE - Card 02, Columns 9 through 12. Contains the integer case 
number of the G»*op window input file (if 0 or blank no check on input 
case No, will be made). 

6; IPCe - Card 02, Golum.ns 13 through 15, Contains an integer that 
spectfies the maximum percent cloud cover to be used in selection 
of ACQUIS I file segments. This number is expressed; in tenths of a 
percent (e.g. , 15.4% v/oul d be expressed as 154). All ACQUISI File 
cloud cover percents must be less than or equal to this value. 
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f. ISANG - Card 02, Column 16 through T9. Contains an integer tnat 
specifies the minimum degrees sun angle to be used in selection 
of ACQUISI file segments. This number is expressed in hundreds 

of degrees and ranges between 0 and 90 degrees (e.g., 12,50 degrees 
would be expressed as 1250). All AGQUISi file sun angles must be 
i-han or equal to the specified angle. 

g. ICOUN - Card 02, columns 20 through 23. Contains 4 hollerith 
characters that identify the country to be selected. 

h. IRE6 - Card 02, Columns 24 through 26. Contains a three digit 
integer that identifies the region to be selected. 

i. IZONE - Card 02, Columns 27 through 29. Gontains a three digit 
integer that identifies the zone to be s elected, a zero entry 
allows all zones in a reeion to be seleeted, 

j. Output options, Columns 30 through 31. Contain a logleal flag 
that speGifies' whether the program reports are to be generated. 

A T specifies print the report; an F speeifies don't print the 
report. 

0 IPRPT-30- Daily Processing Load Report 
© ICRPT-31- Crop Window Report (sepent acquisition) 

■ 3. 6. 1. 2 Card Formats 

See Table 3. 6-1. 





T=YES, 'PRINT 
F=NO, DON'T PRINT , 
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3. 6. 1. 3 Deck Setap 

Each of the two input cards must be provided and must be supplied 
in card number order {see columns 79 and 80). 

3. 6. 1.4 Rules for Entering Data on Cards 

See Section 3. 6. 1. 2. 




m 
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3, 6, 2 Sample Car^ Inputs 


reproducibility op Tffi 

ORIGINAL pfyjp 


c*se I 


NomIMaI. CASa fOK DFUUUGIKt SACS CAS£ I 

hkaoek 




hOMTNAL 

CA4>E 

FOH 

neRuUGiNG Sag 

CASE 1 

^ ^ - 

' 

ItAiiE it 


:CW 

IPCC 

ISANG 

ItOu« 

19£G 

J2CNE IPKPf 


i ■ ■■■ 

■ 1 ■■■ ■ 

1 

90 0 

350 

USA 

1 

i> 1 T~ • 


3, 6. 3 output Report Dafea DgfMitiona 
See Table 3, 6-!'2, 

3. 6, 4 Sample Pages from Ea eh Output Report 

1, See Figure 3. 6-1 for Sample Segment A equisition. Data Repc 

2. St j Figure 3. 6-2 for Sample Daily Froeessing Load Report 
3, 6, 5 File Requirementa 


Input: GRGFW - Crop W indow File 

SEGREF - Segment Referenee File 


Output: ACQUIS - Data Acquisition File 
3, 6. 6 Error and ReGQYery 

3. 6, 6. 1 General 


All fatal errors detected by this program shall result in termination 
of program execution. The detected errors are listed in Table 3. 6-3. The 
processing logic is described as follows: 

a. Each time an error is detected, the routine A BARF is called 
with the error ID. 


b. The routine prints the error and forces the.progralii to 

terminate exeGution if the error is fatal or returns control 
if the error is non-fa tal. 


Non-fatal errors shall be printed and execution Continued. 
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Figure 3, 6>-l, Sample Segment Acquisition Data Reiport 
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Table 3, 6-3, Error Messages 



■•Non-fatai error message 





3.7 L,EM 
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The LEM executive controls the subprograms: 
Generator, YES, CAMS and CAS. 


Segment Truth 


Operational A ssumptions 

• Only 1 case may be run at a time. 

• Only 1 country may be considered in a case. 


• A maximum of 99'? Monte Carlo trials may be run in a case 
and a maximum of 100 trials may be executed on any given run. 

• A maximum of 4 crop Calendar windows and 14 additional 
prediction points can be processed. 

• A maximum of 10 regions per country can be processed. 


a A maximum of 50 zones per country can be processed. 

• A maximum of 20 strata per zone can be processed. 

e A maximum of 325 strata per country can be processed, 

• A maximum of 60 substrata per strata Can be processed. 

0 A maximum of 3200 substrata per eountry can be processed. 

e A maximum of 4000 segments per eountry can be processed. 

• AH control card input data will be echo printed. 

• All control card input data will be cheeked for errors before 
any error will cause the processing of a case to terminate. 

0 in a repetitive Monte Carlo trial case, normally the individual 
subprogram reports will be allowed to print during the first 
and last trial only. An option will exist to eliminate all reports or 
allow printing of all reports for each trial or for just the last trial. 


0 All input data files will be checked for correct ease nurnbers. 

0 The program will require iess than 20, 000 words pf storage 

in the CPU of the UNIYAG 1110, 

0 The program will have a restart capability which win allow the • 
program to continue with additional Monte .Carlo iterations start- 
ing from the last iteration of the previous run. 
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• AH files will be in country, region, zone, stratum, substratum, 
and segOaent order (to whatever level that is appropriate). 

3. 7. 1 Input Card Data. 

3^ 7. 1, 1 LEM 

3. 7'. 1. 1. 1 Input Data Description 
See Table 3, 7« 1. 

3. 7. 1. 1. 2 Card Formats 
See Table 3,7-2. 

3. 7, 1. 1, 3 Deck Setup 

Each of the four LEM control cards is required and they must be 
in card number order. In addition, control cards are always required 
for the GAMS module even if this module is not used. If CAMS is 
skipped (e, g, , ICAMS = 2), then the corresponding control cards must 
still be included. Finally, the control cards for CAS must always be 
specified. The contents and format of the CAMS and CAS control cards 
are specified in Sections 3. 7. 1. 2 and 3. 7. 1. 3, respectively, 

3, 7. 1, 1, 4 Rules for Entering Data, oh Cards 

1. Integers must be right justified. 

2. Alphanumeric fields should be left justified. 

3. Real (i, e, , floating point) fields must have the decimal point 
present. 




reproducibility of the 

GRIGINAIi PAGE IS POOll - Table 3.7-1. LEM lnpu:t Data Description 


v> ti 

Col. 

Name 

Dimension, 

Value 

Range 

1-40 

- 

TITLE 

10 

blanks 

• 


1-4 

ICASE 

i 

0 

0-9999 

6-11 

CUNTRY 

1 

-- 

— 

13-15 

NTRIAL 

1 

1 

1-999 

(NTRIAL- 
RSTART ^ 
100) 

17^19 

RSTART 

1 

0 

0-999 

20-22 

IPRINT ' 

■ 

1 

0 

0-3 

24-2.6 

■ 

STARTR 

1 

0 

0-999 

28-30 

: STARTZ 

1 

0 

0-999 

3.2-34 

ENDR 

1 • 

0 

0-999 

36-38 

ENDZ 

1 

0 

0-999 

39-41 

. ISTG 

1 

0 

0-3 

■■ 

i ■ ■ 

ICAMS 

1 

■ 

■■0 ■ 

0-3 

•t 


De sgription 




Problem header to be 
(iormat 10A6) 


printed at the top of each output page. 


Case number 

Country (must agree ■with all input files). 

Total numiber of Monte Carlo iterations at the end of the current run 
(including previous runs if this is a restart). NTRIAL must be ^ 100 
if CAS distribution file is to be generated (i. e. , if DISTFF ^ 0 in CAS 
input). 

Restart Flag; 

= n ^ 0 to restart after n Monte Carlo iterations, 0 if this is not a 
restart. 

Iteration Print Flag (for Segment Truth Generator) 

0 to print first and last iterations) 1 to print each iteration; 2 to 
print last iteration only; 3 to suppress printing. 


Starting region number. 
Starting zone number. 
Ending region number. 
Ending zone number. 


Both zero or both non-zero 


Both zero or both non-zero 


Segment Truth Flag: 

0 to vary error statistically, 1 to hold error constant using results 
from the first iteration only, 2 to hold error constant using a 
previously generated Segment Truth file, 3 to eliminate the Segment 
Truth error (error is zero). 

CAMS Error Flag: 

Usage is similar to the usage of the Segment Truth flag described 
above except that for the case in which ICAMS = 2, the CAMS Output 
file is used. 
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Taijle 3. 7- 1. LEM Input Data Deflcription (cont'd) 

ORIGINAL PAGE m POOR 




YES Error Flag; 

Usage is similar tP the usage of the Segment Truth flag described 
above except that for the case in which lYES = Z, the YES Output 
file is used. 

Segment A e^uisition Flag: 

0 to include segment acquisition conditions* 1 to eliminate segment 
acquisition conditions. 

Classification Error Flag; 

0 to vary classification error in CAMS, 1 to hold the classification 
error constant, 2 to eliminate the classification error (set it to zero). > 

* t • 

Signature Extension Error Flag; . 

0 to vary Signature Extension Error, 1 to hold Signatiure Extension 
Error eonsitant, 2 to eliminate the Signature Extension Error (error 
is zero). 

Segment Crop Calendar Error: 

0 to vary the error, 1 to hold the error constant, 2 to eliminate the 
error (error is zero). ‘ 

Flag for Group II 

"Moist Recent Non-Epoch Year" 

Historical Proportion of Wheat: 

Usage is similar to Signature Extension Error Flag described above. j 

Flag for Group III 
Multi-Y ear Proportion of Wheat: 

U sage is similar to Signature Extension Error Flag described above 

Iteration print flag for CAMS: 

U sage is similar to that of IPRINT, 

Iteration print flag for YES: 

Usage is similar to that of IPRR^T. 

Iteration print flag for CAS: 

Usage is similar to that of IPRINT. 
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Table 3. 7->l. L«EM Input Beta Description (cont'd|~ 


I Card 
* Col. 


Name 


ICSESG 


ICSECW 


I Nominali 
: Value I Banee 


Description 


.0-9999 I Case number for Segment ID file, 
0-^999 Case number for Crop Calendar file. 


ICSESH 


0-9999 


Case number for Substrata Historical file. 


ICSECE 

ICSEYM 


0-9999 Case number for CAMS Error Model file. 
0-9999 Case number for YES Error Model file. 


ICSESE 

ICSEAC 

ICSEST 

ICSECO 

IGSEYS 

RSEEDl 


RSEED2 


’ 0-9999 Case number for Signature Extension file, 

0-9999 Case number for Data Acquisition file. 

0-9999 Case number for Segment Truth file. 

0-9999 Case number for CAMS Output file, 

0- 9999 Case number for YES Output file. 

1- Initial random number, seed for Segment Truth Error (odd positive 
999999999999 integer in F -format). 


999999999 intei 
1 - I Initi 


Initial random number seed for Classification Error. 


99999999999 


RSEED3 


RSEED4 


RSEED5 


RSEED6 


RSEED7 


1- I Initial random number seed for Signature Extension Error, 

999999999999 

1- { Initial random number seed for Segment Crop Calendar Error. 

<'99999999999 

1- I Initial random number seed for Yield Error. 

999999999999 


I 1- I Initial random number seed for error in Group 11 "Most Recent Non- 

999999999999 Epoch Year" Historical Proportion of Wheat, 


I ' 1- I Initial random number seed for error in Group III "Multi-Year 

999999999999 Proportion of Wheat. " 
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Card 1 


Card 2 


Card 3 


Card 4 


Table 3.7-2. LEM Card Formats 


Header card: Alphanumierle problem header entered in Columns 1«60. T.EM is entered 
in Columns 75-77 and 01 in Columns 79-80. 
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3, 7, I, 1. 5 Re strigtions 
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1, In general, the various error sources may be independently 
varied, held constant, or eliminated by specifying appropriate 
values for the input paranieters ISTG, ICAMS, lYES, ICIoASS, 
ISEXT, ISCC, XCAS2, andICASS. However, the user may not 
specify ISTG - 0 unless ICAMS is also zero. The reason for 
this restriction is that if ICAMS = 1 or 3 for exanrpie, the CAMS 
Output file will be generated on the first iteration and then used 
on all subsequent iterations. If ISTG were zero indicating that 
the Segment Truth error was to be varied, the program would 
be in trouble beeause to vary the Segment Truth error means 
that the results on the GAMS Output file must also be variable 
even if the CAMS errors are constant or zero. The way 
around this restriction is to specify both ICAMS = Q and 

ISTG = 0 and then to specify non-zero values for ICLASS, 

ISEXT, and ISCC. In this manner the GAMS Output file will 
be written on each iteration even though the CAMS errors are 
really constant. 

2, The variable possible combinations of the input parameters 
ICAMS, ISTG, lYES, IGLASS, ISEXT, ISCC, 1GAS2, and 
ICAS3 as well as other options (CAMS eiassifieation model, 
multi- temporal sampling, and acquisition effects) are presented 
in Table 3. 7-3. 

3, On a restart run the input ease number ICASE must agree with 
the case number on both of the following files which may be input . 
to LEM: 

• CAS Cumulative Output file 

e CAS pis tribution Output file 
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Oh a restart run the CAS Cumulative Output file and the CAS 
Distribution Output file must always be specified. In addition, 
the Segment Truth file, the CAMS Output file and/or the YES 
Output file should be specified if the input flags ISTG, ICAMS, 
and/or lYES are set to 1, 2, or 3, Note, however, that if the 
GAMS Output file is specified, then it is not necessary to specify 
the Segment Truth file. 
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3. 7. 1.2 CAMS 

Some data needed by CAMS is included on the LEM control card. 

See Section 3.7, 1. 1. CAMS also requires: 

1. A control card, specifying options 

2. A multi -temporal matrix 

3. Crop calendar error coefficients 

Besides the LEM card, CAMS requires a total of 13 cards, which fanust 
be in order. See Figure 3. 7-4 for the deck setup. 

3, 7, 1. 2. 1 List of Data Quantities and Formats 

a. LEM control card, see LEM Section 3, 7. 1, 1. Data relevant to 
CAMS includes: 

ISEXT Signature extension error option, = 0, 1 simulate error 

= 2 bypass error 

ISCC Crop calendar error option, = 0, 1 simulate 

= 2 bypass 

ICLASS Classification error option, = 0, 1 sirnulate 

= 2 bypass 

lACQ Acquisition file option, = 0 include file 

= 1 no acquisition file 

iCAMS CAMS error option, = 0^ 1 simulate . 

= 2 read previous file 
- 3 bypass all errors 

SEED2 Random no. seed for classification error 

SEED3 Random no. seed for signature extensioh error 

SEED4 Random no. seed for crop calendar error 

IPRCAM Print flag for CAMS >!• this controls if a report is 
printed - the flag, ■ IREP, on the CAMS eontrol 
card, controls what is printed. 
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b, CAMS control card, see Table 3. 7>4 for the format and list of 
quantities. 

c. Multi-temporal sampling matrix. The multi-temporal 
sampling model describes the effect of the a iquisition of a 
sample segment in more than one bio^window. There are 

15 possible non-zero acquisition states for a sample segment. 
These states are shown below. 


Windows 

included 

Group 


E 

fzls' 

4 I 1 , 2 

1.3 

1,4 

2, 3|2,4 

3,4 

|l,2, 3 

h 2, 4 

1.3, 4 

2,3,4 

1.2, 3, 4 

1 

2 3 4 5 

6" 

7 

00 

10 

11 

12 

13 

14 

15 


The effect of acquisition conditions Gorresponding to any 
particular state is modeled by a weighting factor, M. For 
ease of handling, three values of M are ehoseii as program 
inputs associated with groupings of the 15 acquisition states. 
These values indicate no improvernent = 1), small knprove- 



of class if ication of the segment with a given acquisition state 
over classification in the present window (last window in the state 
definition) alone. 

The values needed are which M (Ml, M2, or M3) to use for each 
of the IS groupings, and the values of M2 and M3 (Ml = 1). Note 
that by definition, for group no, 1, 2, 3, and 4, Ml = 1, and 
Ml = 1 > M2 >M3. See Table 3, 7-S for the data and format. 

A total of eight cards are needed for the matrix, each with the 
above information, since the acquisition conditions depend On 
a, wheat type - winter or spring, b. the model - 1 or 2 (modei 1 
requires three cards, one for wheat, mixed, and other components; 
model 2 requires only one card). 

If model 1 is. being ed, the cards fOr mode.i 2 Must be pf esent 
but rnay have blank fields except for the IB (CAMS) and sequence 
number, and vice versa if model 2 is being used, since the 





T^-~^ '■:pfij.'\ ■• v-i,. 
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f 

'i J 




• ' r" ■ 

":J . ■ 

Name 

Dimension 

Nominal 

'V'alue 

Range 

Units 

Description 

i 

i IMODEL 

’ 


1-2 


Flag = 1 use model 1 , complex model 
= 2 use model 2, simple model 

IMULTI 


0 

! 

; . 

0-1 


Flag = 0 includ'e multi -temporal sampling error 
4 0 bypass multi -tern poral sampling error 

ISIGEX 


0 

0-1 


Flag 4 0 use multiplicative model of signature extension 
= 0 use additive model of signature extension 

■ 1 ISKlP 

• ! 


0 

: 0-1 


Flag = 0 skip if cannot correlate ordinary 

4 0 classify as training with training segment 

1 ITMAX 

1 



0-99 

Days 

Maxim uni no. of days between training and ordinary segment 
acquisition dates for successful correlation. 

j IREP 

i 

! 

G 

!' 

0-1 


Flag = 0 include error breakdown factors in estimate report 
4 0 print estimate report only 

1 IWIND 


4 

1-4 


From which window to take the probability of classifying as wheat 
given mixed to Galculate the proportion of pure wheat pixels; if 


in|m!fe Data - CAMS Control Card Format 


CAMERR inptit file; see file descriptions, Seetion 


c. c. 


1 

2 

3 

4 

: :5 
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M : 
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Table 3. 7 -4. Inpuit Data - CAMS Control Card Quantities 
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i ' 

m 









j 




Name 

nimension 

Nominal 
; ‘ Value 

Range 

Units 

De scription 

IGROUP 

15 

; 


1-3 


Whieh value of M to use for each acquisition state^ see 2. 1. 1, c. , 
for which windows are in each state 

= 1 use Ml (= 1), no improvement 
= 2 use M2,, small improvement 
= 3 use M3, large improvemeiit 
Resitrietion: IGROUP(l), lGROUP{2), xUROUP( 3) and IGR0UP(4) 
are always - 1, by definition* and so need not be 
inputted. 

M2 


1 ' ■ G 

i . 1 

. 0<M2<1. 1 

9 . 
1 

Value of M2, small improvement 

M3 


; , (j 

).0<M3<1. 

[ ■ 

0 

Value of M3, large improvement 
Restriction: M3 < MZ < 1, by definition* 

ISEQ 



2-9 


Sequence no^ - the matrix requires eight cards, each with the 
same format. See Deck-Setup, Figure 4, for the order and 
description, ^ 

■ ^ 
-J 


Inftiit Qata - Multi -Temporal Matrix Format 
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va.lues are not used« If data for both winter and spring is not 
available, the data may be left blank (except id and sequence 
number), but the cards must be present, and if CAMS tries to 
use the missing data, an error will be reported. If the multi^ 
temporal error bypass is specified (IMULTI* 1), then all the 
cards must be present but all the data but the id and sequence 
number may be left blank. 

d. Crop calendar coefficients. The e^eet of cr vp calendar errors 
on segment classification, particularly in an analytic sense, is 
not well established at this time. For this reason, a simple 
generic model was chosen to represent this effect. This model 
generates a bias (B) and standard deviation (b) from a quadratic 
function with user input eGeffieients. 

B = Gj (At) + G 2 (At)^ 

a = Hj (At) + (At)^ 

The value of the coefficients are to be determined from off-line 
analysis, curve fitting, etc. , to represent the observed effects. 

The values for Gl, G2, Hi, and H2 are needed. See Table 3,7-6 
for a deseriptiGn of the quantities and format. Model 1 requires 
these four values for the three eotoponeiiits., wheat, miaced, and 
other, a total of 12 values. Model 2 requires only the four values, 
since the mixed crop effect is not present. Sincf these values 
may be different for winter and spring wheat, two sets must be 
inputted. Thus, four cards are always needed: 

1, Spring wheat - model 1 - 3x4 values 

2, Spring wheat - model 2-4 values 

3, Wintex Wheat - model 1 - 3x4 values 

4, Winter wheat - model 2-4 values 


Gl 

G2 

m 

H2 

ISEQ 


Nominal 

IDim ension : ‘ ‘V'aluie 



+9.999 
i9'9. 99 
+9.999 
+99;. 99 
10-13 



*• 








ij.jiiij III! 

, Units . : ■ 


Description 




Crop calendar error coeificiient 
Crop calendar error coefficient 
Crop calendar error coefficient 
Crop calendar error coefficient 

Seqnenee nuimber - the calendar requires four cards. See Deck 
Setnp, Figure 4. . 


input Data - Crop Calendar Fornaat 
Model 1: 
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Table 3. 7-6. Input Data - Crop Calendar Coefficients Quantities 
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n model 1 is used, the cards for model 2 must be present but 
may have blank fields except for the ID, CAMS, and sequence 
number, and vice versa for model 2. If data for both winter 
and spring is not availablej the data may be left blank except for ID 
and sequence number, but the cards must be present, and if CAMS 
needs the missing data, an error will be reported. If the crop calendar 
error bypass is specified (ISCC = 2), then all the cards must still be 
present, but all the ID and sequence number fields may be left blank. 


3, 7, 1. 2, 2 Peek Setup | | 

See Figure 3.7-1, GAMS requires 13 card inputs. 

3.7. 1,2. 3 Rules for Entering Data, | I 

I V 

See Section 3, 7. 1. 1 for general rules for entering data. 

3. 7. 1.3 CAS 

3. 7. 1. 3. 1 List of Data Quantities 

See Table 3.7-7. 1| 

3.7, 1. 3,2 Card Formats 

"CAS'* is entered in Columns 75-77 of each control card, and a 
sequence number is entered in Columns 79-80, See Table 3. 7^8. | 

» ■■••■V 

3. 7. 1. 3. 3 Deck Setup | 

Each of the three GAS control cards is required (even if only seven 
or less prediction points are specified), and they must be in the proper 

order. Furthermore, the GAS control cards must follow the LEM control 

' ■ ■ ' ............ , 

cards and the CAMS control cards as specified in Section 3, 7. 1. 1. 3. 


3. 7. 1, 3, 4 Rules for Entering Data on Cards 

1 , Integers must be right justified, 

2. The prediction dates must be entered in the format 

7 (312, IX) 

with a maximum of seven dates per card. 

















NHISTY 


3-99 


TOPT 


AUNITS 


DISTFF 


IWIND 


W PRIOR 


A PREP 


IPRD 


3, 14 


>64 year 
01-12 montl 
1-31 day 


M = Number of historical years for Group III ratio 
calculations. 

H = Minimum namber of segments required for applying 
S2 regression equation. 

T - option flag: 

= 0 to set T = 0, 

= 1 to calculate T 

where T is the second terna of the variance equation for 

a * 


Units Option: 

= 1 to print area in hectares and production in metric tons, 
= 0 to print area in acres and production in bushels 

CAS distribution file flag: 

= 0 to generate CAS distribution file, 

= 1 otherwise 

Prediction bio -window flags: 

IWIND (n) =1 to process bio -window n, 

= 0 Oitherwise 

B io -v/indow pr io r itie s : 

List of bio-windows in decreasing order of priority, 
e. g. , • 4, 1, 3, 2 or , 

3. 1, 0. 0 

Print option for area and production summary report: 

= i to print report, 

= 0 otherwise 

Predietion dates (up to 14 dates): 
lPRD (l,n) - year - 1900 
IPRD (2, n) = month ^ 

IPRD (3, n) = day 
The prediction dates must be in ascending order. The first 
zero date terminates the list. 


28234-6025-RU-00 
Page 144 



■ ' P 8* gas Card Formats 

■ rr’iAI.PAGK-ISPOOR 

10 13 16 19 22 25 28 3 1 34 37 40 











' ^ ■ 

■ 

• 

' 








. 

*1*^ ■ 

■ 



to 






CO 



fS3 . 

‘ ro 



P? 

•ff! 

5 

PU 

'H ^ 



H 



■Q i 

ca 

Q 

o 

S 

o 

Q 

2 

u 

CO 

Sd : 

(1, 

^ . 

H 

CO 


. Is: 
>— 1 ; 

z 

M 

s ■ 

HH 

M 

ei 

h-J 

f4 

H4 

p. 

ffi 

: ‘2 . 


o 

: E>\. 

< 

Q 

: aq ■ 

1 

■ffl 



Pk 

> 

Pn 


1 3 5 


10 12 15 17 19 22 24 26 29 31 33 36 38 40 43 45 47 



1 

15 

: 

I 


iMigil 




Card 3 




CAS 03 


A4 12 » 

























3. 7. 2 Sample Card Inputs 


HJihDOVtR ThST CASe CCAhS t'hiRURS = 0 * VaKY STbtYtS) 

■■■ "1 G“'d’Nr'i''ir”o"T c~a“ 




LPP SlJiJLATiON 


*s POOR 

CASE 1 page 


TiTUf 

HANDOVER ft' ST Case 


(CAMS ERROR'S = ft* VAKY SlGtYES) 


UEH 01 


(Case. . CUNijf k V ’M(.' r AL HSTArT iPRiMl SIarTR' SIARTZ tN.i>R t'MDZ (Si'S ICAKS IYES IACQ ICUASS ISExT 
1 USA ft , -Cl * *■'’ ^ ® ® ® 

ISCC I(;A52 ^ICASS IPRCAI' ^iPRYlS JHJCC^S LArEL 

■ 0 • " o ' — g •' “ J i . J LEM (ik 

icsfsr. irsFCK_ icsfsh i„csEf b.„icsEYiL_ icsese icsea c i cspst ics>;co_,_icsErs 

1 1 " 1 ■'? , 2 2 J I ■ UU I 

' HSH.ni ItPf'H/? LAHEL 

, i ft fj ft V ftu O'O 0 -;io+ ft ft I Lsl'L0l*’£O 0 (MjO o f l-.E.tL.y,5 


R.SEE0J RSEtrift RSEtOS HSi EDA RSEE07 I.Atttl. 

, 1-0'P 0 0 (J 000 0 00+ 0 ftl . T 000 Ojl 0.0 0 0 00 + 001 l onOO oft 0 (10 00 + 00 1 . 1000 000 0 n 0;00+_OJ) I ^.10000 IV() 000.00+ 00 .L.Em ,0ft , 


HAHDOVK 

u It 

51 GASE 



CCAmS 

EKRypS. = V* VARY SlUtYES) 

IPP 

S1*(ULAT10n 

CASE 

l PAGE 



G A 

K S 

J 

N V U 7 

C Ij 

fJ ) 

K D L C A K D S 







•ODEL 

IKUUTI 

IS Kir. X 


15KIH 

itmax 

IREP I wind 







2 

0 


0 


1 

10 

1 

3 Cams i 







IROUPC 

b) 

C 6.) 

( H c 

»0 ( V) 

no) (in 

(1?) (13) Ubl 


M3 






2 

? 

2 

2 

a 

3 

3 

3^33 

*600 

• /1 0 0 

L A MS 





2 

a 

2 

a 

a 

a 

2 

2 a a a 

* voo 

• 70 0 

GAflS 

3 




3 

3 

3 

3 

■3 

3 

3 

3 3 3 3 

• 300 

• luo 

CAihS 

6 




2 

2 

2 

3 

3 

3 

3 

3 3 3 3 

• bOy 

• 200 

C AM S 

5 




2 

2 

2 

2 

3 

3 

3 

3 3 3 3 

• 

• ObS 

L i MS 

6 




a 

? 

3 

3 

3 

3 

3 

3 3 3 3 


*7b3 . 

. CAMS 

7 




2 

2 


t; 

3 

3 

i 

3 3 3 3 

.ibb 

• IbS 

C A MS 

S 




a" ■' 


a 


a' 

■ 3 

’i 

5 3 3 3 

*SSb 

• aas" 

CAMS ' 

9 " 


.... 

.1 

r-2 




H? 

(fi 


_Hl H2 

_ ni 

G2 


HI 

H2 


SftO, 

\ , 

fi6 ' 

' i . 0 h i) 


?,0U 

^ .'so 6 


1 * n o i .tJcTo " "2.0 0 

'.SOy”^ 

lo 

ou i 

.oao" 

2 * 0 0 

'(.a'ksiv 

bftu 

t . 

Oij 

.boo 


] .nij. 









CAMS! 1 

Sft'j 

1 « 

00 



?,na 

.sgj^ 


X. fl 0 A -. o ? 0 2.0 0 

0__ 

1. 

00 _ 1 

• ano_ 

J a .00^ 

CAMSl d 

sno" 


0 0 ^ 

1 .on a 


'a. 0 0 









CAMS 13 
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■Bk 


3.7. 2 Sample Card Inpufts (c&nt'd) 


HftMDOYtB TcST CASE 


(GAmS thRURS s 0» VARY STUtYES) LPP SIMULATION CASE 1 PAGE 3 


CAS V N P l.i I C U N I H 0 L f. A K U S 


NHISTV HH roPT AUNITS DISUK Hk!nI>C1) U) ■ (3) (^0 xPHluKO) (2) (5) (4) APREp 

fi e t I .1 1 I 1. 0 _ ,3.. I 2 -0 . __1 CAS 1 


IPIvRt 1) tpRUC p) IPiim i) TkRjC 'O IPkOI S) TPROT b) IPRDC 7) 

7b/ 1/1 s j_ , j Ay 16 r J>y'r. 'V r:^ z'^./ r, 9/ r y ..t 5 / - o /- o 


IPK%( «). ^PRI'C <1) IPKH(tO) IPKtHlU IPRIHViij IP«U(13J IPRRtlA) 

? 0 /■r’O/zA..^ . ~ r* /•*;0 ■> u/~ n /- » ^o/zS} ~i*/'-o /- v /-.(T.. cas' 3 


ft _ ^/or;i: A i AiL T’E IfeCJ-LiL Jf.I!l 


. 0 HAr/ii. ERROHS I'tTr-CTHO ttrl CUNIrOI eAMUS 




28234-6025 -RU -00 
Page 147 


28234-60^5 -RU -00 
Page 148 

3,7.3 Oatpat Report Data Pefinitions 

There are no standard output reports for the LEM control pirogram, 

3, 7. 3, 1 Segment Truth 
See Table 3,7-9. 

3. 7. 3. 2 YES 

See Table 3. 7-10. - 

3. 7. 3. 3 CAMS 

CAMS prepares two report formats; 

1, Ordinary segments data 

2, Training segments data 

Each report may be prepared with or without the individual error Gomponent 
contributions by input option (IREP), The definition of report parameters is 
given in Table 3. 7-11. The parameters in parentheses are keyed to the out- 
put format given in Figure 3. 7-2. Note that when TID = the ordinary 

segment gouM hot be properly correlated with a training segmeht. In this 
case the ordinary segment is classified as though it were a training segment. 
The resulting output has the ordinary segment titles with the training seg- • 
ment data parameters. In particular, data under Z1 is crop calendar data 
and data under Z2 is the mUlti-teropGral sampling factor and is of the form: 



2^1 

^2 

w 

(DELTA) 

(MULT) 

M 

(CROPD) 

(MULT) 

O 

0,00 

(MULT) 

Also: 

|Ti 

lAINR) = 

0. 0.0 

(TRAIND) = 
(A LOCAL) = 

3. 7. 3. 4 CAS 

0.00 

(PES) 


CAS prepares two reports: . 

1. Area and Production Summary Report 

2, Country Summary Report 

The definition of report parameters is given in Table 3. 7!- 12. 




Table 3.7-9. Segment Truth Output Report Data Definitions 


Report Name 


Syjoabgj- 


Range 


Units 


Report 


Description 


Substrata 

True PW for siub strata 


True PW for segment 


Average PW 
Error in PW 


True PM for segment 
Iteration 


PW, 


pw 


k. 

1 


PM. 


k. 

1 


1 - 9999 
0 . - 100 
1 - 9999 
0 - 100 


0 - 10:0 
4 100 


0 - 100 

1 - 100 


% 

% 


% 

% 

% 


. Tr. 

Seg. Tr. 
Seg. Tr. 
Seg. Tr. 


Seg, Tr, 
Seg, Tr. 
Seg. Tr. 
Seg. Tr. 


Substrata no. 

True proportion of wheat for a substrata 
Segment no. 

True proportion of wheat for a segment 


Average true proportion of wheat for all 
segments in a substrata 

Differenee between true PW for substrata, 
average PW for substrata 

True proportion of mixed crops for a 

S'( 


Monte Carlo iteration no. 
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OST 



Table 3.7-11 





Report Name 1 

Symbol 

Range 

Country 


- 

- 

Region 


- ■ 

1 - 999 

Z one 


- 

1 . 999 

Strata 


- 

1 - 9999 

Substrata 


-• 

1 - 9999 

Training s egment 


1 - 9999 

Cr op windo w 
window 



1-4 

Aeqi. date 
Mo/Da/Yr 


“■ ■ 

-■ 

EiSitim. Prop. 
Total error 

(PES) 

PW 

0 - 100 

(TOT) 



+ 100 

(A LOCAL) 

- 

+ 100 

Error 




Wheat 




Tot 

(ERTOT) 


+ 10 

Bias 

(ERBIAS) 


+10 

R and 

(ERRAND) 

- 

+ 10 

Mixed 




Tot 

(ERTOT) 


+ 10 

Bias 

(ERBIAS) 

- 

+10 

Rand 

(ERRAND) 

- 

+ 10 

Other 




Tot. 

(ERTOT) 

- 

+10 

Bias 

(ERBIAS) 

- 

+ 10 

Rand 

(ERRAND): 


+10 







CAMS Output Report Data Oefinitions 




Description 

4 character country ID 
Region no. within a covintry 
Zone no. within a region 
Strata no. with 
Substrata no. within a strata 
No. of training segment 
Crop window no. 

Segment acquisition date 

Estimated proportion wheat for the segment 

Differemee est. and true prop, wheat for seg. 
Local est. for ordinary segment 


0 mission error' 
for pure wheat i 
pixels ( 

Bias component I 

Random comp. 


Mod. 1 


Total, Bias, 
and random 
error com- 
ponents 


Mod, 2 


Error for mixed 
Bias component 
Random comp one 


wheat pixels 


Model 1 only- 


Commission error for non- 
wheat pixels 
Bias component 
Random component 


fJZG9?/&£Z8Z 





ul- - w;.r -■. ■: -: :.»: ; -- ■ -.r.-;-- ^ ' V ., ^ 

w ., O- - - '.; 

REPEODUdBILIlY OP THE 
OEIGINAij PAGE IS POOR 

country XXXX. JREGION XX. ZONE XXX, STRATA XXX, SUBSTRATA XXXX> TRAINING SEGMENT XXXX 
TRUE PROPORTION WHEA T = XXX. XX ^ 


A.CQ DATE ESTIM. TOTAL 

CROP WINDOW MO/DY/YR PROP. ERROR 

XXXXXXXXXXXXXXXX XX/XX/XX (PES) {TOT) 

XX/XX/XX XXX. XX XX. XXX 


ERROR 

WHEAT mixed OTHER 


C LASSIE. ERROR CROP 

WHEAT MIXED OTHER CAL 


/TOT • (ERTOT) (ERTOT) (ERTOT) (GLTOT) (CLTOT) (CLTOT) (DELTA) W 
BIAS (ERBIAS) (ERBIAS) (ERBIAS) (CLBIAS) (CLBIAS) (CLBIAS) (CROPD) M 
RAND (ERRAND) (ERRAND) (ERRAND) (CLRAND) (CLRAND);(CLRAND) O 

X.XXXX X. XXXX X.XXXX X.XXXX X.XXXX X. XXXX .XX 


MULTI- 

TEMP 

(MULT) 
(MULT) 
(MULT) 
.XX • 






If only 
estimate 
report 
* wanted, 
fromQ) 
to the right 
omitted. 


COUNTRY XXXX, REGION XX, ZON2. XXX, STRATA XXX, SUBSTRATA XXXX, ORDINARY SEGMENT XXXX 
TRUE PROPORTION \YHEAT = XXX. XX ^ 

ACQ DATE ESTIM. TOTAL 
CROP WE-fDOW MO/DY/YR PROP, ERROR 

XXXXXXXXXXXXXXXX XX/XX/XX , (PES) (TOT) 



ERROR 

WHEAT MIXED OTHER 

WHEAT 

SIG. EXT. 
mixed OTHER 


Zl 

Z2 

TRAIN 

SEG. 

TOT 

(ERTOT) (ERTOT) (ERTOT)' 

(G.ETGT) 

(CLTOT) (CLTOT) 

w 

(Z) 

(Z) 

(TID) 

BIAS 

(ERBIAS) (ERBIAS) (ERBIAS) 

(CLBIAS) 

(CLBIAS) (CLBIAS) 

M 

(Z) 

(Z) 

(TRAINA) 

RAND 

(ERHAN r5(ERRA ND)i(ERRAND)(GLRAND) (C LRA ND)(CLRAND). 

O 

(Z) 

(Z) 

(TRAIND) 


Note:: If model 2, MIXED and OTHER quantities will appear as aeros* 

If unable to correlate ordinary segment acquisition with training segment and so treated as training segment, 
TRAIN SEG, column will appear as out of range (ail ^'s} and SIG.. EXT. Z l will have the CROP CAE data, and 
SIG. EXT. Z2 the MU DTI- TEMP da ta. 


Figure 3.7-2, CAMS Estimate and Error Reports 
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Country 

Biowindow 


iteration 


Zone 


Strata 


True WA 


Es*. WA 


Mo. in subistrata group 
1 2 3 


No. of s eginenfes of gr oup 
I 2 


CV area est. pot. true 


CV error pet. true 


True yieid 


Estk yield 


Std> dev. pot. error 






•4 

- 100 

- 999 

- 999 

- 9999 

- 999999 

- 999999 

- 9999 

- 9999 

- 999. 99 

- 999. 99 

- 999. 99 

- 999. 99 


1 1000 Hectares 
or 10000 acres 


% 

% 

Quintals / 
HeGtare or 
Bushels/ 
acre 

Quintals/ 
Heetare or 
Bushels/ ■ 
acre 

Quintals / 
Hectare or 
Bushels/ 
acre 


1. , 2 . 4 character country designaiion 

1. Bio window no. 

1. , 2 . Monte Carlo iteration no. 

1. Region no. within country 

1. Z one no. within region 

1. Strata no. within zone 

:1. True wheat area 

1. , 2. Estimated wheat area 

1. ; Cihg. no. of substrata belonging to group 

1, 2 and 3 respectively 

1. The no. of segments belonging to Group 1 

and Group 2 substrata 

1. , 2. Std. deviation of estimated area error in % 
of true area 

1. , 2. Std. deviation of Monte Carlo area error in 
% of true area 

1. True yield 

1. , 2. ' Estimated yield 

1. ,-2. Standard deviation yield percent error 
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.. „ :.. M, 


Report Name 


Table 3. 7-12, GAS Oatpiit Report Data 
Symbol I Ranee I Units 


Definitions (cont’d) 




True prod. 


PRD 


0 - 999999 


Eat, prodw 


CV PRD est. pet. true 


0 . 999999 

0. 999 


CV error pet. true 


Pr edietion point 
GV anal. WA (peti 


0 - 999 


CV anal. prd. (pot. true) 


0 _ 999 


Historioal area 


area 


0 - 999999 


True area 


area 


0 - 999999 


True yield 


0 - 999 


True produGtion 


PRD 


0 - 999999 


REPRODUCIBfLnY OF 
OEi-GHHAL.' PiOF IS 


poo: 


[THE 


lOOOGO bush- 
els or 1000 
metric tons 


1 . 


100000 bush- 
els or lOO'O 
m e trio ton s 


% 


1 , , 2 . 
1 . , 2 . 


1. , 2. 


2 . 

2 . 


% 


2. 


1000 Hectares! 
or 10000 
acres 


1000 Hectares 
or 10000 
acres 


2 . 


Quintals / 
Hectare or 
Bushels/ 
acre 


2 . 


100000 bush- 
els or 1000 
metric tons 


2 . 


True production of wheat 


Estimated production of wheat 


Std. deviation of estimated production error 
in % of true production 


Std, deviation of Monte Carlo production error?; 
in % of true production :;j 


Bio window no. or prediction date 


Std. deviation of analyidc within county area 
variation in % of true area 


Std. deviation of analytic within county 
production error in % of true production 

Historical wheat area of the country 


True wheat area for the country 


True yield for the country 


True production for the country 


+0 

if ,00 
TO 

►- • 

«JI O' 

ui O I 

rvi * 


§ 

I 

O 

O 


^ - . 7^ — rr 


>r K ' ii iwtitfjtyirftTa aa 








Table 3. 7-12. CAS Oatpufc Report Data Definitions (cont'd) 


Report Name 


Symbol 

Range 

Units 

Repp3?t 


0 - 999 

% 

2. 

- 

0 - 999 

% 

2. 

- 

0 - 999 

A 

% 

2. 


0- 999 

% 

2. 

- 

0 - 999 

% 

2. 

-• 

0 - 999 

% 

2. 

* 

0 - 999 

% 

2, 


0 - 999 

% 

1 

1 

2, 

> 

\ 

: 

■w 




n 


Area Gonfidence levels 


True/WC 

Production confid. levels 
True /Err or 


True/WG 


Confidence level about die true WA using 
the Monte Carlo area error variance 

Confidence level about the estimated WA 
using the estimated area error variance 

Confidence level about the true WA using 
the estimated area error variance 

Confidence level about the true WA using 
the analytic within county area variai^ce 


Confidence level about the true production 
using the Monte Carlo production error 
variance 

Confidence level about the estimate pro- 
duction using the estimate production 
error variance 

Confidence level about the true production 
using the estimated production error 
variance 

Confidence level about the true production 
using the analytic within couuEy production 
error vzirianGe 
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3. 7, 4 Sjarnple Pages frqm Qatgu^t Reports 
3« 7, 4, 1 Segment Triith 
See Figure 3. 7-3. 


3. 7. 4. 2 YES 

See Figure 3, 7-4. 

3. 7.4. 3 CAMS 

1. See Figure 3.7-5 for Proportion Estimate Error Report - 
Ordinary Segments. 

2. See Figure 3.7-6 for Proportion Estimate Error Report - 
Training Segments. 

3. 7.4.4 CAS 

1. See Figure 3. 7-7 for Area and Production Summary Report. 

2. See Figure 3. 7-8 for Country Summary Report, 

3, 7, 5 File Requirements 

The files, input, output or intermediate, are listed by subprogram. 
Refer to Section 2. 2. 7 for a LEM overview of the file requirements. 

3. 7. 5. 1 Segment Truth 


iO File 


Input: 

SUBMST 

SEGID 

Output; 

SEGTRU 

. 2 YES 


Input; 

YESERR 

Output: 

YESOUT 


HAJiOtjivtH Ttsr easfc 


(CAHS f.KBURS = Of VA»Y ST;,*rES) LPP Sl'lULATlON 

irp”u~R^T" ' 



loLir. Pv, ^nl< 


TfrUU fi.'R 

AVitRAOE 

tPKflH 

IRiJE Pm M?h‘ 

SUBSlHAf A 
■ 13 

■ ' S^UIVSTnA-T a"” 

d\ . 

sfetKeftr 

to 0.2 

■ StKMfciJT 
23. 7200 

PW 

23. 7208 

1 N Pk 
‘2.6903 

SEGrlEMt ITtRATlCi 

10.5790 

2t’ 

20,7^06 

1 U6<. 


20.1 ‘>72 

-.6230 

9-. 099 0 


a^Vnisis 

13 i h 


as • 3976 

• V916 

12.9333 



til 7 

2A.3'760 


-IpOllS 

19.5190 

2 b 


r/?5 

1P.V86/J 

' 19 99B'6<U 

,6662 

9*Sb44a 

63 



20.V906 

20,9‘nio 

.761S 

ia*ti79 

ii 

e'O, 5^116 

1/27 

23.02/0 

ai.oava 

2.6858 

12*S99b 


16.0315 

"■")52(r 


is.Rvavj 

-2.6799 

0 * 1 Cl ? 6 



lbP9 

1/.62/0 


,/9bHi 

8. 0960 

-rs 

ap.TSVf 

r»t> 

2 u." '12 60 

22.0070 


«,7faab 

• 


t /29 

21. Ob.!? 


-•9?/6 

12*0991 



1 30.. , 

2o.2‘iai 


1 pSl 70 

lUioOti 



t/3t 

”3 1 . 0 ‘nr/’ 


-l,33w« 

liU0aV9 


/«.».0;aij_... 

t/3? 

32.0787 

32.0/67 

3.99/5 

12,34S77 



I f) .H676 

t /33 

^u.oo/i 

16.^04407 

1.13'V7 

13,2649 



t73i| 

1 1 


-2.6')26 

8. <1999 



J/35 



-1.16/7 

8.M 1 / 


a? • ivm 

i/3,. 


20.9978 

2. /l*l« 

11.2836 

^ 

71 

£0,Sf19 

lb3o 

_ 

^U. 1)^19 

19.9860 

-.aioo 

10.9106 : 



lb3] 

1 V • 9 " ? : 


-•i/4 

9.6088 




Figare 3. 7-3. Sample Segment Truth Output Report 


SSI 

G0-nH*‘SZ09-?'£Z8Z 


CAiE 


1 HAU£ i!9 


HAKOOVtH IhST eA&f 


ttAMS EWfttJPS = 0» VARY j>U.»VFS) 


LPP SiMULATlOK 


Vl,S”VlFLb FiliriMML U‘A'|a PtPljHl -tlEr<Ailt'N lu, 1 


EWlN t-RV' uii'A IvFGI t^iW 1 . /'!•>*!■ io 

PRili K. t .r>ftlK TRUf y Thl p 

' W/GV/VR r AK 

1 ^/ \/fb ft.'M 

sfTiATiJ-^ 

F b T l f' - V lELri 

PLRUFOJ 

SlANOARy UEV, , , 

Cl A=K 

'^U ^6 

l.?J 

’ in>Tf>i*./MEL 1 ar 
3 * 36 

C)« 6 V: . . . . . _ . ^ ^ — 

i/ 1/76 

6 ,‘?i 

. 6.«3 
6 ,*? 7 

'fia.Vti 

7.15 

b. 1.7 

, 6 

UVj 

* •• 6 A 

~ V/lb /76 

h,'li 

6*7 3 

e ?.91 

, OU 

GOI»MTkY UlrA UFUTOM i YOivh 30 

iu M b :0 




T Lfiir 1 f P 1 z'' 

F S 1 \ h , Y 1 > l.p 

P^RbrisiT 

SIAiNPaRU UfcV. 

iMJ/oy / iH 

~ 7 yu i IV V /Nr .r \ AK 

fJU i n./HLC l AK 


* J U i T’! • / H L C T 7 i H 

V/ 7 b 



6,73 

:S. 3 fa 

• i / t / f ^ 

^ 0 i > 

'4 • 4 f! 

S t . b V 

?-69 - 

. 1 - * ^ ' / V 

* ^/iS /76 

7 . Of. 

6 • t 6 

Tit? *01 

^•Oi: 

i/ 1/76 

. 7 , 0 o 

6 , 1 / 

1 , 6 4 

1 • 3 ti 

S/lSv' 7 t> 

7 . 0 >J 

S • 7 fa ^ 



. - — 

•l/lEZ/t. 

• 

’” 6 , 7.4 

4 . 7 fa 

• 0 u 
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Figure 3. 7-5. Sample CAMS Proportion Estimate Error Report - Ordinary Segments 
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Figure 3. 7-8. Sample GAS Country Summary Report 
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Figure 1. 7-8. Sampl-e CAS Country Summary Report (cont'd) 
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3. 7. 5. 3 CAMS 

Input: GAMERR - CAMS Error Data File 

SIGEXT - Signature Extension Data File 

CROPW - Crop Window Data File 

SEGTRU - Segment Truth Data File 

ACQUIS (Optional) - Acquisition Data File 

Output; CAMSF - CAMS Output File 

Intermediate: TAGQ - Training Segment Save Data Storage 
3, 7. 5.4 CAS 

Input: SUBHST - Substrata Historical File 

CAMSF - GAMS Output File 

YESOUT - YES output File 

Output: CASF - CAS Cumalative File 

CASDIS - CAS DistributiGn File 

Intermediate: GASDSF - CAS First Pass Work File 





28234-6OE5-RU-O0 
Page 166 


3, 7, 6 Error and Regovery 
3. 7. 6. 1 LEM 
3, 7. 6. 1. 1 gajieral 

The program will attempt to find as many errors as possible daring 
the proG‘^ssing of the inpat control cards. The program will continae 
checking for additional inpat errors if any input error is detected. There 
are two levels of error. These are; 

Level 1 - non -fatal, continae processing. 

Level 2 - job fatal. Terminate job after processing all input 
control cards. 

When a level 1 error is detected, the program will print an informative 
message and continae processing. When a level 2 error is detected, the 
program will print an informative message, set a fatal error flag, and 
continue proeessing. When ail control cards have been processed the 
program will continae executing if no fatal errors were found or will 
return control back to the operating system if at least one fatal error is 
detected. 

The errors which may be detected by the LEM control program 
itself are described below. Any error eonditions which are detected by 
CAMS, YES, or CAS will be described separately in Sections 3. 7, 6, 3, 

3. 7. 6. 2, and 3. 7. 6. 4, respectively. 

3.7. 6. 1,2 Input Errors Detected by LEM 

1. Message; 

TOO MANY MONTE CARLO TRIALS REQUESTED, NTRlAL=n, 
RSTART - r. MAX. NO. OF TRIALS PER RUN IS m. 

Meaning; 

On the LEM control cards the user has specified n - r Monte Carlo 
trials for the current run but the program permits a maximum of 
m trials for any single run. 

Remedy ; 

Fatal error -- the user should cheek NTRIAL and RSTART and 
be sure NTRiAL-RSTART does not exceed the maximUra allowable 
value. 
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2. Message: 

RSTART = r MUST BE LESS THAN NTRIAL = n. 

Meaning: 

RSTART, Che final iteration number from the previous run 
from which the user is trying to restart, must be less than 
NTRIAL, the total number of iterations desired at the end of 
the current run. 


Fatal error -- RSTART is fixed. Hence, NTRIAL must be 
increased on the LEM control cards. 

3. Me ssage: 

STARTR = nj MUST BE BETWEEN 0 AND ENDR = n^. ENDR 
MUST BE . LE. m. 

Meaning: 

The starting region and the ending region must satisfy 
the inequalities 

0 ^ n^ n 2 ^ rn 

where m is the maximum region number* 



Fatal error -- the user should cheek STARTR and ENDR on the 
LEM control cards to be sure they satisfy the above inequalities. 

4. Message; 

STARTZ = nj MUST BE BETWEEN 0 AND ENDZ = n2* ENDZ 
MUST BE . LE. m. 

Meaning: 

The starting zone and the ending zone n 2 naust satisfy the 
inequalities 

0 ^ n j ^ n2 ^m 

where m is the maximum zone number. 
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Fatal error -- the user should check STARTZ and ENDZ on 
the LEM control cards to be sure they satisfy the above 
inequalities. 

Message; 

ISTG = n^. ICAMS - AND lYES = MUST ALL BE 0, ifz, 
OR 3. 

Meaning: 

One or more of the parameters ISTG, ICAMS, and lYES have 
an illegal value specified. The only allowable values are 0. 1. 
2, or 3. 


Fatal error -- specify the proper value (s) for the offending 
parameter(s) on the LEM control cards. 


Message: 


IF ICAMS IS NONZERO, THEN ISTG MUST BE NONZERO. 
I.E. IF the cams errors are HELD CONSTANT, THEN 
SO MUST THE SEGMENT TRUTH ERROR. 


] 


Meaning; 

Self-explanatory. The user cannot vary the Segrient Truth error 
while holding the CAMS errors constant by setting ICAMS 4 0. 
However, it is possible to hold the Segment Truth error constant 
while varying the CAMS errors. 



1 


Fatal error 
cards. 


change either ICAMS or ISTG on the LEM control 


4 
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Message; 

CASE NUMBER = OR COUNTRY C^ FROM SEGMENT ID 
FILE DOES NOT AGREE WITH INPUTS ICSESG = n2 AND 
CUNTRY = C^. 

Meaning; 

Possibly the wrong Segment ID file has been specified, or the" 
wrong values have been specified for the parameters ICSESG 
and CUNTRY on the LEM control cards. 

Remedy : 

Fatal error -- mount the proper Segment ID file Or specify the 
correct values for ICSESG and CUNTRY on the LEM control 
cards. It might be necessary to dump the header record of the 
Segment ID file. 

Me s sa ge : 

CASE NUMBER = n^ OR COUNTRY Cj FROM CROP WINDOW 
FILE DOES NOT AGREE WITH INPUTS ICSECW = n^ AND 
CUNTRY = G^. 

Meaning; 

Possibly the wrong Crop Window file has been specified, or the 
wrong values have been specified for the parameters ICSECW 
and CUNTRY on the LEM control cards. 

Remedy ; 

Fatal error -- mount the proper Crop Window file or specify the 
correct values for ICSECW and CUNTR-Y on the LEM control cards. 
It might be necessary to dump the header record of the Crop 
Window file. 
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9. Message; 

CASE NUMBER = OR COUNTRY Cj FROM CAMS ERROR 
FILE DOES NOT AGREE WITH INPU TSICSECE =: AND 

GUNTRY = C2- 

Meaning; 

Possibly Qie wrong CAMS erro? file has been specified, or fke 
wrong values have been specified for the parameters IGSECE 
and CUNTRY on the LEM control cards. 

Remedy ; 

Fatal error -- mount the proper GAMS error file or specify 
the correct values foriCSEGE and GUNTRY on the LEM control 
cards. It might be necessary to dump the header record of the 
CAMS error file, 

10. Message; 

CASE NUMBER = n^ OR COUNTRY C^ FROM THE SICaSFATURE 
EXTENSION FILE DOES NOT AGREE WITH INPU TS ICSESE = 
AND CUNTRY = C^. 

Meaning; 

Possibly the wrong signature extensiGn file has been specified, 
or the wrong values have been speeified for the parameters 
ICSESE and GUNTRY on Uie LEM control cards. 

Remedy ; 

Fatal error -- mount the proper signature extension file or 
Specify the eorreet values for IGSESE and CUNTRY, It might 
be necessary to dump the header record of the signature 
extension file. 


a^id iMCW i sjTi i L 

■ ,■ -'X ■ 
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11. Message; 

CASE NUMBER = Hj OR COUNTRY Cj FROM THE DATA 
ACQUISITION FILE DOES NOT AGREE WITH INPUTS ICSEAC = 
and CUNTRY = 

Meaning; 

Possibly the wrong Data Acquisition file has been specified, or 
the wrong values have been specified for the parameters ICSEAC 
and CUNTRY on the LEM control cards. 

Remedy ; 

Fatal error -- mount the proper Data Aequisitipn file or specify 
the correct values of ICSEAC and CUNTRY. It might be necessary 
to dunap the header record of the Data Acquisition file. 

12. Message; 

CASE NUMBER = OR COUNTRY Cj FROM THE YES ERROR 
MODEL FILE DOES NOT AGREE WITH INPUTS ICSEYM= 

AND CUNTRY =: C^. 

Meaning: 

Possibly the wrong YES Error Model file has been specified or 
the wrong values have been specified lor the parameters IC SEYM 
and CUNTRY on the LEM control cards. 



Fatal error -- mount the proper YES Error Model file or specify 
the correct values for IGSEYMand CUNTRY. It might be necessary 
to dump the header record of the YES Error Model file. 

13. Me ssage; 

CASE NUMBER = n^ OR COUNTRY Cj FROM THE SUBSTRATA 
HISTORICAL FILE DOES NOT AGREE WITH INPU TS ICSESH = 

AND CUNTRY = C^. 
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Meaning; 

Possibly the wrong Substrata Historical file has been specified 
or the wrong values have been specified for the parameters 
ICSESH and CUNTRY on the LEM control cai ds. 


Remedy ; 

Fatal error -- mount the proper Substrata Historical file or 
specify the correct values for ICSESH and CUNTRY. It might 
be necessary to dump the header record of the Substrata 
Historical file. 

14. Message : 

CASE NUMBER = OR COUNTRY = C^ FROM THE SEGMENT 
TRUTH FILE DOES NOT AGREE WITH INPUTS ICSEST = n^ AND 
CUNTRY = C2- 

Meaning: 

Possibly the 'vrong file has been specified as the Segment Truth 
file, or the wrong values have been specified for ICSEST and 
CUNTRY on the LEM control cards. 

Remedy : 

Fatal error -- mount the proper Segment Truth file or specify 
the correct values for ICSEST and CUNTRY. It might be necessary 
to dump the header record of the Segment Truth f.le. 

15, Message ; 

CASE NUMBER - n^ OR COUIITRY = Cj FROM THE CAMS 
OU TPU T FILE DOES NOT A GREE WITH INPU TS ICSECO = n2 
AND CUNTRY = C 2 . 

Meaning; 

Possibly the wrong file has been specified as the CAMS Output 
file, or the wrong values have been specified for ICSECO and 
CUNTRY on the LEM control cards. 


Remedy ; 
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Fatal error -- mount the proper CAMS Output file or specify 
the correct values for ICSECO and CUNTRY, It might be 
necessary to dump the header record of the CAMS Output file. 


16, Message: 

CASE NUMBER = n^ OR COUNTRY = C^ FROM THE YES 
OUTPUT FILE DOES NOT AGREE WITH INPUTS ICSEYS = 

AND CUNTRY = C^. 

Meaning; 

Possibly the wrong file has been specified as the YES Output 
file, or the wrong values have been specified for ICSEYS and 
CUNTRY on the LEM control cards. 

Remedy; 

Fatal error -- mount the proper YES Output file or specify the 
correct values for ICSEYS and CUNTRY. it might be necessary 
to dump the header record of the YES Output file. 

17. Message: 

CASE NUMBER = OR COUNTRY = Cj FROM THE GAS CUM 
OUTPUT FILE DOES NOT AGREE WITH INPUTS IGASE = 

AND CUNTRY = C^. 

Meaning; 


On a restart run the ease number and country of the CAS Cum 
Output file must agree with the parameters iCASE and CUNTRY 
on the LEM control cards. Possibly the wrong file has been 
mounted, or the wrong values have been specified for IGASE and 
CUNTRY on the LEM control cards. 

Remedy ; 

Fatal error -- mount the proper CAMS Cum file or specify the 
proper values for ICASE and CUNTRY. It might be necessary 
to dump the header record of the GAMS Cum file. 


18. 
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Message: 

CASE NUMBER - OR COUNTRY = FROM THE CAS 
DIST Or^PUT FILE DOES NOT AGREE WITH INPUTS 
ICASE = AND CUNTRY = 

Meaning: 


On a restart run the case number and country of the CAS 
Distribution Output file must agree with the parameters ICASE 
and CUNTRY on the LEM control cards. Possibly the wrong 
file has been mounted, or the wrong values have been specified 
for ICASE and CUNTRY on the LEM control cards. 

Remedy : 

Fatal error -- mount the proper CAS Distribution file or specify 
the proper values for ICASE and CUNTRY. It might be necessary 
to dump the header record of the CAS Distribution file. 

19. Message! 

IMPROPER HEADER LABEL ON FILE filename. 

LABEL = label. 

Meaning: 

The file label specified in the first two words of the header 
record of the file "filename" does not correspond to the expected 
label. Possibly die wrong file has been mounted. 

Remedy : 

Mount the proper file. It might be necessary to dump the header 
record of the designated file. 

20, Message: 

IMPROPER LABEL and SEQUENCE NUMBER ON A LEM 
CONTROL CARD. LABEL AND SEQ. NO. = . 

Meaning; 

The four LEM control cards are supposed to have LEM Oi 
entered in Col'^J^^s 75-80 (where i = 1,2, 3, or 4). Possibly • 
the control cards are out of order. 
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Be sure the LEM control cards are in the proper order and 
that the label and sequence numbers are entered properly. 


21, Me s sage; 

ITERATION NUMBER NT FROM filename FILE = n DOES 
NOT AGREE WITH RSTART = r FROM LEM CONTROL CARD, 
where filename = CASGDM or CASDIS 

Meaning: 

The iteration number NT from both the GAS Cumulative file 
(CASCUM) and the CAS Distribution file (CASDIS) must be con- 
sistent with the value of RSTART specified on the LEM control 
cards. Possibly the wrong file has been specified or else 
RSTART is specified ineorreGtly. 

Remedy; 

Fatal error -- mount the proper file or specify the correct value 
for RSTART. 
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3.7.6. 1.3 Processing Erroys 

£ach of the modules STG, CAMS, YES, and GAS performs its own 
error checking during execution. The error messages for CAMS and 
CAS are described separately in Sections 3. 7. 6. 3 and 3.7, 6,4. The error 
messages for the Segment Truth Generator and the YES module are described 
below. 


, 6. 1. 3. 1 Processing 

Errors 

Detected by the Segment Truth Generator 

Message: 



THE SEGMENT ID FILE AND THE SUBSTRATA HISTORICAL FILE 

ARE INCONSISTENT 

« 


SEGID 

SUBHST 

REGION 

rl 

r2 

ZONE 

zl 

z2 

STRATA 

si 

s2 

SUBSTRATA 

kl 

k2 


Meaning: 

The region, zone, strata, and/or substrata ID*s do not agree between 
the SEGID file and the SU BUST file. Perhaps the wrong file has been 
specified for one or both files. This error message was intended 
primarily for checkout purposes. The error should not occur during 
production usage. 

Remedy : 

Non-fatal error -- the segment from the SEGID file will be dropped 
and execution will continue. However, the user should cheek both files 
to be sure the proper files have been specified, it may be necessary 
to dump part or ail of one or both files. 

2. Message : 

SEGMENT IS NOT IN IDSEG FROM SUBHST FOR REGION r, ZONE z, 
STRATA s, SUBSTRATA k 
SEGMENT WILL BE DROPPED. 
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Meaning; 

The indicated segment ID cr from the Segment ID file was not found in 
the array IDSDG read from the SUBHST file. Apparently the SBGID 
file and the SUBHST file are inconsistent. 



Non-fatal error — the indicated segment will be ignored and execution 
will proGeed, However, the user should cheek both files to be sure the 
proper files have been specified. 


Message: 

ERROR RETURN FROM BETAD ROUTINE. 
lER = f . PWK = PW, SIGMA = CT PWKI = PW. 


Meaning: 

An error return from the Beta Distribution routine has occurred in 
STG. The error flag f indicates the nature of the error. 


f = I 
£ = 2 


XBAR - PWK (X = PWj^) is nOt in i)ie range 0 ^ X ^ 1. 
X was reset within BETAD, 


O’ not in the range 0 ^ <y ^ X 


” 5T-J 


X 


X+ C 


where C = 10 

or was reset Within BETAD, 


f - 3 the random number could not be found within 35 

iterations via the inverse incomplete Beta function 
method, so XI set to XBAR. 

Remedy : 

The error is non-fatal and the program continues execution. 


Message: 
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ERROR RETURN FROM BETAD ROUTINE. 
lER = f . PMMEAN = PMj^ SIGMA = f 

MeaniTig: 

Similar to error message 3 described above except that here 
XBAR = PMMEAN 
(X = PMj^> 


Remedy : 

If lER = 1 Or 2, the error is non-fatal atid the program Gontinues 
execution. 

V 

If lER = 3, then the user should eheek the values DELTPM, and 

CVj from the SUBHST file. 

Me § sa ge ; 

WARNING... NSEGS = n . NE. NSEG 

(FROM SUBHST) = m 

SEGMENT IDS MAY BE INCORRECT. 

Meaning; 

The number of segments read from the SEGID file does not agree with 
NSEG, the number of segments on the SUBHST file. Po ssibly the two 
files are inconsistent. The error is considered to be nc n-^f atai by the 
program, but the error may be quite serious and the results should be 
regarded with, suspieion. 

Remedy: 

Cheek to be sure the SEGiD and SUBHST files are consistent^ It might 
be necessary to dump all or part of both files. 

Me s sage : 

WARNING. . . NO SEGMENTS PROCESSED BY SEGMENT TRUTH 
GENERATOR. 


Meaning; 
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For some reason the Segment Truth Generator failed to process any 
segments. This is a very degenerate situation and should not normally 
be encountered in production runs. The most likely explanation is that 
all of the substrata processed by ST G had zero segments, 

Remer 'y: 

Non-fatal error -- however, this condition will probably cause an abort 
in CAMS, Cheek the SU BUST file. Also check STARTR, STARTZ, ENDR, 
ENDZ, 

3.7. 6, 2 YES 

3. 7.6. 2.1 General ” 

There are two possible errors besides system error s in YES, which 
generate the following messages: 

YES INPUT FILE (YESERR) - BEGINNING REGION AND ZONE 

NOT FOUND 

YES INPUT FILE (YESERR) - ENDING REGION AND ZONE 

NOT FOUND 

if the beginning region and zone are not found, this i. a fatal error, and 
causes return of control immediately to LEM. If the ending region and 
zone are not found, (his generates a warning, but LEM will continue, 
having processed all records from the beginning region and zone to the 
end of file. The beginning and ending regions and zones are speciiied on 
the LEM control cards. 




i.:s vr 


3. 7. 6. 3 CAMS 
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3. 7. 6. 3, 1 General 

See Section 3. 7. 6. 1, 1 for a surrimary of ov^erall error handling. 

3. 7. 6. 3, 2 In:pyit Data Errors 

1. CAMS XX MODEL NOT 1 OR 2 - X 

Fatal error - model number on CAMS control card, CAMS 01, 
is not 1 or 2. See Table 3. 7-4. 

2. CAMS XX ITMAX NOT BETWEEN 0 AND 99 - XX 

Fatal error - ITMAX value on CAMS control card bad. See 
Table 3. 7-4. 

3. CAMS XX IWIND NOT BETWEEN 0 AND 4 - X 

Fatal error - IWIND value bad {note that 1-4 good values, 0 = 
default value, set to 4) on CAMS control card. See Table 3,7-4. 

4. BAD CAMS ID OR SEQUENCE NO. -XXX XX 

Fatal error - CAMS control cards bad, perhaps out of order, 
or missing one. See Figure 3, 7-1. 

5. CAMS XX CROP CALENDAR COEF. OUT OF RANGE - XXXXX 

Fatal error - crop calendar coefficient should bd between 
+9.999 or +99. 99. See Table 3. 7-6. 

6. CAMS XX BAD MULTI-tEMPORAL MATRIX VALUE M(X) ^ XXXX 

Fatal error - M2 not in range M3 ^ M2 < 1.0 or M3 not in range 
0S M3^ M2. See Table 3, 7-5. 

7. CAMS XX BAD MULTI- TEMPORAL MATRIX VALUE IGROUP (XX) - X 
Fatal ei*ror - IGROUP value not 1, .2, or 3. See Table 3.7-5,. 




t.' ■ . 
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3. 7. 6,3,3 Processing Errors 

1. CAMS INPUT TAPE SEGTRU - BEGINNING REGION AND ZONE 
NOT FOUND 

Fatal error - start region and zone specified on LEM control 
card, LEM 02, is not present in file. 

2 . CAMS INPUT TAPE SEGTRU - ENDING REGION AND ZONE.. 

NOT FOUND 

Warning - end region and zone were not found, so CAMS processed 
all records until EOF (end-of-file). 


3. CAMS INPUT TAPE XXXX - MISSING RECORD 

Fatal error - input tape does not eorr elate correctly with key 
tape SEGTRU - perhaps wrong file mounted for SEGID or this 



4. BETA DISTRIBUTION ERROR - FLAG = X 

Warning - Beta distribution subroutine, BETAD, returns error 
for mean productioh error: 

a, FLAG = 1 mean not in range 0 ^ mean ^ 1 so if mean > 1, 
mean set to 1; mean < 0, mean set to 0. 


b. FlaAG = 2 sigma not in range 0 ^ sigma ^XBAR 

so was reset within BETAD, 

c, FLAG = 3 the random number could not be found within 3S 

iterations via the inverse incomplete Beta function 
method, so XI set to XBAR. 


1 -XBAR 
XBAR+10 


-4 
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3. 7. 6. 4 CAS 

3. 7. 6, 4, 1 General 


The program will attempt to find as many errors as possible during 
the processing of the input control cards. The program will continue 
checking for additional input errors if any input error is detected. There 
are two levels of error. These are: 

Level 1 - noh>>fatal, continue processing. ' 

Level 2 - job fatal. Terminate job after processing all input 
control cards. 

When a level 1 error is detected, the program will print an informative 
message and continue processing. When a level 2 error is detected, the 
program will print an informative message, set a fatal error flag, and 
continue proGessing. When all control cards have been processed, the 
program will continue exeGUting if no fatal errors were found or will 
return control back to the operating system if at least one fatal error is 
detected. 

The errors which may be deteeted by the CAS Simulator are 
described below. 

3. 7, 6, 4. 2 Input Error s Petected by CAS 
1 , Me ssage; 

IMPROPER LABEL AND SEQUENCE NUMBER ON A CAS CONTROL 
CARD. LABEL AND SEQ. NO. = 

Meaning! 

Fatal error -- the three CAS control cards are supposed to have 
CAS Oi entered in Columns 75-80 (where i = 1, 2, or 3). Possibly 
the control cards are out of order. The CAS control cards must 
always be preceded by the LEM control cards and the CAMS control 
Gards. 




Retnedy: 


( 
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i. V 


Be sure that the LBM, CAMS, and CAS control cards are all present, 
and are in the proper order. Chech the label and sequence numbers 
to be sure they are entered properly. 

2, Message; 

NHIST = m IS OUT OF RANGE. 

(1 .L.E. NHIST ,LE. 20 ) 

Meaning; 

Fatal error — NHIST must satisfy 1 ^ NHIST ^ 20;. 

Rernedy: 

Change the input value of NHIST or change the limits of NHIST within 
the CAS Simulator (subroutine CASIN). 

3, Message : 

HH = b IS OUT OF RANGE 
( 3 . LE. HH . LE. 99 ) 

Meaning: 

Non-fa tal error -- HH must be within range 
3 £ HH s 99 

2 

in,order to apply the regression relation for S . 

Remedy: 

Change the input value of HH or change the limits for HH within the CAS 

Simulator (subroutine CASiN). The program will set HH = 99999 so 

2 

that the first formula for S^ will always be used. 

4, Me s sage: 

ILLEGAL WINDOW SPECIFIED IN W PRIOR = W,, W_, W-, W,. 

12 3' 4 

(EACH WINDOW MUST BE 1-4 OR 0) 


. /V -r*'-: - * 
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Meaning; 

Fatal error — an imiproper value has been specified for one or more 
of the windows in the array W PRIOR, The only allowable values are 
0, 1, 2, 3, or 4, 

Renvedy; 

Correct the offending values. 

5, Message: 

ALL ENTRIES IN WPRiOR ARE ZERO 
Meaning; 

Fatal error - - each entry in the array WPRIOR is zero, but at least 
one window 1-4 must be specified. . 

Remedy; 

Specify at least one non-zero window number in the array WPRIOR. 

6, Message: 

ILLEGAL PREDICTION DATE yy/mm/dd SPECIFIED. YEAR = yy 
MUST BE , GE, 64, MONTH = mm MUST BE 1-12, DAY MUST BE 

1-31 

Me^'ning; 

Fatal error -- an illegal predietion date has been specified in the 
array IPRD, The prediGtion date must satisfy 
year ^ 64 
1 ^ month ^ 12 
1 ^ day ^31 

Note: Dates such as Feb. 30 or Sept, 31 will be accepted by the 

program without being recognized as being in error. 


CprreGt the offending dates. 




:t: 
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7. Message; 

PREDICTION DATES NOT IN ASCENDING ORDER OR DUPLICATES. 


Meaning; 

Fatal error — the prediction dates entered in the array IPRD must be 
in ascending order with no duplicates. 



Enter the prediction dates in ascending order and eliminate any 
duplicates. 

3, 7, 6.4. 3 Processing Errors Detected by CAS 
1. Message; 

DIVISION BY ZERO NOT ALLOWED 
EQN. (n), symbol = 0, 

Meaning: 

Fatal error -- the program detected a zero divisor in attempting to 
compute equation (n). The offending zero divisor is indicated 
symbolically by "symbol. " 


Remedy; 


The user should attempt to discover why the indicated quantity was 
zero. Usually potential zero divisors were supposed to be anticipated 
during the analysis leading to the coding of the CAS Simulator. The 
program logic should avoid the calculation of zero divisors. 

2. Message: 

IF NT = 1. VARIANCE ERRORS AND CONFIDENCE LEVELS CAN 
NOT BE COMPUTED AND WILL ARBITRARILY SET TO ZERO. 

Meaning: 


Non-fa tal error -- on the first Monte Carlo iteration it is not possible 
to compute the variance errors VEA^, VEP^, and VEY^ and the 
confidence levels CLWA, CLPRD,' etc. These values will arbitrarily 


set to zero. 
Remedy: 


Not required. 
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3, Message: 

STARTING REGION r ^ NOT FOUND ON filename FILE, 
(where filename i3 YESOUT, SUBHST, or CAMSF) 

Meaning; 


Fatc.l error -- the starting region r^^. 


specified by STARTR on the 


CAS control cards was not found on the indicated file. Either 


STARTR is incorrect or something is wrong with the indicated file. 


Remedy: 

Change STARTR or mount the correct file. It might be necessary 
to dump part of the file to determine the starting region and zone 
numbers. 


4. Message: 

STARTING ZONE z , NOT FOUND ON filename FILE 

St 

(where filename is YESOUT, SUBHST, or CAMSF) 

Meaning; 

Fatal error -- the starting zone z specified by STARTZ on the 

S V 

CAS control cards was not found on the indicated file. Either 
STARTZ is incorrect or something is wrong with the indicated file. 


Remedy; 

Change STARTZ or mount the correct file. It might be necessary to 
dump part of the file to determine the starting region and zone numbers. 

5* Message: 

ENDING REGION r^^^ NOT FOUND ON filename FILE 
(where filename is YESOUT, SUBHST, or CAMSF) 

Meaning; . 

Non^fatal error -- the ending region specified by ENDR on the 

CAS control cards was not found on the indicated file. Either ENDR 
is incorrect or something is wrong with the indicated file. The pro- 
gram will use ail regions up to the end of data on the file. 
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Rennedy: 


Change £NDR. Zero is a permissible vaiae indicating to use aU 
regions up to the end of data. 

Message: 

ENDING ZONE z , NOT FOUND ON filename FILE 
end 

(where filename is YESOUT, SUBHST, or CAMSF) 


Meaning: 

Non-fatal error -- the ending zone specified by ENDZ on the 

GAS control cards was not found on the indicated file. Either ENDZ 
is incorrect or something is wrong with the indicated file. The pro- 
gram will use all zones up to the end of the last region or the region 
indicated by ENDR. 

Remedy; 

Change ENDZ. Zero is a permissible value indicating to use all 
zones of the final region (ENDR), 


Message: 

ZERO PREDICTION DATES ON YESOUT FILE FOR REGION r, 

ZONE z, STRATUM s (DATA RECORD n) 

Meaning; 

Fatal error -- all six prediction dates from the YESOUT file are 
zero for the indicated region, zone, and stratum. Thus the program 
cannot determine which value of estimated yield to use. Something 
must be wrong with the YESOUT file. 

Remedy; 

Dump out part of the YESOUT file to check the prediction dates and 
yields. In particular record n+1 should be checked. 

Message: 

ILLEGAL GROUP NUMBER g FROM SUBHST FOR REGION r, ZONE z, 
STRATUM s, SUBSTRATUM k (DATA RECORD n) 
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Meaning; 

Fatal error - - the group number g from the data record of the 
Substrata Historical File SUBHST is not 1, 2, or 3, indicating 
something is wrong with the SUBHST file. 

Rernedy; ’ 

Dump out the n+1 record of file SUBHST to check the group number 
and other substrata data. 

Message: 

NAGR = n OR NA = m FROM FILE SUBHST ARE ZERO. GROUP 
NUMBER g IS CHANGED TO 3. 

Meaning: 

Non-fatal error -- one or both of the quantities NAGR and NA from 
the Substrata Historical File are zero. Thus the group number g 
was changed to 3 by the CAS Simulator. 


The input to the LUMP program, which generated SUBHST. If 
NAGR or NA are zero, then the group number should be 3. 

Message; 

iNCONSISTENCY BETWEEN YESOUT AND SUBHST FILES. 



RECORD 

REGION 

ZONE 

STRATA 

YESOUT 

"l 

^1 

®1 

"1 

SUBHST 

”2 

^2 

^2 

®2 

Meaning : 





Fatal error 

-- the region, 

zone, and 

strata from the YESOUT 


SUBHST files do not agree. Agreefnent was supposed to be assured 
by the YES module. 
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f 

V > 



^ ) 


Remedy; 

This error should never occur in production. The logic of the YES 
and CAS modules should be carefully checked. Also it might be 
necessary to dump portions of the two files. 

11. Me St sage? 

• INCONSISTENCY BETWEEN SUBHST AND CAMSF. 


SUBHST 

RECORD REGION 

“1 'I 

ZONE 

^1 

STRATA 

SUBSTRATA 

‘'l 

CAMSF 

“2 

'2 

®2 

®2 

^2 

Meaning; 
Fatal error 

the region, 

zone. 

strata, and substrata from, the 


SUBHST and CAMSF files do not agree. Agreement shpuld have been 
assured by the CAMS Enodule, Actually, the CAMS module uses the 
CROPW file rather than SU BHST,; but the two files should agree with 
each other and thus with CAMSF. 

Remedy; 

This error should never oeeur in production. The logie of the CAMS 
and GAS modules should be carefully cheeked. Also, it might be 
necessary to dump portions of the two files. 

12, Message; 

ERROR RETURN FROM BETA DISTRIBU f JON SUBROUTINE. 

ERROR FLA G = n. 

Meamng; 

Non-fatal -- an error was detected by the BE TAD subroutine while CAS 
was attempting to compute the production wheat for the most 

recent non-epoch year. 

The meaning of the error flag is as follows: 
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1. 

2 . 


X not within range 0 ^ X ^ 1 so was reset within BETAD. 
or not within range 

0 « a 5 X 


vs 


€ 

where C - 10 

-^o O was reset within BETAD. 


the randpns number cpuld not be found within 35 
iterations via the inverse incornpiete Beta function 
method, so XI set to XBAR. 


X is the mean value PW 


K 


a is the standard deviation 


a = CVj =!* PWjj 


Message.; 

NO SEGMENTS IN SUBSTRATA k, STRATA s, ZONE z, REGION r 
(SU BUST RECORD n). 

Meaning; 


Fatal error -- the program detected a group I substrata with no 
segments. Only group II or group IH substrata with no segments 
are 


Cheek record n+1 of the SUBHST file (in partieular check GRPNO and 
NSEG). This error should never occur during production. 

Message: 

ZERO OR NEGA TIVE DIVISOR IN COMPU TING TAU2A, SIGM2S 
tEQS. 93D-93F) 

Meaning: 


1/2 


Fatal error -- the denominator DENOM = IIWA12 + RN2(v) * (MY V 12) 

2 2 

in Eqns. 93d, 93e, and 93f in the Galculation of Tg andCTg is zero or 
negative (subroutine DSIO). This probably indicates that the group 1, II 


historical wheat area and the multiyear variance are zero. The pro- 
gram logic should never reach this point (see message IS), 
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Retnedyt 

Modify tile input to the LUMP program so the historical wheat area 
is nOn-zero. 

Messaget 

WARNING. .. HIST PW - PWFOR SUBSTRATA k, STRATA s, 
ZONE z, REGION r. 

GROUP NO. CHANGED TO 3. 

Meahihg: 


Non-fatal -- the program, will not accept a group I or group n sub- 
strata with a zero (or negative^ value of historical PW from the 
SUBHST file. 

Remedy: 

None required -- the program will automatically change the group 
number to 3 and proceed. The user may wish to enter a non-zero 
value of HIST PW in the LUMP input data. 

Message: 

INPUT PREDICTION DATE (m) = d 

. LT. ALL PREDICTION DATES ON YESOU T FILE FOR STRATA s. 
ZONE z, REGION r (RECORD n). 

Meaning: 

tb. til 

N on-fatal -- the m Zulu prediction date (obtained from the m 

prediction date on the CAS input control card data) is less than all 

prediction dates on the YESOU T file for the indicated stratum on the 

n*^*^ YESOUT data record). 

Remedy ; 

The error is non -fatal. The program will drop the indicated stratum 
and proceed. However, the user may wish to check the prediction 
dates entered on the CAS control cards. 
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Message; 

TOO MANY MONTE CARLO ITERATIONS FOR THE CAS DISTRIBUTION 
FILE. 

Meaning.: 

Fatal error --a maxiraum of 100 Monte Carlo iterations is allowed^ 
if the CAS distribation file is to be generated. 

RenftedY; 

Specify NTRIAL ^ 100 in the LFM control card data or specify DISTFF = 0 
in the CAS control card data. If more than 100 iterations are required 
and if the CAS distribution file is desired, then the dimensions of the 
arrays CASDSR and BUFFR in common block /GASGNM/ naay have to be 
increased. Also the routine RWDISF would have to be modified* 

Message; 

Symbol IN EQ. n = a 
REF. value = b 


Meaning; 

Non-^fatal -- in subroutine YSUBj which calculates a quantity Y, the 
argument a for the square root is negative, which could cause trouble. 
To avoid the problem, the program resets 
a = 0 . if a < 0 

and prints a warning if 

jaj s b X 10"*^ 


Renaedy! 

None required the error is non-fatal, and execution will continue with 
* 30 * 

a = 0 and Y = 10 . However, if jaj is signifiGantly large, the user 

should investigate why. 

Note: To prevent excessive amounts of printout, this message will be 
printed a maximum of five times per iteration. 
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3. 8 POUT 

Operational A s suxnptiens 

• Only one major type of printed report wiU be produced per run. 

• Population reports will be produced by option for either zone, 
region or country level, 

• For Population or Monte Carlo reports any one, a set, or all 
of the parameters can be selected in one run. 

• 1%e Substrata Historical File will also be required to produce 
the first three Population reports. 

• Printed report control will be by card input. 

• All control card input data will be echo printed. 

• All control card input data will be checked for errors before 
any error will cause the processing to terminate in the middle 
of a cuse. 

• J^ternal print units will be in English, an Optional override is 
available for metric units via the control card input, 

• The Header Card and the four Data Cards are required input. 

• All data to be processed will be obtained from disk (or tape) 
files. 

• All input data files will be checked for correct case numbers. 

3. 8. 1 Input Card Data 

3. 8, 1. 1 List of Data Quantities 

See Table 3. 8-1. 

3, 8, 1. 2 Card Formats 

"POUT” is punched in card columns 75-78 of all cards. A sequence 

number is pundied in card columns 79-80. See Table 3. 8^2. 


■ J , -I V A'.. ■ H -If’ 


Table 3, 8-1.' POUT Input Data Description 


Symbol 


Nominal . 

Dimension | 'V'alue | Range 


Units 


Description 


IHEADR 


sRPTYPE 


Blanks 


gAUNITS 


72 character case header which prints out at 
the top of every page 

Major type of report identification 

01 - Substrata Reference Data Report 

02 - Population Mean, Standard Deviation and 

Histogram Report 

03 - Histograms of Monte Carlo Statistics 

Report 

04 - Frequency of Samuple Segment Acquisitions 

Reports 

External print units flag 
0 - English units 

Wheat area value in 10, 000 acres 
Production value in lOO, 000 bushels 
yield value in bushels/aere 
^0 - Metric units 

Wheat area value in 1000 hectares 
Production value in lOOO metric tons 
Yield value in quintals /hectare 


pSTART 


ItNTVl,! 


?fNTVB2 

I ^ ^ 

fBREAK2 

II 


-500,0 


100.0 


- 100 . 0 


100.0 


-999.9 to 
999.9 

0. 1 to 

100. 0 

-999. 9 to 

999. 9 

0. 1-100.0 

-999.9 to 
999.9 


Histogram interval start value 


Histogram interval value in percent 


Histogram breakpoint to change interval value 


Histogram interval value in percent 
Histogram breakpoint to change interval value 
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Table 3.. 8-1, POU T Input De 


Nominal 

Symbol : Dimension I Value 


Range 


IPARMTR 


100. 0 


0 . 1 - 100 . 0 


STOP 


500. 0 


-999. 9 to 
999. 9 




Histogram interval value in percent 

His:togram interval terminal value 

A maximum of 51 range intervals is allowed 
although less may be used. Intervals may not 
overlap. 

Report parameter type option flags for 
RPTYPE = 02 or 03. =0 - no report; ??0 - 

produce report. 

PARMTR(1)V0 RPTYPE = 02 Produce 

Population Sampling Error 
Report 

RPTYPE = 03 Produce Monte 
Carlo Area Error Report 
PARMTR(2H0 RPTYPE = 02 Produce 

Population CAMS Error Report 
^0 RPTYPE = 03 Produce Monte 

Carlo Production Error Report 
PARMTRf3i).#0 RPTYPE = 02 Produce 

Population Yield Error Report 
^0 RPTYPE = 03 Produce Monte 

Carlo Yield Error Report 
PARMTR(4)i«0 RPTYPE = 02 Produce 

Population Area Error Report 
#0 RPTYPE = 03 Produce Con- 
fidence Level Report 
PARMTR(5)^fb RPTYPE = 02 Produce 

Population Production Error 
Report 

U RPTYPE = 02 or 03 and all PARMTR value 
are zero, then all reports of that RPTYPE wi 
be produced. 
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Table 3, 8-1. POUT Inipat Data Deseription (coat’d) 


Nominal 



Unite 


Description 


CASIN 


fiEVEL 


F ■ 

r . 


0. 


0-3 


0-9999 


Si. . • 


Case aamber identification associated with tiie 
required input file. A case number must be 
input for every report produced. Relationship 
between RPTYPE, PARMTR and ICASIN 
follows. 

ICASiN(l) = case no. RPTYPE = 01 

= 02 PARMTR(I|/0 
= 03 PARMTR(1)^0 
= 04 

IGA^(2) = case no. RPTYPE = 02 PARMTR(2)^0 

= 03 PARMTR(2)70 

iCASlN(3) = case ho. RPTYPE = 02 PARMTR(3)/0 

= 03 PARMTR(3)/0 
IGASIN(4) = case no. RPTYPE = 02 PARMTR(4)^0 

= 03 PARMTR(4)jJ0 
ICAS1N(5) = case no, RPTYPE = 02 PARMTR(5)^0 


If all PARMTR values are zero for a RPTYPE = 
02 or 03, then all the ICASIN values (case 
number) niust be entered for that RPTYPE value. 


Parameter report level indicator for 
RPTYPE = 2, =1 - reports produced at zone 

level; =2 - at region level; =3 - at country 
level. • 


Case number identification for the SUBHST 
input file. Required only when RPTYPE = 2 
and PARMTR = 1, 2 or 3. 


p> 
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os 


00 
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Table 3, 8-1, POUT Input Data Description (cont'd]) 


Name 


Symbol 


Dimension 


Nominal 


Value 


Sange 


Unite 


Description 


BIOWD 


0 , 1 


f 


I 


WPRTY 


0-4 


IPRD 


3. 14 


year: >64 
month: 
01-12 
day: 01-31 




'To 


^ ' 


for 


Prediction bio-window flags: 

BIOWD(n) =1 to process bio-window n 
=0 otherwise 
RPTYPE = 02 and 

PARMTS.(2)?fO Selects by bio-window in 
record. 

Uses last yield date in 
record. 

Selects by bio -window in 
record. 

Selects by bio-window in 
record. 


PARMTR(3)?^0 


or 

for 


PARMTR(4)?^0\ 
PARMTR(5)i«0/ 
RPTYPE = 03 


Bio -window priorities: 

Lrist of bio -windows in decreasing order of 
priority, e. g. , 4, 1, 3, 2 or 3, 1, 0, 0, Used only 
by RPTYPE = 02 and PARMTR{2)^0 in conjunc- 
tion with the prediction dates below. 


Prediction dates (up to 14 dates) for the selection 
of data for RPTYPE = 02 or 03. 

IPRD{1, n) = year - 1900 k 

iPRD;(2, n) - month S 

IPRD(3, n) = day ^ 

I 

The dates must be in ascending order. The firsg 
zero date terminates the list. n 

For RPTYPE = 02 and . V' 


PARMTR(2)i^0 Prediction date selected as ^ 
a function of priority above | 


and latest date less than or S 


PARMTR(3)^0 


PARMTR(4) or 
and all RPTYPE = 03 






.equal to the input prediction 
date. . 

Selects the latest yield date^ 
less than or equal to the ® 

input prediction date. ^ 

the input prediction 
dates must match with 
file dates; otb.erwis.e a 
niessage is printed 
and that date nkinned. 
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Table 3, 8-2.. POUT Data Card Formats 
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Table 3, 8-2, POU T Data Card Formats (cont'd) 
DATA CARD 3 





75 


79 



DA TA CARD 4 


l^iv 


.-VJ 

|:J 


1 

3 

5 

8 

10 

12 


15 

17 

19 



1 









F? 
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GCj 



i s? 
Q: 
pi 


1 


p 

pi 
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3, 8. 1. 3 Deck Se tap 

1. Header Card - sequence 01 

2. Data Card 1 - sequence 02 

3. Data Card 2 • sequence 03 

4. Data Card 3 - sequence 04 

5. Data Ca-rd 4 - sequence 05 


3. 8. 1; 4 Rules for Entering Data on Cards 
3.8. 1.4. 1 General 

1. Integers must be entered right- justified. 

2. Alphanumeric names must be entered left- justified. 

3. F format numbers must have the decimal point present* i. e. > 
F6*l iXXX.X 

4. The card sequence numbers in CC. 79-80 must be present on 
all cards. 

3. 8. 1.4. 2 Specific Fields 

The correspondence between the RPTYPE, PARMTR and ICAS£M 
input values, and the required input files versus the produced reports is 
shown in Table 3.8-3 below. 

Histogram ranges* PARMTR, BIOWD, WPRTY, IPPD input values 
are' required only if RPTYPE=02 and ^03. LEVEL is required only if 
RPTYPE-02. 




28234^6025-RU-00 
Page 201 


Table 3, 8-3, Input Flag Value Definitions 


RPTYPE ! where I = 


where I = 



Req. Input 
Files 

Report 

SUBHST 

Substrata Reference Data 

SEGTRU 

Population Sampling Error 

CAMSF 

Population CAMS Error ■ 

YESOUT 

Population Yield Error 

CASF 

Population Area Error. 

CASF 

Population Production Error 

CASDIS 

Monte Carlo Area Error 

CASDIS 

Monte Carlo Production Error 

CASDIS 

Monte Carlo Yield Error 

CASDIS 

Confidence Level 

ACQUIS 

Frequency of Sample Segment 
Acquisition 


SSUBHST file also required for each report. 


41 














3. 8. 2 Card Inputs 


Substrata Reference Data Report 


TEST GASt i o HPTYPF.SI,. SUOSTRaTA HEfERENCt OAtA REPORT COnTVAC HAND OFF LPP SlMUUATlON PAGE I 


KrTYHK 

PaRMT 

igastn 




lEVEL 

iCSESH 

1 

lone ! 

>: I 

1 

-fl 

-u 

-0 1 

1 


tNTPRyAL STAR! TNTWH BREARt 1 NTvL 2 BRFAKa INTVL3 


0* 0*0*0 •o*o*o*o 


0 * 0-0 *o-o-o‘ - 0 - 0*0 -o**o-u - 0 -n-o -o-o-o 


Population Histogram Report 

TfcST CASt 2 - HPTYKFar, PAKMTR(3)=fl PAKWTRC5)oi YIELD AND PROD REPORTS LPP SlhOLATlON PAiiE 


KPTVPr^ 

parb r 

TEAS IN 



LEVEL 

TCSESH 


0 0 10 1 

I 

1 1 

i- 1 

J 1 

1 


interval start 


BHEaKI INTVL2 BR€AK2 INTVL3 


1 I 1 0' 3 1 ^ a 


76 U5 
•*0— Oi— 0 


6 

0-:0«h0 


• 0 - 0 * 0 *0 - 0 - 0 -0 - 0- 0 
















mmmm 


3 . 8 » Z Sajaaple Cftrd Inpiats (cont' d) 


Monfee Carlo Histograxn Report 


UST cask 3 - RPT YPFS3 * P AKMTR ( t . 3) ■! 


MUNTE CAR(,0 reports 


RPTYiPF PARrtT ICasIN LEVEL ICSESH 

3 I 0 I 0 0 1 1 1 1 1 1 1 


INTERVAL 5TARI TNTVLl RREaKI liNTVti BREAK2 XNTVL3 STOP 

~ SO..O lO.U yln.itf .5 10.0 10.0 SO.O 


BTOWQ PPPTY 
"i i 10 3 1 i» 0 


IPHD 

76 U5 76 2 1 76 2i5 

pQ^OwQ 


<"0^U*0 »0*0»t> <i*0"0«0 

*^*?Q*0! iPpyOyO. 


Acquiaitioxi Data Report 


G 

2 6 
E5 


TEST CASE 1 » ftiPTYPEs« ACQUXsTlON DATA REPORT CUNT VAC HAND OFF) 


Kf’TYPE PARmT TeAsTN LEVEL ICSESH 

Ji .o-O-O-O-O ’ i "O «0 «0 -0 *0 *0 




INTrR'VAL start TmTVLI HHFAK i INTVL2 . BHFAK2 ^IlNTyLS^^l'^ 
-.,0„ »;« iU -« .U •‘t.O ^ ZjJL 


oa — 


m 


Bl ut u. fjWI V 

• 0- U ♦ 0« 0 - u ^ 0 -- 0- 0 


iPhiiJ 

»o«o»o 


‘0»o^y 


•O-CI-iO 

*0*-0»iQ. 


•O' 0-0 


>0-0-0 

>0*0*0 


*U— 0-— 0 *0>*0— 0 

*0#0*0 -fl*0»0 




UPP SIHULATIUN PAUE 1 


LPP SIMUlATIUN PACE 1 


oo-nuj-szo^-msz 
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3, 8, 3 Output Report Data Def initio ns 


See Table 3. 8-4. 

3. 8. 4 Sample Pages from Ofitput Repoi’ts 

1. See Figure 3. 8-1 for Sample Substrata Reference Data Report. 

2. See Figure 3. 8-2 for Histogram Output Report. Note that all 
the five Population Reports and the four Monte Carlo Reports 
use this same format. Only the title changes. 

3. See Figure 3. 8-3 for Frequency of Sample Segment Acquisitions 
Report* 

3. 8. 5 File Reguirements. 

Input: AU files are optional and files required are specified via 

the report options chosen. 


SUBHST 

SEGTRD 

CAMSF 

YESOUT 

CASF 

CASDIS 

ACQUIS 


Substrata Historical Data 
Seg;ment Truth Data 
CAMS Output File 
YES Output File 
CAS Output File 
CAS Distribution File 
Acquisition Data File 


Output: None 


3. 8. 6 Error and Reeoyery 
3. 8, 6. 1 Input Errors 


, AH input errors in t-Ms program will be fatal. 

1. Label and sequence ehecking on control cards. Message 

*** IMPROPER LABEL AND SEQUENCE NUMBER ON 

POUT CONTROL CARD NO. LABEL AND SEQ. NO 

2. Range testing on RPTYPE value. Message 

RPTYPE VALUE OUT OF RANGE. RPTYPE = 
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= su .n; 


'V n?; 







■4 

fei 


■f . . 1 


Kn'' . V 

?&■'-' .a 


I 




f' 

V,., 


Table 3. 8-4. POUT Outpat Report Bata Definitions 


' Repor t Name 

Symbol . 

Range 

Units 

Report ; 

Description 

Country 

- 

- 

- 

1# t 2, f 3. 

4 character country ID 

Case Number 

- 

1 - 9999 


1. , 2. , 3. 

Case ID no. 

Strata 


1 - 9999 

- 

1..2 

Strata no. within zone 

Zone 

- 

1 - 999 

- 

1. , 2. , 3. 

Zone no. within region 

Region 

; ■ 

: 

1 - 999 

- 

1. , 2. , 3. 

Region no. within country 

Hist WA 

^ WA 

0 - 999999 

10000 acres 

1. 

Historical wheat area at various ID levels 




or 1000 
Hectares 



True WA 

WA 

0 - 999999 

10000 acres 

1. 

Tree historical wheat area at various ID 



or 1000 


levels 




Hectares 



No. in substrata group 
I 2 3 

- 

“ 


1. 

The number of substrata belonging to Group 1, 
2, and 3, respectively 

No. of segments of group 




1. 

The number of segments belonging to Group 1 

1 2 “ 3 





and Group 2 substrata 

Input bias % Ms,te 

- 

0 - 9999 

% 

1. 

PCT deviation between the true WA and the 
historical value computed from input bias data 

Simulatipn bias % hiB-ts 

- 

0 - 9999 

i % 

1. 

PCT deviation between the true WA and the 
historical value used in the simulation 

Input CY % true 

- 

0 - 9999 

% 

1. 

i Aggregated within county area std. deviation 
in PCT of true WA 

Simulation CV % true 

- 

0 - 9999 

■ 

1 

% 

1. 

Std. deviation of difference between substrata 
true WA and historical w‘A in % of true area 






t) ~ ^ 

m 00 : 

■ • ■ ■ ' 

■4/^y . ■■-i* ,-- i 






'S' . 1 } . '■ 




ts) 1 

. o o>- 






ui o 

tNi). 

UI* 



% ■' 


1 

1’ 

» i 
a 

, i * ' 


: ^ ■ -| 



1 

1 

o 


1 

: 

9 

• 


p 















, POUT Output Report Data Definitions (cont'd) 


Unite 


Report 


2. , 3. 
2 . 


Description 


Name of input file used to generate report 

A table illustrating' for any parameter the 
ramge of values that would have to appear 
in any 1 of 51 slots or buckets - histogram 


Name of the variable for which the histo- 
gram is being printed 

Mean value of the data readings of the 
variable under study 

Standard deviation of the data readings of 
the variable tmder study. 

Value used to normalize the range of the 
reported variable 

Expected value of the population sum 
Standard deviation of population sum 


Minimum data reading of variable 


Maximum data reading of variable 


No. of data readings of variable 
Month, day and year of prediction point 


No. of times that a d.-t:: -"eading falls in a 
; slot divided by the total no. of readings. ® 

: 51 slots and cum percentage from the first w 
bucket to the bucket of interest t 

Total no. of segments that were tallied 9 

within an ID level 


;% of segments that have at least one acquired 
observaiion in a particular bio window combina 
tion 










TfST 1 • HPiVPf:-l9 aUHJVrHATA RerpKENCE DATA Ri;PORI (UNIVAC HAND OfP LPP SIMULATTPH PADE B 


i I 


% 


t'.' 

■ 

K; ■- 


■I iV 


I ■" . 



* 


S U H S 1 K A 

T A i 

« e F 

E RE N C 

E l> A T A 




gounvhv 

USA 

CASE MU Hum 

1 











A f 

1 € A 




• •■ ■ - ■■ 

■...,. ...... 






(THOOSaMD he CT awe 8) 




H 

■ 







Hf). IN 


NO. 

OF 


i 







5l/ft^STkA1 A 


SfeGHEwTS 

INPUT 

simulation 

INPUr SlHOLATlsON 




MIST 

TliUt 

afTuUP 


OE GROUP 

blAS 

BIAS 

CV2 IVI 



STRATA 

HA 

ha 

1 2 

3 

1 

2 

FCT HIST 

fgt hist 

FCT true FCT TRUfc 



- 20 

Jf|‘0 « • 7 


0 5 

0 

0 

0 

.05(|!Q 

.0357 

',.084:9 ,0420 



50 

Tbl.9 

16S.1 

2 0 

C 

2 

0 

.0500 

♦ 0202 

♦0707 P0099 



70 


5 u n . 3 

1 0 

0 

2 

0 

. 05 OiO 

.2396 

,100>0: ,0000 



9(i 

«as.6 

59P.7 

0 1, 

5 

0 

0 

♦d5a 0 

,1175 

.08 0:9 .0268 



* * 4 ♦ ^4 4 ♦ ♦ ♦ +4 39 » % « *4 ^ * 4 * ♦ 444 4 * ♦ ♦♦44 4 4444 4 444 4 4 444 444 4444 4 44 4 444 44 4 44 4 444 44444 44 4 4 ♦ ♦444444 



70ftt 









,05U .0275 



ii 

1133.1 

T^^7.^ 

5 6 

S 

■ 4 

0 

,05di0 

..0630 

‘ 


♦ ^ ♦44i4 ♦44 4 ♦ ♦♦ 4444 444 444444 4444 444 444 4444444444444444444444444^ ♦4444A 44 44 4 444444 M 444 4444444444V ♦444444* 



?0 

BO.*? 

/3.0 

0 0 

J 

0 

0 

/ pU5 y f> 

-.0075 

, tOOiO ,0000 



SO 

60.9 

96.0 

0 0 

1 

0 

0 

,0500 

.lB6o 

♦iooo ,oood 



f^O 

a 0.9 

BH,:P 

0 0 

1 

0 

0 

♦ 0 5.00 

♦ 0976 

,10 00 .OiUOO 



♦ ♦ ;l ^ ♦ ♦♦ ♦ ♦^ ^♦^♦♦4 ♦ ♦♦ ♦ *:4^.4* v^44 ♦44444 ♦♦4444 444 4*444 44 444 ♦ 444 44444444 4444444444 4*4444 *44444 444 4 4444 44* 



z.a^it 

10 


?’j7.'J 

0 0 

3 

0 

0 

.0500 

♦ 0622 

,1000 ,0000 



♦■♦44444444V44 •♦4 144 ♦*♦*4444444* ♦4 4i*4*444*44444* ♦4444*444 * 44 V *4 **4444 **♦**♦ *44!44*4*44*4 44*44444 444 44 V*#** f 



f 0 

ttiO 9 . (J 

P-itv.O 

1 6 

3 

1 

2 

.05l):0 

• 0329 

,0317 .0291 




7 aw , /4 


9 0 

0 

18 

0 

• 0500 


♦0337 ,0510 



50 

■’ Tv«7,^ 

70 3.2 

fl 0 

0 

18 

0 

.o5os6 

,0860 . 

♦035d o0257 



SO 

To 

eo9,4 

fi/7,B 

/I b 

0 

6 

1 

♦ 0500 

,0805 

♦0317 ♦0169 



«0').7 


3 1 

1 

3 

0 

Q O SOif) 

*09 39 

.0849 .0871 



RO 


597-0 

5 2 

0 

6 

d 

,0500 

-♦0336 

, 03/8 , 0217 ; 

&} 00 

crq 

— Sg.'-OI 


90 

566,6 

S96.6 

7 0 

0 

7 

0 

♦ 0500 

,0530 

,0379 .0326 



Figure 3, 8-1* Sample Substrata Reference Data Report 
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FigaTe 3, 8-2. Histogram Output Report (cont'd) 
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Figure 3. 8-2, Histogram Output Report (cont* 
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If ST CASE 3 - RPTYPES39 
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Figure 3. 8-2. Histogram Output Report (cont'd) 
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Figure 3, 8-3. Frequency of Sample Segment Acquisitions Report 
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Range testing on Histogram interval values 
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Messages 


*** START VALUE GREATER THAN BREAK I VALUE. 
START = BREAKl = 


*** BREAKl VALUE GREATER THAN BREAK2 VALUE. 
BREAKl = BREAK2 = 


BREAK2 VALUE GREATER THAN STOP VALUE. 
BREAK2 = STOP = 


*** THE SEGMENT BETWEEN START AND BREAKl WILL 
NOT DIVIDE INTO EQUAL UMTVLl INTERVALS. 
BREAK 1 - START = . INTVLl = 


*** THE SEGMENT BETWEEN BREAKl AND BREAK2 WILL 
NOT divide INTO equal INTVL2 INTERVALS, 

BREAK2 - BREAKl = . INTVL2 = . 


*** THE SEGMENT BETWEEN BREAK2 AND STOP WILL 

NOT Divide into Equal intvls intervals 

STOP - BREAK2 = ’ . INTVL3 = 


REQUESTED NO. OF HISTOGRAM INTERVALS EXCEEDS 
MAXIMUM OF 51. 

REQUESTED NO. OF INTERVALS THAT WOULD BE 
GENERATED = 


4. Case number value missing from input. Message 


*** ICASIN(i) VALUE MISSING FOR REPORT TYPE. 

ICASIN(i) = RPTYPE = PARMTR(i) = 


3, 8. 6.2 ProGessing Errors 


1. «««>}= INPUT PREDICTION DATE DOES NOT MATCH DATE 

ON FILE. 

RPTYPE = PARMTR{i) = . INPUT DATE = 


Message in non -fatal. The input date is skipped. 


2. Case number on control card does not match case number on 
file. Message 


*** ICASIN(i) DOES NOT MATCH CASE NUMBER ON FILE. 
ICASINCi) - FILE CASE NO. = ■ . FILE NAME = 


The report associated with the missing file is skipped and the 
next report tried. 






1 
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4. 0 NOTES ON PROGRAM USE 


4.1 GENERAL 


EPHEMS will only have to be run about once a year once a particular 
orbit has been modeled correetly. 

Th^ utUity, GRID, to generate the INDEX matrix should only-need,to be ' 
run once. 

SEE will require rerunning any time an error model Will change. It 
would be expected, however, that it would be desirable to have several 
different sets of models available. Consequently, a group of runs may be 
needed but once these runs are completed, another run should not be needed 
for a long time. Since the files generated by SEE are small, it may be 
desirable to keep them on disk as cataloged files. 

LUMP will have to be run for each country and each time the data 
base is changed. The three output files from each LUMP run will have to 
be kept on tape. 

SAGE will only be run normally when LUMP has been run so that the 
SEGREF file should be kept on tape. Since this is the only use of SWATId 
and SWATHR files, they should be kept on tape also, 

SACS can be run quite frequently if it is desired to change the coe.- 
ditions of segment access. SAGS is a fast running program and the ACQUIS 
file is small compared to SEGREF. Whether to rerun SAGS or try to use 
an existing ACQUIS file from a library depends on how many file s have to be 
saved and accounted for. 

LEM is the main LACIE simulation application program. Thus it will 
receive the most usage, if Monte Carlo iterations are not being used, the 
runs will be reasonably fast and inexpensive. Additional time can be saved, 
if desired, by using saved SEGTRU, CAMSF and YESQUT files where possible. 
When running in this mode, consideration must be given to the number of 
saved tape files that can accumulate vs. the run time required to regenerate 
the files. GASF and GASDIS are the only required output files. The pro- 
gram will run very slow if all subprograms are executed for more than just 
a few iterations. « 






28234- 6025 -RU -00 
Pa ge 217 


POUT will usually be run when more detailed statistical analysis is 
required of LEM output data. This program generates no output files and, 
for any one output report, will run fast. 


4. 2 RESTART 

No provision for restarting any program in the middle of a run due 
to an error has been made. The frequency of use and the run times of all 
programs except for LEM make this feature unnecessary. Restart is 
provided in LEM between each Monte Carlo iteration. However, the way 
the program is implemented, restart can be used only in the non-errpr 
mode. If EXEC VUI as used at NASA has a user error wrapup exit, 
then the wrapup code from LEM could be duplicated in an error routine. 

H this routine is called in case of a processing fatal error, then thO last 
random number seed could be printed and the CASDIS and CASE files can 
be closed, 

It should be noted that since extensive input data error checking is 
performed at the beginning of all programs, the only job terminating error 
that is apt to occur is a system /hardware error. 

Noter On restart it is necessary to input the final active random 
number seeds (RSEED(i)) which are printed out at the end of the previous 
LEM run. Care must also be taken to enter the appropriate values of 
RSTART and NTRIAL. For example, if a restart is desired after five 
trials to continue for another five trials 

RSTART - 5 (number of previous interations) 

NTRIALS = 10 (total number of iterations) 

4.3 USAGE OF LEM 

The LEM program contains the Segment Truth Generator, CAMS, 
YES, and CAS subprograms, and provides the Monte Carlo Gontroi for 
the execuidon of these modules. The LEM control coupled with the GAMS 
options allow the user to tayior error model configurations to his specific 
needs. The options available ie the user have been previously defined. 

This section will illustrate the use of these options to cons tr act meaningful 
cases by use of specific examples. 
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Figure 4. 3*1 shows the LFM control data flag set up for ten specific 

cases running from no errors to all error sources with various combina- 

tions. These are but a few of the possible configurations that the user may 

construct. Case 1 gives the sample error oniitting within county variance 

or cloud cover effects. Cases 2 and 3 add cloud cover, effects and Within 

county Vg riance, respectively^ Cases 4 and 5 examine YES and CAMS 

errors, respectively. Note that there are three different CAMS errors 

that may be controlled separately or in combination with each other. 

In addition to this control, t^e specific form of the errors can be con* 

trolled by the CAMS input options. Case 7 examihes the. effect of 

varying the multi-year historiGal statistics used for GAS aggregation 

over a number of iterations P'1). Cases 8 and 9 give examples of hold- 

s fc 

ing data constant from previous iterations (1 ) or from previously 

generated files. The errors or combination of errors so treated may 
be varied in accordanGe with the user's desires, in the latter ease, care 
must be taken to mount the p:roper files to be used. Case 20 is a Monte 
Carlo case. This example shows only the classification error varying; 
however, any error or combination of errors may be varied. The restart 
option allows the user to run a subset of the Monte Carlo trials, examine 
Ihe results, and continue with the run at a later time. 

The CAMS subprogram provides options for different classification 
models, sigriature extension models and the use of multi-temporal sam^ 
pling. The classification model flag, JMODEIj, maybe set to; 

1 =1 - Complex omission/ commission model 
= 2 simple error model 

Care must be taken to ihsure that the proper error data is defined in SE^ 
for the model chosen. The signature extension modei for ordinary seg- 
ments includes additive and multiplicative bias and randonU errors. The 
option of selecting the additive or multiplicative case for the random 
component can be exercised by setting; 

- additive random component 
- multiplicative component 


ISIGEX 


= 0 
3^0 




CASE DESCRIPTION 

1. No errors, no acq, data (samiple error) 

2. No error s, use acq. data (cloud cover elEects) 

3. Seigment truth variation (within county variancei) 

4. YES error only 

5. CAMS errors only 

6. All error sourees 

7. Vary CAS multi-year effect 

S*t 

8. Mold I iteration CAMS error constant and vary 
YES error 


9» Read CAMS errors from e3d.stin>g file and vary 
YES error 


10. Monte Carlo classification error 




/ ^ LEM CONTROL DATA 
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1 2 


3 1 

1 1 

3 1 

1 0 

3 1 

1 0 

0 1 

1 0 

3 1 

1 0 

0 1 

1 0 

0 

0 0 

0 1 

1 0 

0 1 

1 0 

1 1 
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Figure 4. 3^1. EEM Control Data Examples 
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The multi-temporal sampling effect Is provided to simulate the effect of 
segment acquisitions in more than one window. The simulation is accom- 


plished by multiplying the error of the present window by a factor which 
depends on the previous windows acquired. The option for using this 
feature may be exercised by setting: 


IMXJLTI 


= 0 - include multi -temporal sampling effect 
^ 0 - bypass multi -temporal sampllhg effect 


4.4 USAGE OF POUT 

The preparation of control data to create the substrata reference 
data report or the frequency of sample segment acquisition report is 
straightforward and will not be discussed here. See Sections 3. 8, 1 and 
3. 8. 2 for details. The preparation of control data for the population and 
Monte Carlo reports is more Gpmplicated and examples are discussed below. 


4, 4. 1 Poputetion. Report Example 


Assume that the' population CAMS error report and the population 
area report is to be geaerated. English units are to be used. Two 
intervals are to be used for the histograms. iO between 0 and. 200 and 25 
between 200 and 475. The reports are to be produced at the region level. 
Beside s the three bio -windows, there are four prediction date s. The files 
SUBMST, GAMSF and CASF must be mounted. The data inputs would be 
as follows: 


RPTYPE = 2 
AUNITS = 0 


START = 0 > 


IWTVL,! = 10 


BREAK 1 = 200 
INTVL2 = 25 




Dp not enter data for INTVIjS or STOP in this case 


BREAK2 = 475 J 

PARMTR(2) - 3, PARAMTR(4) = 6 
ICASIN(2) = 1001, ICASiN(4) = 2005 


LEVEE 2 


BIOWD(l) = 1, BIOWD(2) == 4, BIOWD{4) = 2 
WPRTY(l) = 2, WPRTY(2) = 4, WPRTY(3) = I 

IPRD(1, 1) = 71. 7. 20, IPRD(1,2) = 71, 9. 15, IPRD(1, 3) = 71, 11, 21, 
IPRD(i,4) = 72,2, 1 
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4* 4. 2 Monte C?arlo Repoyt Option 

As some that the Monte Carlo production error report is to be 
produced. Metric units are to be used. The nominal irterval and start 
and stop values are to be used for the histograms. The report is to be 
produced at the zone level. Two bio-windows are to be used and two 
prediction points. The only file required is the CASDIS file. The control 
card data would be as follows: 

RPTYPE = 4 

AUNITS = 2 

PARMTR(2) = 7 

1CASIN(2) = 7021 

BIOWD(l) = 5, BIOWD(2) = i 

IPRD{1, 1) = 74, 4, 17, 1PRD(1, 2) = 74, 7, 1 


